Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,681,446 G→T 6.3% A157A (GCG→GCT folM → dihydromonapterin reductase, NADPH‑dependent; dihydrofolate reductase isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,681,4460GT 6.3% 233.3 / 2.7 63A157A (GCG→GCTfolMdihydromonapterin reductase, NADPH‑dependent; dihydrofolate reductase isozyme
Reads supporting (aligned to +/- strand):  ref base G (25/34);  new base T (2/2);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

ATACCCCATACCTGCTCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCA  >  NC_000913/1681307‑1681581
                                                                                                                                           |                                                                                                                                       
aTACCCCATACCTGCTCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTCGCGACAAACATATTGCGTATGCTGCAAGCAAAGCg                                                                                                                                         >  1:1537209/1‑140 (MQ=255)
aTACCCCATACCTGCTCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCg                                                                                                                                         >  1:1537021/1‑140 (MQ=255)
   ccccATACCTGCTCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCa                                                                                                                                      <  1:21185/140‑1 (MQ=255)
   ccccATACCTGCTCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCa                                                                                                                                      <  2:1460914/140‑1 (MQ=255)
   ccccATACCTGCTCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCa                                                                                                                                      <  2:1402903/140‑1 (MQ=255)
   ccccATACCTGCTCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCa                                                                                                                                      <  1:996511/140‑1 (MQ=255)
      cATACCTGCTCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTg                                                                                                                                   >  1:519271/1‑140 (MQ=255)
       aTACCTGCTCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTgg                                                                                                                                  <  1:536298/140‑1 (MQ=255)
       aTACCTGCTCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTgg                                                                                                                                  <  1:77174/140‑1 (MQ=255)
       aTACCTGCTCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTgg                                                                                                                                  <  1:77165/140‑1 (MQ=255)
       aTACATGCTCAACCATGCGCTGGAAAGATTTCTGAGTGGGCACGGACACGCAGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAAAATATTGCGTATGCTGCAAGCAAAGCGGCACTgg                                                                                                                                  <  1:1666550/140‑1 (MQ=255)
        tACCTGCTCAACCATGCGCTGGAAAGTTTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCTCTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAAGATATTGCGGATGCTGCAAGCAAAGCGGGACgggt                                                                                                                                 >  2:1636333/1‑139 (MQ=255)
             gCTCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGataata                                                                                                                            <  1:1304174/140‑1 (MQ=255)
             gCTCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGataata                                                                                                                            <  1:969183/140‑1 (MQ=255)
              cTCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCTGCACTGGATAatat                                                                                                                           >  1:821695/1‑140 (MQ=255)
              cTCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCTGCACTGGATAatat                                                                                                                           >  1:823070/1‑140 (MQ=255)
              cTCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAatat                                                                                                                           >  1:161520/1‑140 (MQ=255)
              cTCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAatat                                                                                                                           >  2:531917/1‑140 (MQ=255)
              cTCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAatat                                                                                                                           >  2:81339/1‑140 (MQ=255)
               tCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATg                                                                                                                          <  1:1022119/140‑1 (MQ=255)
               tCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATg                                                                                                                          <  1:1021763/140‑1 (MQ=255)
               tCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATg                                                                                                                          <  1:1022152/140‑1 (MQ=255)
                     aTGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGc                                                                                                                    >  1:1100150/1‑140 (MQ=255)
                     aTGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGc                                                                                                                    >  1:624484/1‑140 (MQ=255)
                      tGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCt                                                                                                                   <  2:519271/140‑1 (MQ=255)
                        cgcTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCg                                                                                                                 >  1:176730/1‑140 (MQ=255)
                        cgcTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCg                                                                                                                 >  1:177811/1‑140 (MQ=255)
                                       tGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTg                                                                                                  >  2:66024/1‑140 (MQ=255)
                                       tGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTg                                                                                                  >  2:190695/1‑140 (MQ=255)
                                       tGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTg                                                                                                  >  2:1643021/1‑140 (MQ=255)
                                       tGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATAGTGCTTATGCTGCAAGCAAAGCGGCACTGGCTAATATGACCCGCTCGTTTGCCCGCAAGCTg                                                                                                  >  2:961489/1‑140 (MQ=255)
                                            gggCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCAcc                                                                                             >  2:365053/1‑140 (MQ=255)
                                            gggCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCAcc                                                                                             >  2:365038/1‑140 (MQ=255)
                                            gggCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCAcc                                                                                             <  1:1174543/140‑1 (MQ=255)
                                              gCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCgg                                                                                           <  1:531917/140‑1 (MQ=255)
                                              gCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCgg                                                                                           <  1:81339/140‑1 (MQ=255)
                                              gCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCgg                                                                                           >  2:1483233/1‑140 (MQ=255)
                                                   gACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtg                                                                                      >  1:944164/1‑140 (MQ=255)
                                                        gccgccAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtgaaagt                                                                                 <  1:1095297/140‑1 (MQ=255)
                                                        gccgccAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtgaaagt                                                                                 >  1:704273/1‑140 (MQ=255)
                                                                      cATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTgcgc                                                                   <  1:998342/140‑1 (MQ=255)
                                                                       aTTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGcc                                                                  <  1:1586278/140‑1 (MQ=255)
                                                                          cACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATc                                                               >  2:1153033/1‑140 (MQ=255)
                                                                           aCTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGACGTGAAAGTGAATTCTATTGCGCCATCg                                                              >  2:948707/1‑140 (MQ=255)
                                                                           aCTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCg                                                              >  2:1570359/1‑140 (MQ=255)
                                                                            cTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAATGCGTCACTGGTTAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGc                                                             <  2:985964/140‑1 (MQ=255)
                                                                            cTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGc                                                             <  2:319286/140‑1 (MQ=255)
                                                                             tttACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCt                                                            >  2:554701/1‑140 (MQ=255)
                                                                                 ccGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATc                                                        <  1:143601/140‑1 (MQ=255)
                                                                                        tgtgGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCTGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTa                                                 <  2:823070/140‑1 (MQ=255)
                                                                                        tgtgGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCTGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTa                                                 <  2:821695/140‑1 (MQ=255)
                                                                                         gtggtgGAGGGCGGTAGCGAAAAACATTTTGGGTATGCTGCAAGCAAAGCGGCACTGTATAATATGACCCGCTCGTTTGCCCGAAAGCTGGCACCGGATGTTAAAGTTAATTCTATTGCGCCATCGCTGATCCTTTTTat                                                <  2:177811/140‑2 (MQ=255)
                                                                                         gtggtgGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGTAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTaa                                                <  2:176730/140‑1 (MQ=255)
                                                                                                aGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACat                                         >  1:1190032/1‑140 (MQ=255)
                                                                                                           gACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGa                              <  2:1043180/140‑1 (MQ=255)
                                                                                                           gACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGa                              <  1:1129962/140‑1 (MQ=255)
                                                                                                           gACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGa                              <  2:776012/140‑1 (MQ=255)
                                                                                                                      tGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGacaac                   <  1:135719/140‑1 (MQ=255)
                                                                                                                      tGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGacaac                   >  1:1188718/1‑140 (MQ=255)
                                                                                                                      tGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGacaac                   <  1:1553713/140‑1 (MQ=255)
                                                                                                                             gctgcAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCt            <  2:1031515/140‑1 (MQ=255)
                                                                                                                                    gcaaAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGGGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGACTATCGTCAACAGGCGCTGAATaaa     >  2:434776/1‑140 (MQ=255)
                                                                                                                                       aaGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCa  <  1:352460/140‑1 (MQ=255)
                                                                                                                                       aaGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCa  <  1:352560/140‑1 (MQ=255)
                                                                                                                                       aaGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCAATCGCTGATCCTGTTTAACGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCa  <  1:1265585/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                       
ATACCCCATACCTGCTCAACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCA  >  NC_000913/1681307‑1681581

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