Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,339,336 T→C 5.0% S3P (TCG→CCG)  pgpB → phosphatidylglycerophosphatase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,339,3360TC 5.0% 246.5 / 3.1 81S3P (TCG→CCG) pgpBphosphatidylglycerophosphatase B
Reads supporting (aligned to +/- strand):  ref base T (37/38);  new base C (2/2);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

TGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCT  >  NC_000913/1339191‑1339483
                                                                                                                                                 |                                                                                                                                                   
tGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGa                                                                                                                                                  >  8:111131/1‑149 (MQ=255)
 ggCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGa                                                                                                                                                  <  5:212001/148‑1 (MQ=255)
  gCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGAtt                                                                                                                                                >  1:135720/1‑149 (MQ=255)
  gCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGAtt                                                                                                                                                >  8:108945/1‑149 (MQ=255)
  gCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGAtt                                                                                                                                                <  2:53642/149‑1 (MQ=255)
  gCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGAtt                                                                                                                                                <  2:174007/149‑1 (MQ=255)
  gCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGAt                                                                                                                                                 <  2:42934/148‑1 (MQ=255)
  gCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGCAGAGGCAATGCGTTCGAtt                                                                                                                                                >  8:353753/1‑149 (MQ=255)
       ttttGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCa                                                                                                                                            <  4:108694/148‑1 (MQ=255)
                   gagaAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGt                                                                                                                               >  6:355227/1‑149 (MQ=255)
                   gagaAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGt                                                                                                                               >  4:158065/1‑149 (MQ=255)
                    agaAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAACTCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTg                                                                                                                              <  7:124850/149‑1 (MQ=255)
                     gaAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTgg                                                                                                                             <  8:54850/149‑1 (MQ=255)
                       aGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGa                                                                                                                           <  7:79296/149‑1 (MQ=255)
                           gTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTg                                                                                                                       <  7:108945/149‑1 (MQ=255)
                                    aTCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTAtt                                                                                                                <  2:135720/147‑1 (MQ=255)
                                    aTCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGc                                                                                                              <  5:63971/149‑1 (MQ=255)
                                     tCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTa                                                                                                                                      <  7:326214/124‑1 (MQ=255)
                                     tCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTa                                                                                                                                      >  8:326214/1‑124 (MQ=255)
                                     tCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCt                                                                                                             >  4:220491/1‑149 (MQ=255)
                                      cACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCtt                                                                                                            >  4:286252/1‑149 (MQ=255)
                                       acgacgCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGa                                                                                                                                                  <  7:293125/110‑1 (MQ=255)
                                       acgacgCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGa                                                                                                                                                  >  8:293125/1‑110 (MQ=255)
                                         gacgCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAACAGGAGAGGCCATGCGTGCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTc                                                                                                         <  5:165057/149‑1 (MQ=255)
                                         gacgCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTc                                                                                                         <  3:321206/149‑1 (MQ=255)
                                         gacgCGTCGATCGTTCACGGAATGACACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTc                                                                                                         <  2:151087/149‑1 (MQ=255)
                                          acgCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCa                                                                                                        <  3:340399/149‑1 (MQ=255)
                                          acgCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCa                                                                                                        <  5:90064/149‑1 (MQ=255)
                                            gcgTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCa                                                                                                        <  3:46745/147‑1 (MQ=255)
                                            gcgTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCa                                                                                                        >  4:46745/1‑147 (MQ=255)
                                              gTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGTGACTTATCAAAAAGGAGAGGCCATGCGTCCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGc                                                                                                     >  3:188400/1‑148 (MQ=255)
                                              gTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGTGACTTATCAAAAAGGAGAGGCCATGCGTCCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGc                                                                                                     >  3:193028/1‑148 (MQ=255)
                                              gTCGATCGTTAACGGAATCGCACAAAACTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTTTTATTGTGGGCCTTAACAAAAAGGAGAGGCCAAGCGTACGATTGCCAGACGTACCGCAGTGCGAGCTGCACTATTGCTTGTCATGc                                                                                                     >  8:330034/1‑148 (MQ=255)
                                                   tCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAg                                                                                               >  4:279234/1‑149 (MQ=255)
                                                      ttCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGAATGCCAGACGTACAGCAGTGGGAGCTGCCCTATTGCTTGTCATGCCAGTAGCCg                                                                                            >  4:259523/1‑149 (MQ=255)
                                                         aCGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTAt                                                                                         >  4:41271/1‑149 (MQ=255)
                                                         aCGGAAGGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTAt                                                                                         <  6:211806/149‑1 (MQ=255)
                                                            gAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGa                                                                                      <  1:32881/149‑1 (MQ=255)
                                                             aaTGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGTGACTTATCAAAAAGGAGAGGCCATGCGTCCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGAt                                                                                     <  4:193028/149‑1 (MQ=255)
                                                             aaTGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGTGACTTATCAAAAAGGAGAGGCCATGCGTCCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGAt                                                                                     <  4:188400/149‑1 (MQ=255)
                                                               tGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGAttt                                                                                   >  3:27632/1‑149 (MQ=255)
                                                                   acaAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGc                                                                                                                      <  8:256884/110‑1 (MQ=255)
                                                                   acaAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGc                                                                                                                      >  7:256884/1‑110 (MQ=255)
                                                                   acaAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTctgg                                                                               <  7:330178/149‑1 (MQ=255)
                                                                            tCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTg                                                                      <  7:133356/149‑1 (MQ=255)
                                                                            tCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTg                                                                      <  7:230180/149‑1 (MQ=255)
                                                                                    tCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCt                                                               <  8:68962/148‑1 (MQ=255)
                                                                                      cTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGga                                                            <  4:284932/149‑1 (MQ=255)
                                                                                       tGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGAGTTCTGGCTGGCGTTGGCAACCTg                                                              >  3:180527/1‑146 (MQ=255)
                                                                                         gtaCTGGCGGTTTCAGGTTTGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGaa                                                         <  7:353753/147‑1 (MQ=255)
                                                                                         gAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAATGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGaa                                                         >  6:179893/1‑149 (MQ=255)
                                                                                         gAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGaa                                                         >  6:97192/1‑149 (MQ=255)
                                                                                              ggCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGCGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAg                                                    >  6:21854/1‑149 (MQ=255)
                                                                                               gCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTgg                                                                                                                             <  3:306307/75‑1 (MQ=255)
                                                                                               gCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTgg                                                                                                                             >  4:306307/1‑75 (MQ=255)
                                                                                               gCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGt                                                   <  4:27632/149‑1 (MQ=255)
                                                                                                cGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGtt                                                  <  3:158065/149‑1 (MQ=255)
                                                                                                cGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGtt                                                  <  3:41271/149‑1 (MQ=255)
                                                                                                     tCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGcta                                             >  5:274837/1‑149 (MQ=255)
                                                                                                         ttATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGctacta                                          >  3:112998/1‑148 (MQ=255)
                                                                                                          tATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTaa                                         >  4:17318/1‑148 (MQ=255)
                                                                                                           aTGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTaa                                         >  6:239075/1‑147 (MQ=255)
                                                                                                                     gACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCttttt                             >  2:137739/1‑149 (MQ=255)
                                                                                                                     gACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCttttt                             <  5:239075/149‑1 (MQ=255)
                                                                                                                       cTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTgg                           <  7:111131/149‑1 (MQ=255)
                                                                                                                       cTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTgg                           <  3:17318/149‑1 (MQ=255)
                                                                                                                        ttATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTggg                          >  4:207854/1‑149 (MQ=255)
                                                                                                                          aTCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGt                         <  3:286252/148‑1 (MQ=255)
                                                                                                                             aataaGGAGCGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTTGCCGTATGTATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGTATACTAAACGCGGCTTTTTGGGTTAc                      <  3:259523/145‑1 (MQ=255)
                                                                                                                             aaaaaGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACt                     >  8:275342/1‑149 (MQ=255)
                                                                                                                             aaaaaGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACt                     >  1:172300/1‑149 (MQ=255)
                                                                                                                             aaaaaGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACt                     >  8:133/1‑149 (MQ=255)
                                                                                                                             aaaaaGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACt                     >  4:101503/1‑149 (MQ=255)
                                                                                                                                   gagaGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACt               >  3:193944/1‑149 (MQ=255)
                                                                                                                                     gagGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGt             <  5:355227/149‑1 (MQ=255)
                                                                                                                                      agGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACATGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGt             >  1:236336/1‑148 (MQ=255)
                                                                                                                                       ggCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCAAGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCa           >  6:143070/1‑149 (MQ=255)
                                                                                                                                           aTGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCa           <  2:343522/145‑1 (MQ=255)
                                                                                                                                           aTGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCa           >  1:343522/1‑145 (MQ=255)
                                                                                                                                            tGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAg      >  8:296016/1‑149 (MQ=255)
                                                                                                                                            tGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAg      >  5:78555/1‑149 (MQ=255)
                                                                                                                                                ttCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCt  >  7:33640/1‑149 (MQ=255)
                                                                                                                                                ttCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCt  >  5:14808/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCT  >  NC_000913/1339191‑1339483

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑