Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,964,041 A→G T4A (ACC→GCC)  rho → transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,964,0410AG98.1% 665.0 / 0.9 212T4A (ACC→GCC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base G (87/121);  total (88/124)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.88e-01

Only 100 of 218 total aligned reads displayed.

TAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGAC  >  NC_000913/3963951‑3964130
                                                                                          |                                                                                         
tAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGcc                                                                                         >  2:75526/1‑93 (MQ=255)
 aaaCAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCg                                                                                        <  2:126544/93‑1 (MQ=255)
   aCAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGaa                                                                                      >  2:2224032/1‑93 (MQ=255)
     aGGCATGGATGATCCTGCCATACCATTCACAACTTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAAtt                                                                                    <  1:4591917/93‑1 (MQ=255)
     aGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAAtt                                                                                    <  2:1464153/93‑1 (MQ=255)
      ggCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTa                                                                                   <  1:873599/93‑1 (MQ=255)
      ggCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTa                                                                                   >  2:3609902/1‑93 (MQ=255)
      atcATGGATGATCCTGCCATACCATTCACAACATGAAGTTCGAGATTTACCCCAAGTTTAATAACTCACACCACTATGAATCCTGCCGAATTa                                                                                   <  2:2472413/91‑1 (MQ=255)
       gCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTaa                                                                                  <  1:3423128/93‑1 (MQ=255)
         aTGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAg                                                                                >  2:3317064/1‑93 (MQ=255)
         aTGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTCTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAg                                                                                <  2:1982580/93‑1 (MQ=255)
         aTGGATGATCCTGCCATACCATTAACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAg                                                                                <  2:2623141/93‑1 (MQ=255)
          tGGATGATCCTACCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGa                                                                               >  2:2078889/1‑93 (MQ=255)
           ggATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGaa                                                                              <  2:3774263/93‑1 (MQ=255)
             atgatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATa                                                                            >  2:4708478/1‑93 (MQ=255)
              tgatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATAc                                                                           >  1:4372432/1‑93 (MQ=255)
              tgatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATAc                                                                           >  2:3557214/1‑93 (MQ=255)
              tgatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATAc                                                                           <  1:1207477/93‑1 (MQ=255)
              tgacCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATAc                                                                           <  2:3048860/93‑1 (MQ=255)
                 tCCTGCCATACCATTCACACCATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGcc                                                                        <  1:4202126/93‑1 (MQ=255)
                 tCCTGCCATACCATTCACAACATTAAGTTCGAGATGTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGcc                                                                        >  1:825976/1‑93 (MQ=255)
                   cTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCgg                                                                      <  1:1506150/93‑1 (MQ=255)
                   cTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCgg                                                                      <  1:1415005/93‑1 (MQ=255)
                   cTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCgg                                                                      >  2:2952005/1‑93 (MQ=255)
                    tGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGt                                                                     <  1:2377832/93‑1 (MQ=255)
                     gCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCggat                                                                    >  1:144773/1‑91 (MQ=255)
                       cATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTc                                                                  <  2:826089/93‑1 (MQ=255)
                        aTACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCt                                                                 <  2:509998/93‑1 (MQ=255)
                         tACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTg                                                                >  1:1232260/1‑93 (MQ=255)
                           ccATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAg                                                              >  1:2592195/1‑93 (MQ=255)
                            cATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGc                                                             <  1:3628916/93‑1 (MQ=255)
                            cATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGc                                                             >  1:2653524/1‑93 (MQ=255)
                            cATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGc                                                             >  2:258521/1‑93 (MQ=255)
                              ttCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATCCGCCGGTTTCTGAGCTg                                                           >  1:4084589/1‑93 (MQ=255)
                              ttCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTg                                                           >  2:2801582/1‑93 (MQ=255)
                              ttCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTg                                                           >  2:177172/1‑93 (MQ=255)
                               tCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGa                                                          <  1:4520170/93‑1 (MQ=255)
                               tCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGa                                                          <  2:4542359/93‑1 (MQ=255)
                                cacaACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGAt                                                         >  1:3674926/1‑93 (MQ=255)
                                cacaACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGAt                                                         <  1:1487518/93‑1 (MQ=255)
                                cacaACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGAt                                                         >  2:487938/1‑93 (MQ=255)
                                 acaacaTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATc                                                        >  2:4240866/1‑93 (MQ=255)
                                  caacaTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCa                                                       >  1:3590301/1‑93 (MQ=255)
                                  caacaTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCa                                                       >  1:956577/1‑93 (MQ=255)
                                   aacaTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCAc                                                      >  1:2725281/1‑93 (MQ=255)
                                   aacaTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCAc                                                      >  2:2522315/1‑93 (MQ=255)
                                   aacaTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCAc                                                      >  2:17477/1‑93 (MQ=255)
                                    acaTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCAct                                                     >  1:18126/1‑93 (MQ=255)
                                     caTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTGAAGAATACGCCGGTTTCTGAGCTGATAActc                                                    >  1:311146/1‑93 (MQ=255)
                                     aaTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCActc                                                    <  2:4738566/92‑1 (MQ=255)
                                      aTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCActct                                                   >  2:3630190/1‑93 (MQ=255)
                                       ttAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACtctc                                                  >  1:2208482/1‑93 (MQ=255)
                                        tAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCg                                                 <  2:2723236/93‑1 (MQ=255)
                                          aGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGc                                               <  2:3364610/93‑1 (MQ=255)
                                           gTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCg                                              >  1:79171/1‑93 (MQ=255)
                                           gTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCg                                              >  2:4144235/1‑93 (MQ=255)
                                             tCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGaa                                            <  2:4013695/93‑1 (MQ=255)
                                              cGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGaaa                                           <  1:3721380/93‑1 (MQ=255)
                                                agaTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAat                                         <  1:3775684/93‑1 (MQ=255)
                                                   tttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATg                                      <  1:4282125/93‑1 (MQ=255)
                                                   tttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATg                                      >  2:2897234/1‑93 (MQ=255)
                                                    ttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAATATATgg                                     >  1:1191674/1‑93 (MQ=255)
                                                     tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATggg                                    <  2:2915261/93‑1 (MQ=255)
                                                     tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATggg                                    >  1:1050905/1‑93 (MQ=255)
                                                     tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATggg                                    <  1:512487/93‑1 (MQ=255)
                                                     tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATggg                                    <  2:567695/93‑1 (MQ=255)
                                                     tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATggg                                    <  1:2141647/93‑1 (MQ=255)
                                                        cccAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCt                                 >  2:1253235/1‑93 (MQ=255)
                                                        cccAAAGTTAAGAACTCACACCACTAATAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCt                                 <  1:1688862/93‑1 (MQ=255)
                                                         ccAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTg                                >  1:551406/1‑93 (MQ=255)
                                                         ccAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTg                                <  2:3660711/93‑1 (MQ=255)
                                                         ccAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTg                                >  2:3210307/1‑93 (MQ=255)
                                                          cAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTgg                               <  2:43539/93‑1 (MQ=255)
                                                           aaGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGa                              >  1:3962921/1‑93 (MQ=255)
                                                             gTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGaaa                            <  1:3273719/93‑1 (MQ=255)
                                                               ttAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAAc                          <  2:66747/93‑1 (MQ=255)
                                                                  aGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTg                       <  2:4494956/93‑1 (MQ=255)
                                                                   gAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTgg                      <  2:3691102/93‑1 (MQ=255)
                                                                   gAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTgg                      <  1:3049344/93‑1 (MQ=255)
                                                                   aatcTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTgg                      <  2:1406865/90‑1 (MQ=255)
                                                                    aaCTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGc                     <  1:1900644/93‑1 (MQ=255)
                                                                    aaCTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGc                     <  2:728580/93‑1 (MQ=255)
                                                                    aaCTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGc                     <  2:4676024/93‑1 (MQ=255)
                                                                      cTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc                   >  1:1860496/1‑93 (MQ=255)
                                                                      cTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc                   >  2:1611927/1‑93 (MQ=255)
                                                                        cacaCCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGt                 <  2:1287625/93‑1 (MQ=255)
                                                                         acacCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGTGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTa                <  2:106800/93‑1 (MQ=255)
                                                                          caccacTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt               >  2:3059618/1‑93 (MQ=255)
                                                                          caccacTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt               <  1:4348487/93‑1 (MQ=255)
                                                                           accacTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATg              <  2:852338/93‑1 (MQ=255)
                                                                           accacTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATg              <  2:1862809/93‑1 (MQ=255)
                                                                            ccacTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGc             >  1:2599796/1‑93 (MQ=255)
                                                                              acTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGt           <  1:3405607/93‑1 (MQ=255)
                                                                              acTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGt           <  1:1437345/93‑1 (MQ=255)
                                                                                tATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa         <  1:1258345/93‑1 (MQ=255)
                                                                                 aTGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg        <  1:2843068/93‑1 (MQ=255)
                                                                                 aTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg        >  1:1577709/1‑93 (MQ=255)
                                                                                    aaTCTTGCCGAATTAAAGAATACGCCGGTTTCTGATCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAg     <  1:2891483/93‑1 (MQ=255)
                                                                                    aaTCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAg     <  2:3360244/93‑1 (MQ=255)
                                                                                      tCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGa   >  1:11378/1‑93 (MQ=255)
                                                                                          |                                                                                         
TAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGAC  >  NC_000913/3963951‑3964130

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: