Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,400,193 A→T 12.1% intergenic (‑28/+25) nuoG ← / ← nuoF NADH:quinone oxidoreductase subunit G/NADH:quinone oxidoreductase subunit F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,400,1930AT12.1% 273.0 / 20.8 98intergenic (‑28/+25)nuoG/nuoFNADH:quinone oxidoreductase subunit G/NADH:quinone oxidoreductase subunit F
Reads supporting (aligned to +/- strand):  ref base A (49/37);  new base T (5/7);  total (54/44)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TCCCAGCGCCGGATGCCAGCAAAAGTAAGGAATATCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCA  >  NC_000913/2400046‑2400342
                                                                                                                                                   |                                                                                                                                                     
tCCCAGCGCCGGATGCCAGCAAAAGTAAGGAATATCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCagagc                                                                                                                                                     >  2:1406914/1‑149 (MQ=255)
   cAGCGCCGGATGCCAGCAAAAGTAAGGAATATCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTg                                                                                                                                                  <  1:318706/150‑1 (MQ=255)
       gcCGGATGCCAGCAAAAGTAAGGAATATCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCTCTGACTGAgaa                                                                                                                                              <  2:121180/150‑3 (MQ=255)
       gcCGGATGCCAGCAAAAGTAAGGAATATCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCTCTGACTGAgaa                                                                                                                                              <  1:1375277/150‑3 (MQ=255)
        ccGGATGCCAGCAAAAGTAAGGAATATCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCTCTGACTGAgaaa                                                                                                                                             <  1:1614173/150‑4 (MQ=255)
        ccGGATGCCAGCAAAAGTAAGGAATATCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCTCTGACTGAgaaa                                                                                                                                             <  1:351137/150‑4 (MQ=255)
         cGGATGCCATCAAAAGTAAGGAATATCAAGGCCCAGAGACACACAAGCTTCCAGCAGGTTGTAAGATCCGTTCACATCGAAATCTTTGCCGTCTACAAGAATTGTAGCCATTAGCATGCTTACAGATTTCTCAGTCAGAGACTGAGCGtt                                                                                                                                            >  2:138685/1‑150 (MQ=255)
         cGGATGCCAGCAAAAGTAAGGAATATCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCTCTGACTGAgaaaa                                                                                                                                            <  1:603988/150‑5 (MQ=255)
          ggATGCCAGCAAAAGTAAGGAATATCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCTCTGACTGAgaaaac                                                                                                                                           <  1:1611651/150‑6 (MQ=255)
              gCCAGCAAAAGTAAGGAATATCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCg                                                                                                                                       >  2:1579546/1‑150 (MQ=255)
              gCCAGCAAAAGTAAGGAATATCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCg                                                                                                                                       >  2:258069/1‑150 (MQ=255)
                  gCAAAAGTAAGGAATATCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAAt                                                                                                                                   >  1:1620513/1‑150 (MQ=255)
                   cAAAAGTAAGGAATATCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAAtt                                                                                                                                  <  2:1548462/150‑1 (MQ=255)
                      aaGTAAGGAATATCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCgg                                                                                                                               <  1:258069/150‑1 (MQ=255)
                        gTAAGGAATATCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGtt                                                                                                                             <  1:1424791/150‑1 (MQ=255)
                             gAATATCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAg                                                                                                                        <  1:1249208/150‑1 (MQ=255)
                             gAATATCAAGGCCAAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTTACCTCGTATTCTTTGACGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAg                                                                                                                        >  1:760318/1‑150 (MQ=255)
                                tatCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAgcgc                                                                                                                     <  1:196202/150‑1 (MQ=255)
                                 atCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCACTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGt                                                                                                                                             <  1:985837/125‑1 (MQ=255)
                                 atCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCACTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGt                                                                                                                                             >  2:985837/1‑125 (MQ=255)
                                   cAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGctct                                                                                                                  >  2:57545/1‑150 (MQ=255)
                                        cccAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCTCTGACTGAgaaa                                                                                                                                             <  1:165804/118‑4 (MQ=255)
                                        cccAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCTCTGACTGAgaaa                                                                                                                                             >  2:165804/1‑115 (MQ=255)
                                         ccAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTcagc                                                                                                            <  1:1025874/150‑1 (MQ=255)
                                         ccAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTcagc                                                                                                            >  2:1285806/1‑150 (MQ=255)
                                                cagacaAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCg                                                                                                     >  2:1564252/1‑150 (MQ=255)
                                                cagacaAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCg                                                                                                     >  2:148791/1‑150 (MQ=255)
                                                cagacaAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCg                                                                                                     >  1:996490/1‑150 (MQ=255)
                                                  gacaAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGc                                                                                                   <  1:1350062/150‑1 (MQ=255)
                                                      aGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGaa                                                                                               >  2:1014946/1‑150 (MQ=255)
                                                       gCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAAt                                                                                              <  1:311780/150‑1 (MQ=255)
                                                         ttCCATCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATcc                                                                                            >  1:653253/1‑150 (MQ=255)
                                                         ttCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATcc                                                                                            <  2:36155/150‑1 (MQ=255)
                                                           ccAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAgc                                                                                                                       >  2:636792/1‑121 (MQ=255)
                                                           ccAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAgc                                                                                                                       <  1:636792/121‑1 (MQ=255)
                                                           ccAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCa                                                                                          >  2:505759/1‑150 (MQ=255)
                                                           ccAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAACAGGTTCGGCTGAATCCCa                                                                                          >  2:26969/1‑150 (MQ=255)
                                                            catcagGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCTTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCAt                                                                                         <  2:924803/147‑1 (MQ=255)
                                                             agcagGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCAtt                                                                                        >  1:530550/1‑150 (MQ=255)
                                                              gcagGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTa                                                                                       <  1:856563/150‑1 (MQ=255)
                                                              gcagGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTa                                                                                       >  1:648980/1‑150 (MQ=255)
                                                              gcagGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTa                                                                                       <  1:198151/150‑1 (MQ=255)
                                                                    tGTCCGCTCCGTTGACCTCCTATTATTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGCATGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGGTCGGCTGAATCCCATGAATCaaa                                                                                 >  2:1392507/1‑150 (MQ=255)
                                                                     gTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGATTTGTTGCCATTAGCATGCTTCCAGGTTTCTCAGTCAGAGCCTGAGCGTTAATCGATATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAAt                                                                                >  2:1138378/1‑150 (MQ=255)
                                                                     gTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAAt                                                                                >  2:690918/1‑150 (MQ=255)
                                                                             cGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAAt                                                                                <  2:1576867/142‑1 (MQ=255)
                                                                             cGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAAt                                                                                >  1:1576867/1‑142 (MQ=255)
                                                                               ttGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCt                                                                      <  1:1058341/150‑1 (MQ=255)
                                                                               ttGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCt                                                                      <  1:505759/150‑1 (MQ=255)
                                                                                   ccTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGATCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGTTATTCCTGAAc                                                                  >  2:277055/1‑150 (MQ=255)
                                                                                    cTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACg                                                                 >  2:1375425/1‑150 (MQ=255)
                                                                                    cTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACg                                                                 <  2:1620513/150‑1 (MQ=255)
                                                                                      cGTATTCTTTGTCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGc                                                               <  1:336182/150‑1 (MQ=255)
                                                                                        tATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTg                                                             >  1:799914/1‑150 (MQ=255)
                                                                                        tATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTg                                                             >  1:1401913/1‑150 (MQ=255)
                                                                                           tCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGttt                                                          <  1:878676/150‑1 (MQ=255)
                                                                                           tCTTTGACGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGttt                                                          <  1:762382/150‑1 (MQ=255)
                                                                                            cTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTg                                                         <  1:1296263/150‑1 (MQ=255)
                                                                                                 ccGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTcc                                                    >  1:1158651/1‑150 (MQ=255)
                                                                                                   gTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCg                                                  >  1:1303082/1‑150 (MQ=255)
                                                                                                          tGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCaaa                                           >  1:797964/1‑150 (MQ=255)
                                                                                                           gAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAAt                                          >  1:348794/1‑150 (MQ=255)
                                                                                                               tGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAAttctt                                      >  2:999541/1‑150 (MQ=255)
                                                                                                                gTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAAttcttc                                     >  1:1311542/1‑150 (MQ=255)
                                                                                                                 tAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTcg                                    >  1:1037857/1‑150 (MQ=255)
                                                                                                                  tgCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTcgc                                   <  1:1359782/149‑1 (MQ=255)
                                                                                                                   gCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTcgcg                                  >  1:1298413/1‑150 (MQ=255)
                                                                                                                   gCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTcgcg                                  <  2:1117689/150‑1 (MQ=255)
                                                                                                                    ccATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCgg                                 >  1:1503583/1‑150 (MQ=255)
                                                                                                                     cATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGa                                >  1:540962/1‑150 (MQ=255)
                                                                                                                           cATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATAtt                          <  2:629663/150‑1 (MQ=255)
                                                                                                                           cATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATAtt                          >  2:894849/1‑150 (MQ=255)
                                                                                                                            aTGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATAttt                         <  1:690918/150‑1 (MQ=255)
                                                                                                                            aTGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATAttt                         <  2:797964/150‑1 (MQ=255)
                                                                                                                              gCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGa                       >  2:1470918/1‑150 (MQ=255)
                                                                                                                              gCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGa                       >  1:1457022/1‑150 (MQ=255)
                                                                                                                               cTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGAt                      <  1:1375425/150‑1 (MQ=255)
                                                                                                                               cTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGAt                      >  1:813224/1‑150 (MQ=255)
                                                                                                                               cTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGAt                      <  1:44352/150‑1 (MQ=255)
                                                                                                                               cTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGAt                      >  2:1432147/1‑150 (MQ=255)
                                                                                                                                ttCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATg                     >  2:957309/1‑150 (MQ=255)
                                                                                                                                ttCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATg                     >  1:673764/1‑150 (MQ=255)
                                                                                                                                  ccAGTTTTCTCAGTCAGAGACTGAGCTTTAATCGAAATTCGGTTTCCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGgc                   <  1:138685/150‑1 (MQ=255)
                                                                                                                                  ccAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGgc                   >  1:630305/1‑150 (MQ=255)
                                                                                                                                    ttaacgctcAGTCTCTGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGgcgc                 >  1:1523130/7‑150 (MQ=255)
                                                                                                                                      aacgctcAGTCTCTGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGctc               >  2:989460/5‑150 (MQ=255)
                                                                                                                                       acgctcAGTCTCTGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGctct              >  2:385080/4‑150 (MQ=255)
                                                                                                                                         tctcAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGt            >  1:949176/1‑150 (MQ=255)
                                                                                                                                         gctcAGTCTCTGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGt            >  2:735426/2‑150 (MQ=255)
                                                                                                                                          ctcAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTa           <  2:540962/150‑1 (MQ=255)
                                                                                                                                          ctcAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTa           <  2:813224/150‑1 (MQ=255)
                                                                                                                                          ctcAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTa           >  1:498194/1‑150 (MQ=255)
                                                                                                                                            cagtcagAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAAc         >  1:381345/1‑150 (MQ=255)
                                                                                                                                            cagtcagAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAAc         >  1:886899/1‑150 (MQ=255)
                                                                                                                                            cagtcagAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAAc         >  1:581787/1‑150 (MQ=255)
                                                                                                                                             agtcagAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACg        <  2:673764/150‑1 (MQ=255)
                                                                                                                                              gtcagAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACgg       <  1:1342763/150‑1 (MQ=255)
                                                                                                                                               tcagAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGc      <  2:1303082/150‑1 (MQ=255)
                                                                                                                                               tcagAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAACAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGc      <  1:26969/150‑1 (MQ=255)
                                                                                                                                                   agaCTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCa  <  2:1457022/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCCCAGCGCCGGATGCCAGCAAAAGTAAGGAATATCAAGGCCCAGAGACAGACAAGCTTCCAGCAGGTTGTCCGCTCCGTTGACCTCGTATTCTTTGCCGTCTACATGAATTGTAGCCATTAGCATGCTTCCAGTTTTCTCAGTCAGAGACTGAGCGTTAATCGAAATTCGGTTACCAGCGCTCTTTCAGCAGGTTCGGCTGAATCCCATTAATCAAATGGGTATTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCA  >  NC_000913/2400046‑2400342

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: ‑