Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,399,880 T→C 8.2% M96V (ATG→GTG) ‡ nuoG ← NADH:quinone oxidoreductase subunit G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,399,8800TC 8.2% 245.6 / 10.2 85M96V (ATG→GTG) ‡nuoGNADH:quinone oxidoreductase subunit G
Reads supporting (aligned to +/- strand):  ref base T (43/35);  new base C (2/5);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

GAGATGAATGGCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTATTGCTTCACCGCACACTGG  >  NC_000913/2399732‑2400028
                                                                                                                                                    |                                                                                                                                                    
gagaTGAATGGCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGtcatc                                                                                                                                                     >  1:481157/1‑150 (MQ=255)
gagaTGAATGGCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGtcatc                                                                                                                                                     >  1:702399/1‑150 (MQ=255)
gagaTGAATGGCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAAAGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGtcatc                                                                                                                                                     <  2:467833/150‑1 (MQ=255)
    tGAATGGCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAAcc                                                                                                                                                 <  1:113531/150‑1 (MQ=255)
    tGAATGGCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAAcc                                                                                                                                                 >  1:1312810/1‑150 (MQ=255)
     gAATGGCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCa                                                                                                                                                <  1:72354/150‑1 (MQ=255)
     gAATGGCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCa                                                                                                                                                >  2:1627217/1‑150 (MQ=255)
      aaTGGCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCaagccc                                                                                                                                               >  2:214261/1‑146 (MQ=255)
        tGGCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTc                                                                                                                                             >  1:1195011/1‑150 (MQ=255)
        tGGCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGATGACCAACCAccc                                                                                                                                             <  2:458352/150‑3 (MQ=255)
         ggCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGATGACCAACCACCCg                                                                                                                                            <  1:1455829/150‑1 (MQ=255)
          gCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGAGTGGTTGATGACCAACCACCCgc                                                                                                                                           <  1:1521467/150‑2 (MQ=255)
          gCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGAGTGGTTGATGACCAACCACCCgc                                                                                                                                           <  1:1515291/150‑2 (MQ=255)
             ccAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCa                                                                                                                                        <  2:638630/150‑1 (MQ=255)
              cAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGATGACCAACCACCCgcacga                                                                                                                                       <  1:1260757/150‑6 (MQ=255)
                    cTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTc                                                                                                                                 >  1:483423/1‑150 (MQ=255)
                         tACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGa                                                                                                                            <  2:827955/150‑1 (MQ=255)
                         tACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGa                                                                                                                            <  2:905430/150‑1 (MQ=255)
                          aCGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGaa                                                                                                                           >  2:1618011/1‑150 (MQ=255)
                           cGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAAc                                                                                                                          >  1:1447780/1‑150 (MQ=255)
                               gggTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGtt                                                                                                                      >  2:19490/1‑150 (MQ=255)
                                 gTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTc                                                                                                                    <  1:1083216/150‑1 (MQ=255)
                                 gTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTc                                                                                                                    <  2:1447780/150‑1 (MQ=255)
                                  tACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCg                                                                                                                   <  1:1090408/150‑1 (MQ=255)
                                  tACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCg                                                                                                                   >  1:1517682/1‑150 (MQ=255)
                                  tACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCg                                                                                                                   <  1:695942/150‑1 (MQ=255)
                                  tACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCg                                                                                                                   <  1:19490/150‑1 (MQ=255)
                                    cGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGct                                                                                                                 >  2:965939/1‑150 (MQ=255)
                                        tGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGcttctt                                                                                                             >  1:720611/1‑150 (MQ=255)
                                        tGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGCTTGGTCATCAACAACTCGACCACGCTTTCACGGAACTGTTTCGcttctt                                                                                                             >  2:283527/1‑150 (MQ=255)
                                          gTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTtcg                                                                                                           >  2:56864/1‑150 (MQ=255)
                                          gTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTtcg                                                                                                           <  2:483423/150‑1 (MQ=255)
                                          gTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTtcg                                                                                                           >  1:1151612/1‑150 (MQ=255)
                                              aaCGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTcgtcgt                                                                                                       <  1:949250/150‑1 (MQ=255)
                                                    aGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATgg                                                                                                 <  2:7098/150‑1 (MQ=255)
                                                      cGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGaa                                                                                               <  1:1239495/150‑1 (MQ=255)
                                                         cGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATa                                                                                            >  2:1624301/1‑150 (MQ=255)
                                                                tgtgTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGc                                                                                     >  1:1639841/1‑150 (MQ=255)
                                                                 gtgtCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGcc                                                                                    >  2:1636252/1‑150 (MQ=255)
                                                                   gtCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCAt                                                                                  >  1:524518/1‑150 (MQ=255)
                                                                     ccGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCg                                                                                <  2:524518/150‑1 (MQ=255)
                                                                     ccGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCg                                                                                >  2:1590952/1‑150 (MQ=255)
                                                                            tcaCAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCg                                                                         >  1:707281/1‑150 (MQ=255)
                                                                            tcaCAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCg                                                                         >  1:582876/1‑150 (MQ=255)
                                                                            tcaCAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCg                                                                         >  1:444985/1‑150 (MQ=255)
                                                                            tcaCAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTATCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCg                                                                         >  2:374879/1‑150 (MQ=255)
                                                                              acaGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGgt                                                                       >  2:1146232/1‑150 (MQ=255)
                                                                                aGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGgtgt                                                                     <  1:965939/150‑1 (MQ=255)
                                                                                aGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGgtgt                                                                     <  2:1295968/150‑1 (MQ=255)
                                                                                 gTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTc                                                                    >  1:618096/1‑150 (MQ=255)
                                                                                  tCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCa                                                                   <  1:1506908/150‑1 (MQ=255)
                                                                                  tCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCa                                                                   <  1:1590611/150‑1 (MQ=255)
                                                                                     tatCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACgg                                                                                                                             <  1:1178137/89‑1 (MQ=255)
                                                                                     tatCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACgg                                                                                                                             >  2:1178137/1‑89 (MQ=255)
                                                                                         gtGAAGAGGGCGGTAACTGCCATCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGa                                                            >  2:773202/2‑150 (MQ=255)
                                                                                          tGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGAc                                                           >  1:1039649/1‑150 (MQ=255)
                                                                                          tGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGAc                                                           <  2:1312810/150‑1 (MQ=255)
                                                                                          tGAAGAGGGCGGTAACTGCCATCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGAc                                                           >  2:521325/1‑150 (MQ=255)
                                                                                           gAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACa                                                          >  2:1233250/1‑150 (MQ=255)
                                                                                           gAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACa                                                          >  1:745456/1‑150 (MQ=255)
                                                                                                tGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCAcc                                                     >  2:1295346/1‑150 (MQ=255)
                                                                                                 ggCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCa                                                    >  1:262323/1‑150 (MQ=255)
                                                                                                 ggCAGTTACCGACCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGGCATCAACCACTCGACCACGCTTTCACGGAAGTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCa                                                    <  1:1058859/150‑1 (MQ=255)
                                                                                                    aGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGc                                                 >  2:1409264/1‑150 (MQ=255)
                                                                                                     gTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCg                                                <  1:56864/150‑1 (MQ=255)
                                                                                                      ttACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGa                                               <  1:1618011/150‑1 (MQ=255)
                                                                                                        aCCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGAcc                                             <  1:1233250/150‑1 (MQ=255)
                                                                                                           gCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGCGTCATACAGGACATCACCAGGCGACCAcg                                          >  2:411347/1‑150 (MQ=255)
                                                                                                            cccTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCAcgc                                         >  2:1653664/1‑150 (MQ=255)
                                                                                                            cccTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTACTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCTcgc                                         >  2:83534/1‑150 (MQ=255)
                                                                                                             ccTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCAcgcg                                        <  2:702399/150‑1 (MQ=255)
                                                                                                             ccTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCAcgcg                                        >  1:49925/1‑150 (MQ=255)
                                                                                                                       tACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCcgc                              >  1:1660480/1‑150 (MQ=255)
                                                                                                                        aCCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCcgcg                             >  2:200769/1‑150 (MQ=255)
                                                                                                                         ccGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGt                            >  2:970274/1‑150 (MQ=255)
                                                                                                                              cAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTgg                       <  1:1627217/150‑1 (MQ=255)
                                                                                                                                gTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGgtg                                                                      >  2:43263/1‑101 (MQ=255)
                                                                                                                                gTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGgtg                                                                      <  1:43263/101‑1 (MQ=255)
                                                                                                                                gTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTa                     <  2:1517682/150‑1 (MQ=255)
                                                                                                                                  cGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTAtt                   <  2:582876/150‑1 (MQ=255)
                                                                                                                                    gtcgagtggtTGATGACCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTATTGc                 >  2:1211136/6‑150 (MQ=255)
                                                                                                                                     tcgagtggtTGATGACCACCCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTATTGCt                >  1:821771/5‑150 (MQ=255)
                                                                                                                                     gCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTATTGCt                <  2:481157/150‑1 (MQ=255)
                                                                                                                                       gggTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTATTGCTTc              <  2:1151612/150‑1 (MQ=255)
                                                                                                                                       gggTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTATTGCTTCa             <  2:1195011/150‑1 (MQ=255)
                                                                                                                                         gtggtTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTATTGCTTCAc            <  1:605676/150‑1 (MQ=255)
                                                                                                                                           ggttggtCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTATTGCTTCACCg          >  2:1514474/1‑150 (MQ=255)
                                                                                                                                             ttggtCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTATTGCTTCACCGca        <  1:1146232/150‑1 (MQ=255)
                                                                                                                                                 tcatcaACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTATTGCTTCACCGCACACt    <  2:647005/150‑1 (MQ=255)
                                                                                                                                                   atcaACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTATTGCTTCACCGCACACTgg  >  1:883148/1‑150 (MQ=255)
                                                                                                                                                   atcaACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTATTGCTTCACCGCACACTgg  >  2:116828/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAGATGAATGGCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTATTGCTTCACCGCACACTGG  >  NC_000913/2399732‑2400028

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: