Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,062,795 T→A 5.1% I556N (ATT→AAT)  PP_1838 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,062,7950TA 5.1% 246.7 / 3.1 78I556N (ATT→AAT) PP_1838membrane protein
Reads supporting (aligned to +/- strand):  ref base T (37/37);  new base A (2/2);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

AGTTCTTCGAGAAGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCA  >  NC_002947/2062648‑2062944
                                                                                                                                                   |                                                                                                                                                     
aGTTCTTCGAGAAGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTc                                                                                                                                                     >  1:1036112/1‑150 (MQ=255)
  ttcttcGAGAAGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAg                                                                                                                                                   <  2:1516239/150‑1 (MQ=255)
    cttcGAGAAGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAg                                                                                                                                                   >  1:815059/1‑148 (MQ=255)
    cttcGAGAAGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAg                                                                                                                                                   <  2:815059/148‑1 (MQ=255)
        gagaAGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAg                                                                                                                                             >  2:1112160/1‑150 (MQ=255)
              gcgcGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTcggt                                                                                                                                       <  1:1112160/150‑1 (MQ=255)
                  gTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTcggtgcgg                                                                                                                                   >  2:150243/1‑150 (MQ=255)
                  gTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTcggtgcgg                                                                                                                                   >  2:693320/1‑150 (MQ=255)
                     aGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGa                                                                                                                                <  1:585919/150‑1 (MQ=255)
                      ggAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGaa                                                                                                                               <  1:26388/150‑1 (MQ=255)
                         gCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTATCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATca                                                                                                                            >  2:1615784/1‑150 (MQ=255)
                            ttACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCacac                                                                                                                         >  1:71821/1‑150 (MQ=255)
                                ttCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTc                                                                                                                                                     <  2:461256/118‑1 (MQ=255)
                                ttCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTc                                                                                                                                                     >  1:461256/1‑118 (MQ=255)
                                ttCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTcggtgcgg                                                                                                                                   <  2:511854/136‑1 (MQ=255)
                                ttCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTcggtgcgg                                                                                                                                   >  1:511854/1‑136 (MQ=255)
                                ttCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTg                                                                                                                     <  2:307868/150‑1 (MQ=255)
                                ttCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTg                                                                                                                     <  2:307923/150‑1 (MQ=255)
                                     ggACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGga                                                                                                                                                 >  1:811299/1‑117 (MQ=255)
                                     ggACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGga                                                                                                                                                 <  2:811299/117‑1 (MQ=255)
                                         acCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAg                                                                                                            >  1:1401630/1‑150 (MQ=255)
                                         acCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGAATGCCCGGGGCAATCAg                                                                                                                                                   >  2:1294924/1‑111 (MQ=255)
                                         acCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGAATGCCCGGGGCAATCAg                                                                                                                                                   <  1:1294924/111‑1 (MQ=255)
                                           ccTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGt                                                                                                          >  1:143371/1‑150 (MQ=255)
                                              gTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGa                                                                                                       <  2:71821/150‑1 (MQ=255)
                                                  gTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGCGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATc                                                                                                   >  2:623115/1‑150 (MQ=255)
                                                  gTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTGCGTTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATc                                                                                                   >  1:370216/1‑150 (MQ=255)
                                                   tCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCg                                                                                                  >  2:86262/1‑150 (MQ=255)
                                                    cATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGt                                                                                                 >  1:526587/1‑150 (MQ=255)
                                                     aTCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGAGCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGGGCGGTGAATCACACCGTGGGTATCGAGGTCGACATCGTg                                                                                                >  2:621310/1‑150 (MQ=255)
                                                     aTCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGaa                                                                                                                                              >  1:475121/1‑104 (MQ=255)
                                                     aTCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGaa                                                                                                                                              <  2:475121/104‑1 (MQ=255)
                                                      tcgtcgGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGc                                                                                               >  2:993712/1‑150 (MQ=255)
                                                         tcgGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGAATGCCCGGGGCAATCa                                                                                                                                                    <  2:40699/94‑1 (MQ=255)
                                                         tcgGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGAATGCCCGGGGCAATCa                                                                                                                                                    >  1:40699/1‑94 (MQ=255)
                                                          cgGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGAc                                                                                           <  1:704823/150‑1 (MQ=255)
                                                             tGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGACTCACACCGTGGGTACCCAGGTCGACATCGTGCCGACcat                                                                                        >  1:681601/1‑150 (MQ=255)
                                                             tGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACcat                                                                                        >  1:1395760/1‑150 (MQ=255)
                                                              gACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACcatc                                                                                       <  2:1092444/150‑1 (MQ=255)
                                                                 cATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATg                                                                                    >  1:741292/1‑150 (MQ=255)
                                                                 cATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATg                                                                                    >  2:1146975/1‑150 (MQ=255)
                                                                    ggCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGgc                                                                                 <  1:966425/150‑1 (MQ=255)
                                                                      cTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGgccg                                                                               >  1:1423473/1‑150 (MQ=255)
                                                                       tttGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTc                                                                                                                                                     >  1:1248784/1‑79 (MQ=255)
                                                                       tttGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTc                                                                                                                                                     <  2:1248784/79‑1 (MQ=255)
                                                                             aaCCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTcg                                                                                                                                         >  1:786609/1‑85 (MQ=255)
                                                                             aaCCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTcg                                                                                                                                         <  2:786609/85‑1 (MQ=255)
                                                                             aaCCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGa                                                                                                                                >  2:73383/1‑94 (MQ=255)
                                                                             aaCCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGa                                                                                                                                <  1:73383/94‑1 (MQ=255)
                                                                               ccaccaGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTgg                                                                      >  1:1554003/1‑150 (MQ=255)
                                                                                   cagcagGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAg                                                                  >  1:805698/1‑150 (MQ=255)
                                                                                    agcagGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAgg                                                                 <  2:1395760/150‑1 (MQ=255)
                                                                                            aCCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGccac                                                         <  1:993712/150‑1 (MQ=255)
                                                                                                  cTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGc                                                   >  2:879974/1‑150 (MQ=255)
                                                                                                     gACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTgg                                                <  1:879974/150‑1 (MQ=255)
                                                                                                      aCCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCCCCAgcactggt                                               <  2:347462/150‑8 (MQ=255)
                                                                                                           ggCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGAGTACCTCACTCATTGTCTTATCTACCTTCATGGGTCGCCTGACTCCCAAGGCTCGCCACCAGTGCTGGGGgcgc                                          >  1:849164/1‑150 (MQ=255)
                                                                                                                cTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCt                                     <  2:454782/150‑1 (MQ=255)
                                                                                                                       gTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGccac                                                         <  1:1343080/123‑1 (MQ=255)
                                                                                                                       gTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGccac                                                         >  2:1343080/1‑123 (MQ=255)
                                                                                                                       gTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAAc                              <  2:741292/150‑1 (MQ=255)
                                                                                                                       gTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAAc                              <  2:1554003/150‑1 (MQ=255)
                                                                                                                          ccGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGa                                        <  2:1580890/137‑1 (MQ=255)
                                                                                                                          ccGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGa                                        >  1:1580890/1‑137 (MQ=255)
                                                                                                                                   ctgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGt                  <  2:143371/150‑1 (MQ=255)
                                                                                                                                   ctgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTAGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGt                  <  1:291767/150‑1 (MQ=255)
                                                                                                                                    tgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTg                 <  2:972528/150‑1 (MQ=255)
                                                                                                                                        tGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGAAAAGCTCGCGACCAGGGTGGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCa             <  1:693320/150‑1 (MQ=255)
                                                                                                                                         gCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAg            <  2:1401630/150‑1 (MQ=255)
                                                                                                                                          ccccGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAgg           <  2:526587/150‑1 (MQ=255)
                                                                                                                                           cccGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAggg          >  2:497029/1‑150 (MQ=255)
                                                                                                                                           cccGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAggg          >  1:1478442/1‑150 (MQ=255)
                                                                                                                                            ccGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGc         <  1:1615784/150‑1 (MQ=255)
                                                                                                                                            ccGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTCGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCCACCTGCTCAACCTGCCCGAGGGTGACCAGGGc         <  1:621310/150‑1 (MQ=255)
                                                                                                                                             cGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATgg                                                                                   <  1:1628393/75‑1 (MQ=255)
                                                                                                                                             cGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATgg                                                                                   >  2:1628393/1‑75 (MQ=255)
                                                                                                                                             cGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGc                          <  1:247275/132‑1 (MQ=255)
                                                                                                                                             cGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGc                          >  2:247275/1‑132 (MQ=255)
                                                                                                                                                   ttCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCCGGGCGTAGGCa  >  1:889679/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGTTCTTCGAGAAGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCA  >  NC_002947/2062648‑2062944

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: