Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 341,717 A→T 6.6% N198I (AAT→ATT)  yahJ → putative deaminase with metallo‑dependent hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913341,7170AT 6.6% 314.2 / 7.6 100N198I (AAT→ATT) yahJputative deaminase with metallo‑dependent hydrolase domain
Reads supporting (aligned to +/- strand):  ref base A (51/42);  new base T (5/2);  total (56/44)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTCGAGCAGAAAATGCTGCCGGAACTGCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGG  >  NC_000913/341570‑341860
                                                                                                                                                   |                                                                                                                                               
aaCTCGAGCAGAAAATGCTGCCGGAACTGCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAAtt                                                                                                                                               >  2:796215/1‑150 (MQ=255)
aaCTCGAGCAGAAAATGCTGCCGGAACTGCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAAtt                                                                                                                                               >  2:790371/1‑150 (MQ=255)
   acGAGCAGAAAATGCTGCCGGAACTGCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGc                                                                                                                                            <  2:847975/149‑1 (MQ=255)
    cGAGCAGAAAATGCTGCCGGAACTGCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCa                                                                                                                                           >  1:1258440/1‑150 (MQ=255)
     gAGCAGAAAATGCTGCCGGAACTGCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAg                                                                                                                                          >  2:950386/1‑150 (MQ=255)
      aGCAGAAAATGCTGCCGGAACTGCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAgg                                                                                                                                         <  2:450471/150‑1 (MQ=255)
       gCAGAAAATGCTGCCGGAACTGCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGc                                                                                                                                        >  2:770786/1‑150 (MQ=255)
        cAGAAAATGCTGCCGGAACTGCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCCCCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCGGAAAAATCTGCAAAATTTGCAGGCg                                                                                                                                       >  2:180657/1‑150 (MQ=255)
        cAGAAAATGCTGCCGGAACTGCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAATTTTGCAgatt                                                                                                                                       <  1:1364951/150‑4 (MQ=255)
           aaaaTGCTGCCGGAACTGCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTg                                                                                                                                    <  2:997429/150‑1 (MQ=255)
                 ctgcCGGAACTGCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGgcg                                                                                                                              >  1:1078795/1‑150 (MQ=255)
                 ctgcCGGAACTGCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGgcg                                                                                                                              <  2:153787/150‑1 (MQ=255)
                   gcCGGAACTGCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGcgcg                                                                                                                            <  1:770786/150‑1 (MQ=255)
                      ggAACTGCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACg                                                                                                                         <  2:894168/150‑1 (MQ=255)
                      ggAACTGCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACg                                                                                                                         >  1:338107/1‑150 (MQ=255)
                      ggAACTGCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACg                                                                                                                         <  1:541461/150‑1 (MQ=255)
                          cTGCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAg                                                                                                                     >  1:705067/1‑150 (MQ=255)
                            gCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGc                                                                                                                   <  2:1078795/150‑1 (MQ=255)
                              aaCCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCgg                                                                                                                 >  1:845332/1‑150 (MQ=255)
                                ccGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGc                                                                                                               <  1:705073/150‑1 (MQ=255)
                                    acacTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTg                                                                                                           >  2:71192/1‑150 (MQ=255)
                                    acacTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTg                                                                                                           <  2:638721/150‑1 (MQ=255)
                                    acacTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTg                                                                                                           <  1:1038977/150‑1 (MQ=255)
                                    acacTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTg                                                                                                           >  2:1038977/1‑150 (MQ=255)
                                    acacTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTg                                                                                                           <  1:1025849/150‑1 (MQ=255)
                                    acacTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGCTGCTGGCGCGACGTCAGGCGGGc                                                                                                               <  1:815795/146‑1 (MQ=255)
                                    acacTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGCTGCTGGCGCGACGTCAGGCGGGc                                                                                                               >  2:815795/1‑146 (MQ=255)
                                     cacTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGa                                                                                                          >  1:368040/1‑150 (MQ=255)
                                     cacTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGa                                                                                                          >  1:752143/1‑150 (MQ=255)
                                     cacTCAGCAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGa                                                                                                          >  1:1152082/1‑150 (MQ=255)
                                      acTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAg                                                                                                         <  2:1365823/150‑1 (MQ=255)
                                         cAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAgtgt                                                                                                      >  1:116717/1‑150 (MQ=255)
                                           ggAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGa                                                                                                    >  2:1121493/1‑150 (MQ=255)
                                                 ggCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAAtcgt                                                                                              >  1:397400/1‑150 (MQ=255)
                                                 ggCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAAtcgt                                                                                              >  2:1043954/1‑150 (MQ=255)
                                                   cAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAAtcgtcg                                                                                            >  1:1471391/1‑150 (MQ=255)
                                                     gAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGcc                                                                                          >  2:602419/1‑150 (MQ=255)
                                                      aaaaaCTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCt                                                                                         >  2:1435678/1‑150 (MQ=255)
                                                       aaaaCTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCtt                                                                                        >  2:226877/1‑150 (MQ=255)
                                                           ttgattgatTTATTTCAGTCGAATGGGAGCTGCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCGGGCCTGTAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTTAAATCGTCGCCTTCCCg                                                                                    <  2:1601278/149‑1 (MQ=255)
                                                            tgattgatTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCgc                                                                                   <  2:256355/150‑1 (MQ=255)
                                                                 gattTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCAc                                                                              >  1:1146288/1‑150 (MQ=255)
                                                                  attTATTGCAGTCGAAAGGCCCCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACg                                                                             >  1:437500/1‑150 (MQ=255)
                                                                      aTTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGttt                                                                         <  2:116717/150‑1 (MQ=255)
                                                                        tGCAGTCGAAACGCACCACCATTGCCCGCAGCCACTGCAAAATCGAACCTGTTTCCGGCCTGAAAAATCTGCAAAATTTGCACGCCCTGCTGGCGCGACCTCAGGCGGGCTTAGAGTGTAAAATCGTCGCCTACCCGCAGCACAGTTtga                                                                       >  2:820442/1‑149 (MQ=255)
                                                                         gCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTtgct                                                                      >  1:1356478/1‑150 (MQ=255)
                                                                          cAGTCGAAAGGCCCCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTtgctg                                                                     <  2:834436/150‑1 (MQ=255)
                                                                          cAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTtgctg                                                                     >  1:121778/1‑150 (MQ=255)
                                                                          cAGGCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGCACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTGTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTtgctg                                                                     <  2:828787/150‑1 (MQ=255)
                                                                           aGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGCCGTCAGGCGGGCTTTGAGTGGGAAATCGTCGCCTTCCCGCAGCACGGTTtgctgc                                                                    <  2:1193917/150‑1 (MQ=255)
                                                                            gTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTgctgct                                                                   >  1:268045/1‑150 (MQ=255)
                                                                            gTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTgctgct                                                                   >  2:530701/1‑150 (MQ=255)
                                                                               gAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTc                                                                <  1:226877/150‑1 (MQ=255)
                                                                                 aaGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGa                                                              >  2:1264862/1‑150 (MQ=255)
                                                                                     caccacCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATc                                                          >  2:639692/1‑150 (MQ=255)
                                                                                       ccaccaTTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGtgtg                                                                                                     >  1:703692/1‑105 (MQ=255)
                                                                                       ccaccaTTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGtgtg                                                                                                     <  2:703692/105‑1 (MQ=255)
                                                                                       ccaccaTTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTg                                                        <  2:752143/150‑1 (MQ=255)
                                                                                         accaTTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGaa                                                      >  1:462449/1‑150 (MQ=255)
                                                                                         accaTTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGaa                                                      >  2:871724/1‑150 (MQ=255)
                                                                                          ccaTTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAAc                                                     <  1:790371/150‑1 (MQ=255)
                                                                                          ccaTTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAAc                                                     <  1:796215/150‑1 (MQ=255)
                                                                                           caTTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAAcc                                                    <  2:619808/150‑1 (MQ=255)
                                                                                              tGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTTAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCttt                                                 <  2:397400/150‑1 (MQ=255)
                                                                                              tGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCttt                                                 <  1:56639/150‑1 (MQ=255)
                                                                                              tGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCttt                                                 >  1:190716/1‑150 (MQ=255)
                                                                                              tGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCttt                                                 <  2:268045/150‑1 (MQ=255)
                                                                                               ggccGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGGGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTa                                                <  1:1054211/148‑1 (MQ=255)
                                                                                               gCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTa                                                >  2:226089/1‑150 (MQ=255)
                                                                                               gCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTa                                                >  1:596368/1‑150 (MQ=255)
                                                                                                cccGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTaa                                               >  2:1295421/1‑150 (MQ=255)
                                                                                                  cGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATg                                             <  1:132027/150‑1 (MQ=255)
                                                                                                    cAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCg                                           >  1:1293235/1‑150 (MQ=255)
                                                                                                    cAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCg                                           <  1:639692/150‑1 (MQ=255)
                                                                                                      gCCACTGCAATATCGAACGGGTGTCCGGCCTTAAAAATCTGCAAATTTTGCAGGCTGTGCTGGCGCTACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCACCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTg                                         <  2:1258440/150‑1 (MQ=255)
                                                                                                      gCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTg                                         <  2:368040/150‑1 (MQ=255)
                                                                                                      gCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCGTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTg                                         <  1:649030/150‑1 (MQ=255)
                                                                                                       ccACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAAGCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGa                                        <  1:558752/150‑1 (MQ=255)
                                                                                                        cACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGaa                                       <  2:300095/150‑1 (MQ=255)
                                                                                                            gCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGAAGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGa                                   >  1:1252118/1‑150 (MQ=255)
                                                                                                              aaTATGGAACGTGTTTCCGGCCTGAAAAATCTGCAAAATTTGGAGGCGGTGCGGGCGCGACGGCAGGTGGGCTTTGAGTGTGAAATCGTCGGCTTCCCGCAGCACGGTTTGCTGCTGACGAAATCTTAACCTTTAATGGGTGAAGCGATg                                 <  2:1036171/150‑1 (MQ=255)
                                                                                                              aaTATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATg                                 >  2:1590359/1‑150 (MQ=255)
                                                                                                              aaTATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATg                                 <  2:1152082/150‑1 (MQ=255)
                                                                                                                 atCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAg                              <  2:845332/150‑1 (MQ=255)
                                                                                                                   cGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGc                            >  1:454944/1‑150 (MQ=255)
                                                                                                                    gAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCg                           >  1:898871/1‑150 (MQ=255)
                                                                                                                     aaCCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCgg                          >  2:45292/1‑150 (MQ=255)
                                                                                                                      aCCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCggg                         <  2:144674/150‑1 (MQ=255)
                                                                                                                         ggTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGgc                      <  1:45292/150‑1 (MQ=255)
                                                                                                                         ggTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGgc                      >  2:239218/1‑150 (MQ=255)
                                                                                                                         ggTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGgc                      <  1:256520/150‑1 (MQ=255)
                                                                                                                                 gCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACg              >  1:1162429/1‑150 (MQ=255)
                                                                                                                                    tGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTcg           >  2:1102841/1‑150 (MQ=255)
                                                                                                                                        ccgccTGCAAATTTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTcggcgg       >  1:867012/5‑150 (MQ=255)
                                                                                                                                        ccgccTGCAAATTTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTcggcgg       >  2:1480311/5‑150 (MQ=255)
                                                                                                                                        ccgccTGCAAATTTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTcggcgg       >  2:1259566/5‑150 (MQ=255)
                                                                                                                                         cgccTGCAAATTTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCggcggc      >  1:1374597/4‑150 (MQ=255)
                                                                                                                                         cgccTGCAAATTTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCggcggc      >  2:970646/4‑150 (MQ=255)
                                                                                                                                             tGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTgg  >  2:752160/1‑150 (MQ=255)
                                                                                                                                             tGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTgg  >  1:401046/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
AACTCGAGCAGAAAATGCTGCCGGAACTGCAACCGTACACTCAGGAGCGGGCAGAAAAACTGATTGATTTATTGCAGTCGAAAGGCACCACCATTGCCCGCAGCCACTGCAATATCGAACCGGTTTCCGGCCTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGG  >  NC_000913/341570‑341860

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