Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 888,802 (A)7→6 coding (241/1686 nt) b0847 ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913888,7960A.96.8% 815.5 / 18.5 190coding (247/1686 nt)b0847putative transport protein
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base . (112/70);  total (116/74)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 221 total aligned reads displayed.

CAGCTTACCTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTAATGCTGAAATGT  >  NC_000913/888704‑888884
                                                                                            |                                                                                        
cAGCTTACCTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGa                                                                                          <  2:4618847/93‑1 (MQ=255)
 aGCTTACCTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaa                                                                                         <  2:1503581/93‑1 (MQ=255)
 aGCTTACCTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaa                                                                                         >  1:3295569/1‑93 (MQ=255)
  gCTTACCTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaaa                                                                                        >  1:792613/1‑93 (MQ=255)
   cTTACCTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaaaa                                                                                       <  2:2510104/93‑1 (MQ=255)
   cTTACCTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaaaa                                                                                       >  1:1880490/1‑93 (MQ=255)
     tACCTAACCCTACGGCGATCACCAGCGCACTGCCCACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaaacaa                                                                                     >  1:1978678/1‑93 (MQ=255)
       ccTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaaaaaagt                                                                                   <  2:2182251/93‑3 (MQ=255)
        cTACCCCTAAGGCGATCACCAGCGCGCTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaaaaaagtt                                                                                  <  1:1705179/93‑4 (MQ=255)
         tAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaaaaaagttc                                                                                 <  1:3679613/93‑5 (MQ=255)
            cccTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAAGTTCgg                                                                              <  2:4116558/93‑1 (MQ=255)
            cccTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGt                                                                             >  2:3227518/1‑93 (MQ=255)
             ccTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAAGTTCGGt                                                                             >  2:3906110/1‑93 (MQ=255)
              cTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTcc                                                                           >  2:3944563/1‑93 (MQ=255)
              cTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTcc                                                                           >  2:1526791/1‑93 (MQ=255)
              cTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTcc                                                                           >  1:4678772/1‑93 (MQ=255)
              cTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGGGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTcg                                                                           >  2:1175848/1‑92 (MQ=255)
               tAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCg                                                                          >  1:2027715/1‑93 (MQ=255)
                 aGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGc                                                                        <  1:2284660/93‑1 (MQ=255)
                  ggCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAAGTTCGGTCCGGc                                                                        <  1:2737428/93‑1 (MQ=255)
                  ggCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCt                                                                       <  2:2980832/93‑1 (MQ=255)
                  ggCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCt                                                                       <  1:1388682/93‑1 (MQ=255)
                  ggCGATCACCAGCGCACTGCCAACCATCACCAGTGCGAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCt                                                                       <  1:2889929/93‑1 (MQ=255)
                   gCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCtt                                                                      >  1:3062474/1‑93 (MQ=255)
                      aTCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGa                                                                   >  1:1950048/1‑93 (MQ=255)
                       tCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGAc                                                                  <  2:300923/93‑1 (MQ=255)
                        cACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGAcc                                                                 >  1:3924936/1‑93 (MQ=255)
                         aCCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGAccc                                                                <  2:3579740/93‑1 (MQ=255)
                         aCCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTGCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGAccc                                                                <  1:4458746/93‑1 (MQ=255)
                            aGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTCTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGa                                                             <  1:2633340/93‑1 (MQ=255)
                             gcgcACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAATTTGGTCCGGCTTCGAGCCCCAc                                                            >  1:3899882/1‑93 (MQ=255)
                             gcgcACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGAc                                                            <  2:1064369/93‑1 (MQ=255)
                              cgcACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACg                                                           >  2:254065/1‑93 (MQ=255)
                               gcACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGc                                                          >  2:2751334/1‑93 (MQ=255)
                               gcACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGc                                                          >  1:2160380/1‑93 (MQ=255)
                                cACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGATCGGTCCGGCTTCGACCCCGACGCa                                                         >  2:4199374/1‑93 (MQ=255)
                                 gcTGCCAACCATCACCAGTGGTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCa                                                         <  2:3824469/92‑1 (MQ=255)
                                 aCTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAg                                                        <  2:1436408/93‑1 (MQ=255)
                                  cTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGa                                                       >  2:667336/1‑93 (MQ=255)
                                  cTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGa                                                       >  1:855701/1‑93 (MQ=255)
                                    gCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCACGACGCAGaaa                                                     >  1:3085132/1‑93 (MQ=255)
                                      cAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGaaaaa                                                   >  2:1577735/1‑92 (MQ=255)
                                       gaCCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAACGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATa                                                  >  2:2881443/2‑93 (MQ=255)
                                       aaCCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATa                                                  >  1:1096014/1‑93 (MQ=255)
                                       aaCCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATa                                                  >  1:715231/1‑93 (MQ=255)
                                       aaCCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATa                                                  >  1:743449/1‑93 (MQ=255)
                                       aaCCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATa                                                  <  2:925838/93‑1 (MQ=255)
                                         ccATCGCCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATaaa                                                <  1:224276/93‑1 (MQ=255)
                                         ccATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATaaa                                                <  1:3372396/93‑1 (MQ=255)
                                          cATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAAc                                               <  2:3283933/93‑1 (MQ=255)
                                          cATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAAc                                               >  2:1708636/1‑93 (MQ=255)
                                          cATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAAc                                               >  1:3772767/1‑93 (MQ=255)
                                           aTCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACa                                              >  1:1209841/1‑93 (MQ=255)
                                            tCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAg                                             >  1:1587100/1‑93 (MQ=255)
                                            tCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAg                                             <  1:1445322/93‑1 (MQ=255)
                                             cACCAGTGCTAACAATAGGTTATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCCGTCCGGCTTCGACCCCGACGCAGAGAATAAACAGc                                            >  2:862628/1‑93 (MQ=255)
                                              cccAGTTCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCa                                           <  2:344218/92‑1 (MQ=255)
                                              aCCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCa                                           >  2:887644/1‑93 (MQ=255)
                                               ccAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCAt                                          >  2:96800/1‑93 (MQ=255)
                                                  gTGCTAACATTAGGTAATTTTTCTCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATaaa                                       <  2:3391867/93‑1 (MQ=255)
                                                  gTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATaaa                                       >  2:2935178/1‑93 (MQ=255)
                                                   tGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAg                                      <  1:4087658/93‑1 (MQ=255)
                                                    gCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGc                                     >  2:603554/1‑93 (MQ=255)
                                                     cTAACATTAGGTAATTTTTCCCATCGCGACAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGCAAATAAACAGCATAAAgac                                    >  2:1878225/1‑91 (MQ=255)
                                                      tAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCa                                   >  1:2917651/1‑93 (MQ=255)
                                                      tAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCa                                   >  1:318459/1‑93 (MQ=255)
                                                       aaCATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCaa                                  >  2:1312629/1‑93 (MQ=255)
                                                       aaCATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCaa                                  >  1:2947185/1‑93 (MQ=255)
                                                        aCATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAg                                 >  2:1136034/1‑93 (MQ=255)
                                                        aCATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAg                                 >  2:1964016/1‑93 (MQ=255)
                                                          aTTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGAt                               >  1:15361/1‑93 (MQ=255)
                                                           ttAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGAtt                              >  1:32445/1‑93 (MQ=255)
                                                             aGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTaa                            <  1:839822/93‑1 (MQ=255)
                                                             aGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTaa                            >  2:1292804/1‑93 (MQ=255)
                                                              ggTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAg                           <  1:396124/93‑1 (MQ=255)
                                                               gTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAgc                          >  1:902107/1‑93 (MQ=255)
                                                               gTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAgc                          <  2:2360622/93‑1 (MQ=255)
                                                               gTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAgc                          >  2:2321029/1‑93 (MQ=255)
                                                                  aTTTTTCCCATCGCGAAAAAAAATGGAAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAgcgc                        >  2:1992761/1‑93 (MQ=255)
                                                                   tttttCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCAt                      >  2:766691/1‑93 (MQ=255)
                                                                     tttCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCg                    >  2:4435240/1‑93 (MQ=255)
                                                                      ttCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCgg                   <  1:2811841/93‑1 (MQ=255)
                                                                       tCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTACGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGgt                  >  1:4236511/1‑93 (MQ=255)
                                                                        cccATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGgtg                 >  2:4242594/1‑93 (MQ=255)
                                                                        cccATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGgtg                 <  2:498262/93‑1 (MQ=255)
                                                                         ccATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGgtgt                <  2:3168632/93‑1 (MQ=255)
                                                                         ccATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGgtgt                <  2:1013619/93‑1 (MQ=255)
                                                                         ccATCCCGTGAAAAAATGGAAAAAAGATCGGTCCGGCTTCGAGCCCGACGCAGGAAATAAACATCCTAAAACCAACATTAAGCACATCTgtag                >  2:1907333/1‑91 (MQ=255)
                                                                          cATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGtt               <  1:2347627/93‑1 (MQ=255)
                                                                            tCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTaa             >  2:1362051/1‑93 (MQ=255)
                                                                            tCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTaa             >  1:907354/1‑93 (MQ=255)
                                                                            tCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTaa             <  1:3521419/93‑1 (MQ=255)
                                                                            tCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTaa             >  1:633886/1‑93 (MQ=255)
                                                                                gAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTAATGCt         >  1:4554675/1‑93 (MQ=255)
                                                                                 aaaaaaaaTGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTAATGCTg        >  1:4745147/1‑93 (MQ=255)
                                                                                   aaaaaaTGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTAATGCTGaa      >  1:316221/1‑93 (MQ=255)
                                                                                    aaaaaTGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTAATGCTGaaa     >  2:1073224/1‑93 (MQ=255)
                                                                                      aaaTGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTAATGCTGAAAtg   >  2:4530924/1‑93 (MQ=255)
                                                                                      aaaTGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTAATGCTGAAAtg   >  1:4160389/1‑93 (MQ=255)
                                                                                       aaTGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTAATGCTGAAAtgt  >  1:4383466/1‑93 (MQ=255)
                                                                                            |                                                                                        
CAGCTTACCTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTAATGCTGAAATGT  >  NC_000913/888704‑888884

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