Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,410,196 T→A 24.5% intergenic (+245/‑84) prmA → / → dusB methyltransferase for 50S ribosomal subunit protein L11/tRNA‑dihydrouridine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,410,1960TA24.5% 323.9 / 108.7 167intergenic (+245/‑84)prmA/dusBmethyltransferase for 50S ribosomal subunit protein L11/tRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base T (45/81);  new base A (25/16);  total (70/97)
Fisher's exact test for biased strand distribution p-value = 6.06e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 172 total aligned reads displayed.

TGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAA  >  NC_000913/3410097‑3410294
                                                                                                    |                                                                                                  
tgtAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGAtt                                                                                                    <  2:2201636/100‑1 (MQ=255)
 gtAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGAtt                                                                                                    >  1:3078425/1‑99 (MQ=255)
  tAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTgg                                                                                                  <  2:601554/100‑1 (MQ=255)
  tAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTgg                                                                                                  <  1:3811208/100‑1 (MQ=255)
   aaaTTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGt                                                                                                 <  2:1972685/100‑1 (MQ=255)
        ttAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCACTTCTTTGATCCATCTCAGAGGATTGGTCAAAg                                                                                            <  2:4049331/100‑1 (MQ=255)
        ttAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAg                                                                                            <  1:1710634/100‑1 (MQ=255)
           gTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATAGGTCAAAGttt                                                                                         >  1:570177/1‑100 (MQ=255)
            tCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTg                                                                                        <  2:104731/100‑1 (MQ=255)
            tCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATAGGTCAAAGTTTg                                                                                        <  1:413771/100‑1 (MQ=255)
             cACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTgg                                                                                       >  1:1291001/1‑100 (MQ=255)
             cACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTgg                                                                                       <  1:3927172/100‑1 (MQ=255)
              aCTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGc                                                                                      <  1:3760028/100‑1 (MQ=255)
              aCTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGc                                                                                      <  1:6577/100‑1 (MQ=255)
              aCTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGc                                                                                      >  2:1618915/1‑100 (MQ=255)
                  ttttCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAgatcggaagag                                                                                  >  1:4111715/1‑90 (MQ=255)
                    ttCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCa                                                                                <  1:1995433/100‑1 (MQ=255)
                    ttCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCa                                                                                <  2:3863287/100‑1 (MQ=255)
                      cTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCAtc                                                                              >  2:3222724/1‑100 (MQ=255)
                          tatGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGt                                                                          <  1:654204/100‑1 (MQ=255)
                           atgatgATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGtc                                                                         >  2:3824186/1‑99 (MQ=255)
                           atgatgATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTg                                                                         <  2:1957902/100‑1 (MQ=255)
                            tgatgaTTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGc                                                                        <  2:1706902/100‑1 (MQ=255)
                              atgatTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCCGAGGATTGGTCAAAGTTTGCCCTTTCATCTCGTGCaa                                                                      >  2:47747/1‑100 (MQ=255)
                              atgatTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCaa                                                                      >  2:360337/1‑100 (MQ=255)
                               tgatTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCaaa                                                                     >  2:1977583/1‑100 (MQ=255)
                                gatTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCaaaa                                                                    <  2:3572379/100‑1 (MQ=255)
                                      tATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCg                                                              <  2:1273292/100‑1 (MQ=255)
                                      tATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCg                                                              <  2:256921/100‑1 (MQ=255)
                                       aTGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGt                                                             >  1:487137/1‑100 (MQ=255)
                                       aTGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGt                                                             >  2:806482/1‑100 (MQ=255)
                                        tGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTa                                                            <  2:413043/100‑1 (MQ=255)
                                         gAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGTGGATTCGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTa                                                            >  2:1500054/1‑99 (MQ=255)
                                         gAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTaa                                                           >  2:2960784/1‑100 (MQ=255)
                                         gAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTaa                                                           >  1:108436/1‑100 (MQ=255)
                                          aGTAATTATCGCACCACGCTCA‑TTTTAAAGGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGCTTGGCCTTTCATCTCGTGCAAAAAATGCGTAat                                                          <  1:2568975/100‑1 (MQ=255)
                                          aGTAATTATCAGCCCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGACGATTGGCCAAAGTTTGGCCTTTCATCTCGTGCAAACAAt                                                                 <  1:2683821/93‑1 (MQ=255)
                                           gTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAata                                                         >  1:3236251/1‑100 (MQ=255)
                                            tAATTATCGCACCACGCTCA‑TTTTAAATGAAATTCTTTGATACATCTCAGAGGATAGGTCAAAGTTTGGCATTACAAATAGAAAAAAAAATGCGTAata                                                         <  1:1575833/99‑1 (MQ=255)
                                             aactATCCCACCACGCTCAGTTTAAAATGCAATGCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAatat                                                        <  2:1419524/97‑1 (MQ=255)
                                              aTTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATAc                                                      >  2:2018382/1‑100 (MQ=255)
                                              aTTAGCGCAACACGGGGA‑TTTTAAAAGAAATTCTTTGATCCAACTCAGAGGATAGGTCAAAGTTTGGCCTTTAAAAAAAAGAAAAAAATGCGTAAtata                                                       <  1:2117804/99‑1 (MQ=255)
                                               ttATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATAcg                                                     >  2:117248/1‑100 (MQ=255)
                                               ttATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATAcg                                                     <  2:3997911/100‑1 (MQ=255)
                                                 aTCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATAcgcc                                                   <  2:4045379/100‑1 (MQ=255)
                                                 aTCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATAcgcc                                                   <  1:280204/100‑1 (MQ=255)
                                                 aTCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATAcgc                                                    >  1:3513812/1‑99 (MQ=255)
                                                 aTCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATAcgcc                                                   >  1:3300053/1‑100 (MQ=255)
                                                 aTCGCACCACGATCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCAACAAAAAAAAAAAATGCGTAATATAcgcc                                                   <  1:1855644/100‑1 (MQ=255)
                                                     caccacGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCt                                               <  2:2864824/100‑1 (MQ=255)
                                                     caccacGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCt                                               >  1:1237013/1‑100 (MQ=255)
                                                        cacGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGc                                            <  2:3497701/100‑1 (MQ=255)
                                                        cacGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGc                                            <  2:1683849/100‑1 (MQ=255)
                                                        cacGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGc                                            <  1:1941013/100‑1 (MQ=255)
                                                        cacGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGc                                            <  2:3210791/100‑1 (MQ=255)
                                                         acGCTCA‑TTTTAAATGCAATTCTTTGATTCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGc                                            <  1:2334936/99‑1 (MQ=255)
                                                          cGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAg                                          >  1:4109626/1‑100 (MQ=255)
                                                           gCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGt                                         <  2:3913519/100‑1 (MQ=255)
                                                           gCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGt                                         <  2:2613724/100‑1 (MQ=255)
                                                           gCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGt                                         <  2:2934262/100‑1 (MQ=255)
                                                              cA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTcac                                      <  2:307407/100‑1 (MQ=255)
                                                              cA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTcac                                      >  2:3458000/1‑100 (MQ=255)
                                                              cA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTcac                                      >  2:1531037/1‑100 (MQ=255)
                                                              cA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTcac                                      <  2:2320244/100‑1 (MQ=255)
                                                               a‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTcaca                                     >  1:2614695/1‑100 (MQ=255)
                                                                 ttttAAATGCAATTCTTTGATCCATCTCAGAGGATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAg                                    <  2:2953020/100‑1 (MQ=255)
                                                                   ttAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTa                                  <  2:3664861/100‑1 (MQ=255)
                                                                    tAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTAt                                 <  1:504086/100‑1 (MQ=255)
                                                                     aaaTGCAATTCTTCGATACAACTTAGAGTATTGGTCAAAGTTTGGCCTTACAAAAAAAAAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATg                                <  1:4075457/100‑1 (MQ=255)
                                                                       aTGCAATTCTTTGATCCATCTCAGAGGATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGt                              >  2:1240208/1‑100 (MQ=255)
                                                                        tGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTc                             <  1:2421415/100‑1 (MQ=255)
                                                                         gccaTTCTTTGATACATCTCACATGTTTAGTCAAAGTTTGGCCGTTCATCTCGTGCAAAAAATGCTTAATATACGTCGTCTTGCAGTCACAGTATGGt                              <  1:2484655/95‑1 (MQ=255)
                                                                         gCAATTCTTTGATCCATCTCAGAGGATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCa                            >  1:755593/1‑100 (MQ=255)
                                                                          cAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGAGTAATATACGCCGCCTTGCAGTCACAGTATgacc                             >  1:850929/1‑95 (MQ=255)
                                                                           aaTTCTTTGATCCATCTCAGAGGATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCAtt                          >  1:2943762/1‑100 (MQ=255)
                                                                              tCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCt                       <  1:3853335/100‑1 (MQ=255)
                                                                               cTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTcct                      >  1:2595466/1‑98 (MQ=255)
                                                                               cTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCtt                      <  1:1696574/100‑1 (MQ=255)
                                                                                tttGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTa                     <  2:2459259/100‑1 (MQ=255)
                                                                                tttGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTa                     <  1:48660/100‑1 (MQ=255)
                                                                                tttGATCCATCTCAAAGGATTAGTCAAAGTTTGGCCTTAAAAAAAAAAAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTc                        <  1:1618911/97‑1 (MQ=255)
                                                                                  tGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAAc                   <  2:3676782/100‑1 (MQ=255)
                                                                                  tGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAAc                   <  1:1007699/100‑1 (MQ=255)
                                                                                     tCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCa                <  2:3650163/100‑1 (MQ=255)
                                                                                     tCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCa                >  2:510852/1‑100 (MQ=255)
                                                                                      ccATCTCAGAGGATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCAt               >  2:1007820/1‑100 (MQ=255)
                                                                                       cATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATg              >  1:3472343/1‑100 (MQ=255)
                                                                                       cATCTCAGAGGATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATg              >  2:3977433/1‑100 (MQ=255)
                                                                                         tctcAGAGGATTGGTCAAAGTTTGGCCTTTCATCTTGTGCAAAAAATGCGTAATATACGGCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATgcg            <  2:2394173/100‑1 (MQ=255)
                                                                                         tctcAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATgcg            <  1:2001418/100‑1 (MQ=255)
                                                                                         tctcAGAGGATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATgcg            <  1:2018378/100‑1 (MQ=255)
                                                                                           tcAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCa          >  2:1348997/1‑100 (MQ=255)
                                                                                            cAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCAt         <  1:1734581/100‑1 (MQ=255)
                                                                                            cAGAGGATTGGTCAAAGTTTGGACTTTAAAAAAAAGAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATgcg            <  1:1731672/97‑1 (MQ=255)
                                                                                            cAGAGGATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGAAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCAt         >  1:2033013/1‑100 (MQ=255)
                                                                                              gagGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCCTAACTCATGCGCATCg       >  1:2495582/1‑100 (MQ=255)
                                                                                              gagGATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCg       >  2:3748735/1‑100 (MQ=255)
                                                                                               agGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCgg      <  2:2258228/100‑1 (MQ=255)
                                                                                                ggATAGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCgg      >  1:1305941/1‑99 (MQ=255)
                                                                                                  aTTGGTCAAAGTTTGGCCTTTAAAATAAAAAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACa   <  1:1532960/100‑1 (MQ=255)
                                                                                                    |                                                                                                  
TGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCA‑TTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAA  >  NC_000913/3410097‑3410294

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