Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 991,560 T→C 6.5% I60T (ATC→ACC)  der → GTPase Der

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947991,5600TC 6.5% 317.0 / 8.7 107I60T (ATC→ACC) derGTPase Der
Reads supporting (aligned to +/- strand):  ref base T (53/47);  new base C (3/4);  total (56/51)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.41e-01

Only 100 of 107 total aligned reads displayed.

CGTCGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCAAGCG  >  NC_002947/991417‑991707
                                                                                                                                               |                                                                                                                                                   
cGTCGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCgg                                                                                                                                               <  2:534587/150‑1 (MQ=255)
   cGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGa                                                                                                                                            <  2:1605417/150‑1 (MQ=255)
    ggCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGAc                                                                                                                                           >  2:360981/1‑150 (MQ=255)
    ggCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGATACCACCGGTgtc                                                                                                                                           <  1:1668811/150‑3 (MQ=255)
     gCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGGGACCGCCAGTATGGTGATGCCAGCAGGGATGGGCGGTCCTTCATCCTGATCGACACCGGTGGTCTCACCGGTGACg                                                                                                                                          >  1:1086581/1‑150 (MQ=255)
       aaaTCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGaa                                                                                                                                        <  1:360981/150‑1 (MQ=255)
         aTCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGt                                                                                                                                      <  2:1833595/150‑1 (MQ=255)
          tCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTg                                                                                                                                     >  1:1169971/1‑150 (MQ=255)
          tCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTg                                                                                                                                     >  2:104453/1‑150 (MQ=255)
             accaTGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGc                                                                                                                                  >  1:532516/1‑150 (MQ=255)
                   ttCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCAcc                                                                                                                                                 >  1:532114/1‑129 (MQ=255)
                   ttCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCAcc                                                                                                                                                 <  2:532114/129‑1 (MQ=255)
                   ttCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGAc                                                                                                                            >  1:1140207/1‑150 (MQ=255)
                   ttCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGAc                                                                                                                            <  1:104453/150‑1 (MQ=255)
                     cAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGa                                                                                                                          >  2:622157/1‑150 (MQ=255)
                       aCCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCAt                                                                                                                                >  2:1474191/1‑142 (MQ=255)
                           ccTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAAt                                                                                                                    <  1:1137935/150‑1 (MQ=255)
                                ccAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTggg                                                                                                                                   <  1:1747777/130‑1 (MQ=255)
                                ccAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTggg                                                                                                                                   >  2:1747777/1‑130 (MQ=255)
                                 cAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGa                                                                                                              <  1:867318/150‑1 (MQ=255)
                                  aaGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAg                                                                                                             <  2:1091427/150‑1 (MQ=255)
                                    gACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCa                                                                                                           >  1:1529646/1‑150 (MQ=255)
                                                cATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGTCACCGGTGATACCACCGGTGTCGatcag                                                                                                                                     <  1:524911/112‑5 (MQ=38)
                                                cATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGTCACCGGTGATACCACCGGTGTCGatcag                                                                                                                                     >  2:524911/1‑108 (MQ=38)
                                                cATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCg                                                                                                       >  2:1774980/1‑142 (MQ=255)
                                                cATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCg                                                                                                       <  1:1774980/142‑1 (MQ=255)
                                                cATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGc                                                                                               >  2:1353113/1‑150 (MQ=255)
                                                cATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGc                                                                                               >  1:1041425/1‑150 (MQ=255)
                                                cATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCCGTGACGAAGTGGGCATGGACGAAAAAATGGCCGCGCACTCGCTCATGGc                                                                                               >  2:548976/1‑150 (MQ=255)
                                                  tcgtcgGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCa                                                                                             >  1:616382/1‑150 (MQ=255)
                                                    gtcgGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCAtt                                                                                           <  2:1165627/150‑1 (MQ=255)
                                                     tcgGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATc                                                                                                                                                    >  2:1137425/1‑92 (MQ=255)
                                                     tcgGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATc                                                                                                                                                    <  1:1137425/92‑1 (MQ=255)
                                                       ggTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTgaa                                                                                        <  2:1140207/150‑1 (MQ=255)
                                                        gTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTgaag                                                                                       <  1:1023671/150‑1 (MQ=255)
                                                        gTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTAATCGACACCGGTGGTATCACCGGTGACGAAGTGAGAAAAAAAAAAAAAATGGCCGAGCAGTCGCTCATGGCCATTgaag                                                                                       <  2:489842/150‑1 (MQ=255)
                                                          gACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGc                                                                                                                 >  2:1529812/1‑122 (MQ=255)
                                                          gACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGc                                                                                                                 <  1:1529812/122‑1 (MQ=255)
                                                             cTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGcc                                                                                  >  2:114890/1‑150 (MQ=255)
                                                             cTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGACAGATGGTCGGGATGTCTCTTCATCCGGATCGATTCCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGcc                                                                                  >  1:36270/1‑150 (MQ=255)
                                                              tGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCg                                                                                 >  2:1318235/1‑150 (MQ=255)
                                                                     ccTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAtgt                                                                          >  2:1290224/1‑150 (MQ=255)
                                                                      cTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCAtt                                                                                           >  2:70036/1‑132 (MQ=255)
                                                                      cTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCAtt                                                                                           <  1:70036/132‑1 (MQ=255)
                                                                        gACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCt                                                                       >  2:895466/1‑150 (MQ=255)
                                                                        gACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCt                                                                       >  1:1115070/1‑150 (MQ=255)
                                                                         aCCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACCAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTg                                                                      >  1:1854283/1‑150 (MQ=255)
                                                                           ccGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGtt                                                                    >  2:648130/1‑150 (MQ=255)
                                                                           ccGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGtt                                                                    >  2:1220143/1‑150 (MQ=255)
                                                                            cGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTc                                                                   >  1:1256993/1‑150 (MQ=255)
                                                                               gACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTg                                                                <  1:201472/150‑1 (MQ=255)
                                                                                aCCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTgg                                                               <  2:36270/150‑1 (MQ=255)
                                                                                  cGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGtt                                                             <  2:1854283/150‑1 (MQ=255)
                                                                                  cGCCAGGATGGTTATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGtt                                                             <  1:298074/150‑1 (MQ=255)
                                                                                    ccAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTGGCTCATGGCCATTGAAGAAGACGACTATGTGCTGTTCCTGGTTGa                                                           >  1:782785/1‑150 (MQ=255)
                                                                                     cAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGAt                                                          >  1:802915/1‑150 (MQ=255)
                                                                                     cAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGAt                                                          <  2:532516/150‑1 (MQ=255)
                                                                                     cAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGAt                                                          >  1:1511043/1‑150 (MQ=255)
                                                                                       gTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGc                                                        >  1:1441364/1‑150 (MQ=255)
                                                                                        tATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGcc                                                       >  2:258004/1‑150 (MQ=255)
                                                                                        tATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGcc                                                       >  1:1671698/1‑150 (MQ=255)
                                                                                            gagATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTg                                                   >  2:1935701/3‑150 (MQ=255)
                                                                                            gTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTg                                                   >  2:645097/1‑150 (MQ=255)
                                                                                              gATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCg                                                 >  1:1625840/1‑150 (MQ=255)
                                                                                               aTGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCgg                                                <  2:1529646/150‑1 (MQ=255)
                                                                                               aTGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCgg                                                <  2:1041425/150‑1 (MQ=255)
                                                                                                 gCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCTGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGt                                              <  1:114890/150‑1 (MQ=255)
                                                                                                  ccAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTa                                             >  2:816841/1‑150 (MQ=255)
                                                                                                  ccAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTa                                             <  1:1318235/150‑1 (MQ=255)
                                                                                                   cAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTAt                                            <  1:197510/150‑1 (MQ=255)
                                                                                                   cAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTAt                                            <  2:1511043/150‑1 (MQ=255)
                                                                                                    aGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGc                                                                                   <  1:445056/110‑1 (MQ=255)
                                                                                                    aGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGc                                                                                   >  2:445056/1‑110 (MQ=255)
                                                                                                     gCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGa                                          <  2:1625840/150‑1 (MQ=255)
                                                                                                            gggTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGc                                   >  1:373142/1‑150 (MQ=255)
                                                                                                              gTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGGATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAg                                 >  1:1894708/1‑150 (MQ=255)
                                                                                                                cGTTCCTTCATCCTGATCGACCCCGGTGGTCTCACTGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGAc                               <  2:802915/150‑1 (MQ=255)
                                                                                                                 gTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGAt                                                          >  1:1918935/1‑122 (MQ=255)
                                                                                                                 gTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGAt                                                          <  2:1918935/122‑1 (MQ=255)
                                                                                                                     cTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGAGGCCCGTGCGGGTATGACTGCGGCAGACCAgat                          <  2:1115070/150‑1 (MQ=255)
                                                                                                                      ttCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTgg                                                               <  2:193286/112‑1 (MQ=255)
                                                                                                                      ttCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTgg                                                               >  1:193286/1‑112 (MQ=255)
                                                                                                                         aTCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCt                                                                 >  1:155703/1‑107 (MQ=255)
                                                                                                                         aTCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCt                                                                 <  2:155703/107‑1 (MQ=255)
                                                                                                                             tGATCGACACCGGTGGTATCACCGGTGACGAAGTGTGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTTCGGCAGACCAGATGATCGCCg                  <  1:548976/150‑1 (MQ=255)
                                                                                                                             tGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCTATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCg                  >  2:1730149/1‑150 (MQ=255)
                                                                                                                             tGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAt                                                                            <  1:684507/92‑1 (MQ=255)
                                                                                                                             tGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAt                                                                            >  2:684507/1‑92 (MQ=255)
                                                                                                                              gATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGa                 <  2:1671698/150‑1 (MQ=255)
                                                                                                                              gATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGa                 <  2:616382/150‑1 (MQ=255)
                                                                                                                              gATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGa                 <  2:1441364/150‑1 (MQ=255)
                                                                                                                                   tcacCGGTGATACCACCGGTGTCGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTAt                                            >  1:1656079/2‑118 (MQ=255)
                                                                                                                                   tcacCGGTGATACCACCGGTGTCGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTAt                                            <  2:1656079/117‑1 (MQ=255)
                                                                                                                                   tcacCGGTGATACCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATc                      >  1:517033/2‑140 (MQ=255)
                                                                                                                                    cacCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCt           >  1:540232/1‑150 (MQ=255)
                                                                                                                                     acCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTg          >  1:1463266/1‑150 (MQ=255)
                                                                                                                                      ccGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTgc         <  1:816841/150‑1 (MQ=255)
                                                                                                                                         gtggtATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAg                            >  1:923244/1‑128 (MQ=255)
                                                                                                                                            gtATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCAAGc   >  2:1819050/1‑150 (MQ=255)
                                                                                                                                             tATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCAAGCg  >  1:1127723/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CGTCGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCAAGCG  >  NC_002947/991417‑991707

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