Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,065,866 T→C intergenic (‑523/‑764) trs56 ← / → b1995 IS5 transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,065,8660TC99.3% 461.9 / ‑4.5 147intergenic (‑523/‑764)trs56/b1995IS5 transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (80/66);  total (81/66)
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.79e-01

Only 100 of 147 total aligned reads displayed.

CATTTTCAGAAGCTCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAG  >  NC_000913/2065777‑2065955
                                                                                         |                                                                                         
cATTTTCAGAAGCTCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACg                                                                                        >  1:424657/1‑93 (MQ=255)
  ttttCAGAAGCTCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGaa                                                                                      <  2:4010048/93‑1 (MQ=255)
  ttttCAGAAGCTCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGaa                                                                                      <  2:367368/93‑1 (MQ=255)
       aGAAGCTCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGc                                                                                 >  1:3013864/1‑93 (MQ=255)
       aGAAGCTCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGc                                                                                 >  2:4465161/1‑93 (MQ=255)
          aGCTCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCa                                                                              >  1:3072963/1‑93 (MQ=255)
          aGCTCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCa                                                                              <  1:3822919/93‑1 (MQ=255)
            ctctCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAgg                                                                            >  2:1921678/1‑93 (MQ=255)
             tctcCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGa                                                                           >  1:1190665/1‑93 (MQ=255)
             tctcCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGa                                                                           >  2:3938816/1‑93 (MQ=255)
             tctcCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGa                                                                           <  2:1243023/93‑1 (MQ=255)
              ctcCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGAt                                                                          <  2:709649/93‑1 (MQ=255)
              ctcCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGAt                                                                          <  2:2915616/93‑1 (MQ=255)
              ctcCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGAt                                                                          <  2:690091/93‑1 (MQ=255)
                ccAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCa                                                                        <  2:2708085/93‑1 (MQ=255)
                 cAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCAt                                                                       <  1:2053090/93‑1 (MQ=255)
                   aGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTa                                                                     >  2:2743959/1‑93 (MQ=255)
                   aGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTa                                                                     <  2:2971898/93‑1 (MQ=255)
                   aGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTa                                                                     >  2:1660220/1‑93 (MQ=255)
                   aGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTa                                                                     >  2:3641279/1‑93 (MQ=255)
                    gAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAg                                                                    <  1:1596539/93‑1 (MQ=255)
                     aaaGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAgt                                                                   >  2:4276151/1‑93 (MQ=255)
                     aaaGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAgt                                                                   <  1:1146060/93‑1 (MQ=255)
                     aaaGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAgt                                                                   <  1:1661996/93‑1 (MQ=255)
                       aGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAgtgg                                                                 >  1:373699/1‑93 (MQ=255)
                        gCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAgtggt                                                                >  2:2460788/1‑93 (MQ=255)
                        gCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAgtggt                                                                >  1:4059410/1‑93 (MQ=255)
                         cGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAgtggtg                                                               <  2:3045058/93‑1 (MQ=255)
                         cGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAgtggtg                                                               >  2:3012428/1‑93 (MQ=255)
                           aaTTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCt                                                             <  1:4388149/93‑1 (MQ=255)
                           aaTTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCt                                                             >  1:3985651/1‑93 (MQ=255)
                            aTTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTg                                                            >  2:1188285/1‑93 (MQ=255)
                            aTTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTg                                                            >  1:147954/1‑93 (MQ=255)
                             ttCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTga                                                           <  2:4152868/93‑1 (MQ=255)
                                ttttCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGaa                                                        <  2:1443581/93‑1 (MQ=255)
                                ttttCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGaa                                                        >  2:2206154/1‑93 (MQ=255)
                                ttttCTCAAATTCTGCATTCTCAGACACAGATTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGaa                                                        >  1:3177608/1‑93 (MQ=255)
                                 tttCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAg                                                       >  1:752348/1‑93 (MQ=255)
                                 tttCTCAAATTCTGCATTCTCAGACACAGCTTTCTACCCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAg                                                       >  2:279825/1‑93 (MQ=255)
                                  ttCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGAGGTGCTGAGAAGt                                                      <  1:2538144/93‑1 (MQ=255)
                                   tctcAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGtt                                                     >  1:3463550/1‑93 (MQ=255)
                                     tcAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTaa                                                   >  2:3931515/1‑93 (MQ=255)
                                      cAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTaaa                                                  >  1:2665118/1‑93 (MQ=255)
                                      cAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTaaa                                                  <  1:431679/93‑1 (MQ=255)
                                         aTTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATaa                                               >  1:1828903/1‑93 (MQ=255)
                                          ttCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAAc                                              <  1:1364320/93‑1 (MQ=255)
                                           tCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAAcc                                             <  1:2630266/93‑1 (MQ=255)
                                            cTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCa                                            <  1:2410533/93‑1 (MQ=255)
                                            cTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCa                                            <  1:3465948/93‑1 (MQ=255)
                                             tGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCAt                                           >  1:459582/1‑93 (MQ=255)
                                             tGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCAt                                           <  1:3410256/93‑1 (MQ=255)
                                              gCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATg                                          >  1:3956473/1‑93 (MQ=255)
                                              gCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATg                                          <  1:2993840/93‑1 (MQ=255)
                                              gCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATg                                          <  2:2893386/93‑1 (MQ=255)
                                               cATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGc                                         <  1:717814/93‑1 (MQ=255)
                                                aTTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCg                                        >  1:617399/1‑93 (MQ=255)
                                                aTTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCg                                        <  2:3944160/93‑1 (MQ=255)
                                                aTTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCg                                        >  2:4637341/1‑93 (MQ=255)
                                                 ttCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCgg                                       >  2:1548187/1‑93 (MQ=255)
                                                 ttCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCgg                                       >  2:2271349/1‑93 (MQ=255)
                                                    tcAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAga                                    >  2:3074607/1‑92 (MQ=255)
                                                    tcAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGt                                    >  1:1830926/1‑93 (MQ=255)
                                                    tcAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGt                                    <  2:503390/93‑1 (MQ=255)
                                                     cAGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTa                                   >  2:3312964/1‑93 (MQ=255)
                                                      aGACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTAc                                  >  2:2225179/1‑93 (MQ=255)
                                                       gACACAGCTTTCTACTCCATGGTAATGGGCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTAca                                 >  1:4547819/1‑92 (MQ=255)
                                                       gACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGACGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTAcc                                 >  2:2790085/1‑93 (MQ=255)
                                                       gACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTAcc                                 >  1:3370901/1‑93 (MQ=255)
                                                       gACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTAcc                                 >  2:3685328/1‑93 (MQ=255)
                                                       gACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTAcc                                 >  1:4030833/1‑93 (MQ=255)
                                                       gACACAGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTAca                                 >  2:3285936/1‑92 (MQ=255)
                                                         cacaGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCAc                               <  2:509474/93‑1 (MQ=255)
                                                          acaGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACt                              <  2:2699672/93‑1 (MQ=255)
                                                          acaGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACt                              >  1:4125696/1‑93 (MQ=255)
                                                           caGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTc                             >  2:2443343/1‑93 (MQ=255)
                                                           caGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTc                             >  1:3953594/1‑93 (MQ=255)
                                                           caGCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTc                             >  2:1869556/1‑93 (MQ=255)
                                                             gCTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCgc                           >  1:1122344/1‑92 (MQ=255)
                                                              cTTTCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGtc                          <  1:2185477/93‑1 (MQ=255)
                                                                ttCTACTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGtcat                        <  2:376223/93‑1 (MQ=255)
                                                                    aCTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACg                    >  1:3747829/1‑93 (MQ=255)
                                                                     cTCCATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGt                   <  2:4177568/93‑1 (MQ=255)
                                                                       ccATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTAt                 >  2:789056/1‑93 (MQ=255)
                                                                        cATGGTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTAtt                <  1:407722/93‑1 (MQ=255)
                                                                           ggTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTtat             <  1:1223970/93‑1 (MQ=255)
                                                                            gTAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTtata            <  2:2382089/93‑1 (MQ=255)
                                                                             tAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTtatat           <  2:1249786/93‑1 (MQ=255)
                                                                             tAATGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCCGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTtatat           <  2:885316/93‑1 (MQ=255)
                                                                              aaTGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTtatata          <  1:2457381/93‑1 (MQ=255)
                                                                              aaTGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTtatata          <  2:3225009/93‑1 (MQ=255)
                                                                              aaTGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTtatata          >  2:1369027/1‑93 (MQ=255)
                                                                                tGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATc        >  1:1333065/1‑93 (MQ=255)
                                                                                tGGCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATc        >  1:122724/1‑93 (MQ=255)
                                                                                 ggCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCa       <  1:4222061/93‑1 (MQ=255)
                                                                                  gCCCAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCAt      <  2:3127566/93‑1 (MQ=255)
                                                                                   cccAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATg     <  1:563838/93‑1 (MQ=255)
                                                                                     cAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGaa   >  2:2030701/1‑93 (MQ=255)
                                                                                     cAAACACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGaa   <  1:3815890/93‑1 (MQ=255)
                                                                                      aaaCACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAg  >  1:1403771/1‑93 (MQ=255)
                                                                                      aaaCACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAg  >  2:2478056/1‑93 (MQ=255)
                                                                                         |                                                                                         
CATTTTCAGAAGCTCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAG  >  NC_000913/2065777‑2065955

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