Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 346,135 T→C 5.9% N572D (AAC→GAC)  PP_0287 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947346,1350TC 5.9% 188.5 / 2.0 68N572D (AAC→GAC) PP_0287hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/32);  new base C (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-01

TACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTTCCACCGGCACCCGCTGGA  >  NC_002947/345986‑346281
                                                                                                                                                     |                                                                                                                                                  
tACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGtt                                                                                                                                                    >  2:917936/1‑150 (MQ=255)
tACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGtt                                                                                                                                                    >  2:719374/1‑150 (MQ=255)
   caccaGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTgc                                                                                                                                                 >  1:882714/1‑150 (MQ=255)
   caccaGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTgc                                                                                                                                                 >  1:883002/1‑150 (MQ=255)
   caccaGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTgc                                                                                                                                                 >  1:422377/1‑150 (MQ=255)
    accaGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTAACGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTgcg                                                                                                                                                >  2:86502/1‑150 (MQ=255)
      caGGCTGCCGCGCAGCTCCTGGAAGGGTATTTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCt                                                                                                                                              <  1:840983/150‑1 (MQ=255)
       aGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTg                                                                                                                                             <  2:422377/150‑1 (MQ=255)
           tGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGt                                                                                                                                         >  2:929408/1‑150 (MQ=255)
            gCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGtt                                                                                                                                        <  1:514295/150‑1 (MQ=255)
               gcgcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAg                                                                                                                                     <  1:819096/150‑1 (MQ=255)
                     cTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTc                                                                                                                               <  2:620985/150‑1 (MQ=255)
                      tCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCg                                                                                                                              >  1:44892/1‑150 (MQ=255)
                         tGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGAcc                                                                                                                           >  2:936565/1‑150 (MQ=255)
                         tGGAAGGGGGTGTATCCACCGCGTGGTTCGCTGCTCAGAGTGTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTTACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGAcc                                                                                                                           >  1:38396/1‑150 (MQ=255)
                          ggAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCa                                                                                                                          >  1:100537/1‑150 (MQ=255)
                                gggTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCtt                                                                                                                    >  1:926344/1‑150 (MQ=255)
                                     tCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                                    <  1:481153/129‑1 (MQ=255)
                                     tCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                                    >  2:481153/1‑129 (MQ=255)
                                       cTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGCTTGAGCGTGGCTCTGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGACGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCt                                                                                                             <  2:738260/150‑1 (MQ=255)
                                       cTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCt                                                                                                             <  1:871584/150‑1 (MQ=255)
                                       cTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTCTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCt                                                                                                             >  2:833106/1‑150 (MQ=255)
                                       cTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTCTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCt                                                                                                             >  2:833077/1‑150 (MQ=255)
                                        ttGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTg                                                                                                            >  1:637499/1‑150 (MQ=255)
                                        ttGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTg                                                                                                            >  1:401269/1‑150 (MQ=255)
                                            cgcgTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTAc                                                                                                        <  2:97781/150‑1 (MQ=255)
                                              cgTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCg                                                                                                      <  1:917936/150‑1 (MQ=255)
                                                  gTTCGCTTCTCATCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTgg                                                                                                  >  2:856558/1‑150 (MQ=255)
                                                  gTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGttg                                                                                                  <  2:401269/150‑3 (MQ=255)
                                                   ttCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACGGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGc                                                                                                 >  1:256547/1‑150 (MQ=255)
                                                     cGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTg                                                                                               >  1:206601/1‑150 (MQ=255)
                                                       ctgctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGAGCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTc                                                                                             <  2:605786/150‑1 (MQ=255)
                                                            cAGCGTCTTGCGATTGAGCGTGGCGATGACCGTGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATc                                                                                        >  1:869990/1‑150 (MQ=255)
                                                                    tGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtca                                                                                >  1:304879/1‑150 (MQ=255)
                                                                     gCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtcag                                                                               <  1:86502/150‑1 (MQ=255)
                                                                     gCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtcag                                                                               >  2:867522/1‑150 (MQ=255)
                                                                     gCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtcag                                                                               <  1:929408/150‑1 (MQ=255)
                                                                     gCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtcag                                                                               <  2:206601/150‑1 (MQ=255)
                                                                        aTTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtcagggt                                                                            <  1:936565/150‑1 (MQ=255)
                                                                          tGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGTGGGTCTACCAGTGCTGTCTGGCTGTTACCGGTGGCTGTCGGATCACTggtcagggtca                                                                          <  1:981985/150‑1 (MQ=255)
                                                                                     aTGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCNTTg                                                               <  2:699136/150‑1 (MQ=255)
                                                                                         cTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACcg                                                           >  1:692009/1‑150 (MQ=255)
                                                                                          ttGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCa                                                                                            >  2:344302/1‑116 (MQ=255)
                                                                                          ttGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCa                                                                                            <  1:344302/116‑1 (MQ=255)
                                                                                             gCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACcggcgg                                                       >  2:2581/1‑150 (MQ=255)
                                                                                              cacaACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCggcggc                                                      >  1:441014/1‑150 (MQ=255)
                                                                                                   ctgctgTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTCTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTg                                                 <  1:833077/150‑1 (MQ=255)
                                                                                                   ctgctgTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTCTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTg                                                 <  1:833106/150‑1 (MQ=255)
                                                                                                     gctgTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGACAAGTTAGTGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGtt                                               >  1:193242/1‑150 (MQ=255)
                                                                                                      ctgTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGATAATAGCGACTTCGATCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTc                                              <  2:869990/150‑1 (MQ=255)
                                                                                                       tgTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTcc                                             <  1:867522/150‑1 (MQ=255)
                                                                                                         ttCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCgg                                                                                                     <  1:1118025/92‑1 (MQ=255)
                                                                                                         ttCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCgg                                                                                                     >  2:1118025/1‑92 (MQ=255)
                                                                                                         ttCCAGGTTCGAATTGACCAGGACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCtt                                           >  2:325706/1‑150 (MQ=255)
                                                                                                             aGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACcggcgg                                                       <  2:549636/134‑1 (MQ=255)
                                                                                                             aGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACcggcgg                                                       >  1:549636/1‑134 (MQ=255)
                                                                                                                 tGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCgg                                   <  2:304879/150‑1 (MQ=255)
                                                                                                                    cATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTg                                          <  2:72875/140‑1 (MQ=255)
                                                                                                                    cATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTg                                          >  1:72875/1‑140 (MQ=255)
                                                                                                                     aTTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCtt                               <  2:882714/150‑1 (MQ=255)
                                                                                                                     aTTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCtt                               <  2:883002/150‑1 (MQ=255)
                                                                                                                       tGAGCAGCGCGCGATCGTTGATGGCGAAGTTTGCGCGGCCGGTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTCGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTtga                             <  2:1285/150‑1 (MQ=255)
                                                                                                                        gACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTtgat                            <  2:1098721/150‑1 (MQ=255)
                                                                                                                             gCACGCCAACGTTGATGGCGAAGTTTGCGTTGCCGTTCAGGGTGTCGACCCGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAgt                       <  2:608285/150‑1 (MQ=255)
                                                                                                                                     tCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTTCCAcc               >  1:606383/1‑150 (MQ=255)
                                                                                                                                             ggCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTTCCACCGGCACCCg       >  2:802305/1‑150 (MQ=255)
                                                                                                                                              gCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGAt                                                                                         <  1:226976/67‑1 (MQ=255)
                                                                                                                                              gCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGAt                                                                                         >  2:226976/1‑67 (MQ=255)
                                                                                                                                              gCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCGGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTTCCACCGGCACCCGc      >  1:985619/1‑150 (MQ=255)
                                                                                                                                                  aGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTTCCACCGGCACCCGCTGGa  <  1:802305/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTTCCACCGGCACCCGCTGGA  >  NC_002947/345986‑346281

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