Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,903,376 T→C intergenic (+2/+1371) ygcE → / ← queE putative sugar kinase YgcE/putative 7‑carboxy‑7‑deazaguanine synthase QueE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,903,3760TC98.5% 230.5 / NA 68intergenic (+2/+1371)ygcE/queEputative sugar kinase YgcE/putative 7‑carboxy‑7‑deazaguanine synthase QueE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (35/32);  minor base G (1/0);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCGATGAATAAAGGCATTTTCAAAGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTATCCCCGCTGATGCAGGG  >  NC_000913/2903247‑2903504
                                                                                                                                 |                                                                                                                                
acgCGATGAATAAAGGCATTTTCAAAGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGcc                                                                                                                        <  1:1308941/140‑1 (MQ=255)
     aTGAATAAAGGCATTTTCAAAGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCGAGAGCTGGTCGAATGCGTAAGGGAGGAAGGCcgggcg                                                                                                                   >  1:510523/1‑140 (MQ=255)
      tGAATAAAGGCATTTTCAAAGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCcgggcgg                                                                                                                  >  1:114141/1‑140 (MQ=255)
              ggCATTTTCAAAGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGACGAATGCGTAAGCGGGGAAGGTCGGGCGGGAAACTGc                                                                                                          >  1:450845/1‑140 (MQ=255)
               gCATTTTCAAAGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGcc                                                                                                         <  2:330889/140‑1 (MQ=255)
               gCATTTTCAAAGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGcc                                                                                                         <  1:1335010/140‑1 (MQ=255)
                  ttttCAAAGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCgg                                                                                                      >  1:773051/1‑140 (MQ=255)
                        aaGACCTAACCAAACATACAGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGaa                                                                                                <  2:1334963/140‑1 (MQ=255)
                             cTAACCAAACATACTGACGTGATCCTGAAAACATCCTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATAc                                                                                           <  2:609917/140‑1 (MQ=255)
                               aaCCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCt                                                                                         >  2:783989/1‑140 (MQ=255)
                               aaCCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCt                                                                                         >  1:1294012/1‑140 (MQ=255)
                               aaCCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCt                                                                                         >  1:30228/1‑140 (MQ=255)
                               aaCCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCt                                                                                         >  1:382728/1‑140 (MQ=255)
                                 ccAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGa                                                                                       >  2:956096/1‑140 (MQ=255)
                                       aTACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGtt                                                                                 <  2:1294012/140‑1 (MQ=255)
                                            gATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATccc                                                                            >  2:1226708/1‑140 (MQ=255)
                                                 gATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTg                                                                       >  2:1164673/1‑140 (MQ=255)
                                                 gATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTg                                                                       >  1:308472/1‑140 (MQ=255)
                                                  aTCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCCGAAGATGCCTGAAAGGTTATCCCCGCTGa                                                                      >  1:434858/1‑140 (MQ=255)
                                                   tCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGAc                                                                     >  2:832630/1‑140 (MQ=255)
                                                   tCCTGAAAAAATCGTATGAAGTGATGCATGGGCAATTGGGTAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGCTATCCCCGCTGAc                                                                     >  1:257961/1‑140 (MQ=255)
                                                     cTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGAcgc                                                                   >  2:346397/1‑140 (MQ=255)
                                                      tGAAAAAATCGTATGAAGTGATGCATGGGGATTTGGGGATTGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGAcgcg                                                                  <  2:593089/140‑1 (MQ=255)
                                                      tGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGAcgcg                                                                  <  1:1437133/140‑1 (MQ=255)
                                                      tGAAAAAATCGTATGAAGTGATCCATGGGGACTTGAGGAATGTGGATTTGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGAcgcg                                                                  <  2:450845/140‑1 (MQ=255)
                                                       gAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACGTACGTGGATGGTTATCCCCGCTGACGCgg                                                                 >  1:568012/1‑140 (MQ=255)
                                                       gAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGAcg                                                                    >  2:898601/1‑137 (MQ=255)
                                                       gAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGAcg                                                                    <  1:898601/137‑1 (MQ=255)
                                                       gAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCgg                                                                 >  2:1405604/1‑140 (MQ=255)
                                                       gAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCgg                                                                 >  1:892771/1‑140 (MQ=255)
                                                        aaaaaaTCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCggg                                                                <  1:832630/140‑1 (MQ=255)
                                                        aaaaaaTCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCggg                                                                <  1:346397/140‑1 (MQ=255)
                                                        aaaaaaTCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCggg                                                                <  2:1194947/140‑1 (MQ=255)
                                                            aaTCGTATGAAGTGATGTATGGCGAATTGGGCAATCTGGTTTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAAc                                                            <  2:400372/140‑1 (MQ=255)
                                                            aaTCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAAc                                                            <  2:1257060/140‑1 (MQ=255)
                                                                       gTGATGCATGGGGAATTTGGGACTGTGGATTCGATCCAGGGCTGGTCGCATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGAc                                                 <  2:510523/140‑1 (MQ=255)
                                                                       gTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGAc                                                 <  2:99681/140‑1 (MQ=255)
                                                                               tggggaattggggaatGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTaaa                                         >  2:1266435/1‑140 (MQ=255)
                                                                               tggggaattggggaatGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTaaa                                         >  2:331342/1‑140 (MQ=255)
                                                                                ggggaattggggaatGTGGATTCGATCCAGAGGTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTaaaa                                        >  2:1057075/1‑140 (MQ=255)
                                                                                 gggaattggggaatGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTaaaaa                                       <  1:804687/140‑1 (MQ=255)
                                                                                  ggaattggggaatGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTaaaaaa                                      >  2:361860/1‑140 (MQ=255)
                                                                                    aattggggaatGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATg                                    >  2:779301/1‑140 (MQ=255)
                                                                                        gaggaatGTGGATTCGATCGAGAGCTGGTCGAATGTGTAAGCGAGGAAGGTCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGTCGCGGGGAACATAAGTTCGACAGCTTAAAAAATGATaa                                <  2:257961/138‑1 (MQ=255)
                                                                                         gggaatGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAAc                               <  1:361860/140‑1 (MQ=255)
                                                                                            aatGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAg                            <  2:773051/140‑1 (MQ=255)
                                                                                             atGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGt                           <  2:785384/140‑1 (MQ=255)
                                                                                              tgtgGATTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTc                          <  1:1262521/140‑1 (MQ=255)
                                                                                                   aTTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGCCAAGTTCAAAAATGATAACAAGTCCGGtt                     >  1:660023/1‑140 (MQ=255)
                                                                                                   aTTCGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGtt                     >  2:1313353/1‑140 (MQ=255)
                                                                                                      cGATCCAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTAt                  >  2:366683/1‑140 (MQ=255)
                                                                                                           cAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTATCCCCg             <  1:366683/140‑1 (MQ=255)
                                                                                                           cAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTATCCCCg             <  1:331342/140‑1 (MQ=255)
                                                                                                           cAGAGCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTATCCCCg             <  1:1266435/140‑1 (MQ=255)
                                                                                                            agagCTGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGAAAACTTAAGAAAGGAGAACAAGACCGGTTTATCCCCGc            >  1:691794/1‑140 (MQ=255)
                                                                                                            agagCTGGTCGAATGCGAAAGCTAGGAAGGCCGGGCGGGTAACTGCCCAGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTATCCCCGc            <  2:638252/140‑1 (MQ=255)
                                                                                                                 tGGTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTATCCCCGCTGATg       <  2:64566/140‑1 (MQ=255)
                                                                                                                   gTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGGGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGACCGGTTTATCCCCGCTGATGCa     >  2:170800/1‑140 (MQ=255)
                                                                                                                   gTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACTAGTCCGGTTTATCCCCGCTGATGCa     >  1:638252/1‑140 (MQ=255)
                                                                                                                   gTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTATCCCCGCTg          >  2:1050287/1‑135 (MQ=255)
                                                                                                                   gTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTATCCCCGCTg          <  1:1050287/135‑1 (MQ=255)
                                                                                                                   gTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTATCCCCGCTGATGCa     >  2:126933/1‑140 (MQ=255)
                                                                                                                   gTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTATCCCCGCTGATGCa     >  2:34713/1‑140 (MQ=255)
                                                                                                                   gTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTATCCCCGCTGATGCa     >  1:1159380/1‑140 (MQ=255)
                                                                                                                   gTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTATCCCCGCTGATGCa     >  2:960334/1‑140 (MQ=255)
                                                                                                                   gTCGAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTCAAAAAATGATAACAAGCCCGGATTATCCCCGCTGATGCa     >  2:605707/1‑140 (MQ=255)
                                                                                                                     ggAATGCGTAAGCGAGTGAGGCTGGGCGGGAGACTACCCTGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTATCCCCGCTGATGCAgg   <  2:690702/139‑1 (MQ=255)
                                                                                                                     cGAATGCGTAAGCGATGAAGGCAGGCCGGAAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTATCCCCGCTGATGCAgg   <  1:605707/140‑1 (MQ=255)
                                                                                                                     cGAATGCGTAAGCGAGGCAGGCCGGGCGGGATAGTACCCAGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTATCCCCGCTGATGCAgg   <  2:691794/140‑1 (MQ=255)
                                                                                                                      gAATGCGTAAGCGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTATCCCCGCTGATGCAggg  <  1:1313353/140‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                
ACGCGATGAATAAAGGCATTTTCAAAGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTATCCCCGCTGATGCAGGG  >  NC_000913/2903247‑2903504

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