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,0410AG99.8% 2000.7 / ‑5.2 615T4A (ACC→GCC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (290/323);  total (290/325)
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.63e-01

Only 100 of 625 total aligned reads displayed.

TAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGAC  >  NC_000913/3963951‑3964130
                                                                                          |                                                                                         
tAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGcc                                                                                         <  1:3350399/93‑1 (MQ=255)
tAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGcc                                                                                         >  1:2918288/1‑93 (MQ=255)
  aaCAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGa                                                                                       <  1:2288190/93‑1 (MQ=255)
  aaCAGGCAGGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGa                                                                                       <  1:8341809/93‑1 (MQ=255)
    aaGGCATGGATGTCCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAAt                                                                                     <  2:2748744/92‑1 (MQ=255)
     aGGCATGGATGATCCTGCCATAACATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACAACTATGAATCTTGCCGAAta                                                                                    >  2:2520868/1‑92 (MQ=255)
      ggCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTa                                                                                   <  2:314099/93‑1 (MQ=255)
         aTGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAg                                                                                >  2:7225216/1‑93 (MQ=255)
           ggATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGaa                                                                              >  1:6991287/1‑93 (MQ=255)
           ggATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGaa                                                                              <  2:5865833/93‑1 (MQ=255)
           ggATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGaa                                                                              <  2:3605121/93‑1 (MQ=255)
           gaatgatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGaa                                                                              <  1:9694003/91‑1 (MQ=255)
            gatgatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAAt                                                                             >  2:6838883/1‑93 (MQ=255)
            gatgatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAAt                                                                             <  2:4450557/93‑1 (MQ=255)
             atgatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATa                                                                            <  2:9542679/93‑1 (MQ=255)
             atgatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATa                                                                            <  2:5355279/93‑1 (MQ=255)
             atgatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATa                                                                            <  2:6775977/93‑1 (MQ=255)
              tgatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATAc                                                                           >  1:6849293/1‑93 (MQ=255)
                atCCTGCCATACCATTCACAACATTAAGTTCGAGGTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGc                                                                         <  2:7678774/93‑1 (MQ=255)
                atCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGc                                                                         <  1:396120/93‑1 (MQ=255)
                atCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGc                                                                         >  2:8604860/1‑93 (MQ=255)
                    tGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGt                                                                     <  1:3432456/93‑1 (MQ=255)
                    tGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGt                                                                     >  2:6984212/1‑93 (MQ=255)
                     gCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGtt                                                                    >  1:756485/1‑93 (MQ=255)
                       cATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTc                                                                  <  2:3321923/93‑1 (MQ=255)
                       cATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTc                                                                  <  2:7921263/93‑1 (MQ=255)
                        aTACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCt                                                                 <  1:8787412/93‑1 (MQ=255)
                        aTACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCt                                                                 <  1:5106546/93‑1 (MQ=255)
                         tACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTg                                                                >  2:5676602/1‑93 (MQ=255)
                         tACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTg                                                                >  1:7690138/1‑93 (MQ=255)
                         tACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTg                                                                >  2:4301318/1‑93 (MQ=255)
                         tACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTg                                                                >  1:818018/1‑93 (MQ=255)
                           ccATTCACAACATTAAGTTCGAGATTTACCCCAAGGTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAg                                                              >  1:6345120/1‑93 (MQ=255)
                            cATTCACAACATTAAATTCGATATTTACTCCAAGTTTAAGAACTCACACCACCATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGc                                                             <  2:337678/93‑1 (MQ=255)
                               tCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGa                                                          >  1:4628745/1‑93 (MQ=255)
                               tCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGa                                                          >  1:8673125/1‑93 (MQ=255)
                                cacaACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGAt                                                         >  2:7573436/1‑93 (MQ=255)
                                 acaacaTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATc                                                        <  1:10006156/93‑1 (MQ=255)
                                 acaacaTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATc                                                        >  2:1845586/1‑93 (MQ=255)
                                 acaacaTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATc                                                        <  1:9183698/93‑1 (MQ=255)
                                 acaacaTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATc                                                        <  1:7997476/93‑1 (MQ=255)
                                    acaTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCAct                                                     >  1:7268890/1‑93 (MQ=255)
                                    acaTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCAct                                                     >  2:1541812/1‑93 (MQ=255)
                                     caTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCActc                                                    <  1:1269297/93‑1 (MQ=255)
                                      aTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCActct                                                   >  2:4599067/1‑93 (MQ=255)
                                      aTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCActct                                                   >  2:2345099/1‑93 (MQ=255)
                                       ttAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACtctc                                                  >  1:3650061/1‑93 (MQ=255)
                                          aGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGc                                               <  1:2406558/93‑1 (MQ=255)
                                          aGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGc                                               <  1:3255191/93‑1 (MQ=255)
                                            ttCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGa                                             >  2:2212192/1‑93 (MQ=255)
                                              cGAGATTTACCCCAAGTAAAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGaaa                                           <  2:6291114/93‑1 (MQ=255)
                                               gagaTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGaaaa                                          >  1:3875400/1‑93 (MQ=255)
                                                 gaTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAata                                        <  1:7854600/93‑1 (MQ=255)
                                                 gaTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAata                                        <  1:3509634/93‑1 (MQ=255)
                                                 gaTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAata                                        <  1:667680/93‑1 (MQ=255)
                                                 gaTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTCAGCTGATCACTCTCGGCGAAActa                                        >  1:1444238/1‑93 (MQ=255)
                                                 aaTTTTCCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCCCTCTCGGCGAAAata                                        <  2:6499923/92‑1 (MQ=255)
                                                  aTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCTAAAatat                                       <  2:1698852/93‑1 (MQ=255)
                                                  aTTTACCCCAAGTTTAAGAACTCACACCACTACGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAatat                                       >  1:8971673/1‑93 (MQ=255)
                                                   tttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATg                                      <  2:6137219/93‑1 (MQ=255)
                                                   tttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATg                                      >  2:1438602/1‑93 (MQ=255)
                                                   tttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATg                                      >  1:5504414/1‑93 (MQ=255)
                                                   tttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATg                                      <  1:1556346/93‑1 (MQ=255)
                                                    ttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATgg                                     >  2:7433124/1‑93 (MQ=255)
                                                     tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATggg                                    <  2:6660011/93‑1 (MQ=255)
                                                      tccccAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATgggg                                   <  1:1832153/92‑1 (MQ=255)
                                                       ccccAAGTTTAAGACCTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGc                                  <  2:8812021/93‑1 (MQ=255)
                                                        cccAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCt                                 >  2:30232/1‑93 (MQ=255)
                                                        cccAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCt                                 >  2:9155742/1‑93 (MQ=255)
                                                         ccAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTg                                <  1:596432/93‑1 (MQ=255)
                                                         ccAAGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTg                                >  2:3190076/1‑93 (MQ=255)
                                                           aaGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGa                              <  1:5744826/93‑1 (MQ=255)
                                                            aGTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCGCTCTCGGCGAAAATATGGGGCTGGaa                             >  1:4906479/1‑93 (MQ=255)
                                                             gTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGgaa                            >  1:94704/1‑93 (MQ=255)
                                                             gTTTAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGaaa                            >  2:4922593/1‑93 (MQ=255)
                                                              tttAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGaaaa                           >  1:4248784/1‑93 (MQ=255)
                                                               ttAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACTCCGGGTTCTGAGCTGATACCTCTCGGCCAAAAGATGTGGCTGCAAAAc                          >  1:5319785/1‑93 (MQ=255)
                                                               ttAAGAACTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAAc                          <  2:1339845/93‑1 (MQ=255)
                                                                    aaCTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGc                     >  1:2706622/1‑93 (MQ=255)
                                                                      cTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAGTATGGGGCTGGAAAACCTGGCTc                   >  2:5501252/1‑93 (MQ=255)
                                                                      cTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc                   <  1:2140883/93‑1 (MQ=255)
                                                                      cTCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc                   <  2:8271898/93‑1 (MQ=255)
                                                                       tCACACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCg                  >  2:1131475/1‑93 (MQ=255)
                                                                        cacaCCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGt                 <  2:9814707/93‑1 (MQ=255)
                                                                         ccaccacTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTa                <  2:2031360/92‑1 (MQ=255)
                                                                         acacCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTa                <  2:8048035/93‑1 (MQ=255)
                                                                         acacCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTa                <  1:9115534/93‑1 (MQ=255)
                                                                          caccacTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt               <  1:1139098/93‑1 (MQ=255)
                                                                             cacTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCg            <  1:586507/93‑1 (MQ=255)
                                                                               cTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa          <  2:4098778/93‑1 (MQ=255)
                                                                               cTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa          <  2:3887550/93‑1 (MQ=255)
                                                                                tATGCATCTTGCAGACTTATAGAATACGCCGGTTTCTGAGCTGATCAATCTCGGCGAAAATATGGGGCTGGAAAAACTGGCTCATTTGCGTaa         <  2:8977891/93‑1 (MQ=255)
                                                                                tATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa         >  1:8475358/1‑93 (MQ=255)
                                                                                  tGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGc       >  2:2924091/1‑93 (MQ=255)
                                                                                  tGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGc       <  1:1706764/93‑1 (MQ=255)
                                                                                   gAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCa      >  1:1974121/1‑93 (MQ=255)
                                                                                     aTCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAgg    <  1:9912109/93‑1 (MQ=255)
                                                                                     aTCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAgg    >  2:9686057/1‑93 (MQ=255)
                                                                                      tCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGa   <  2:5048018/93‑1 (MQ=255)
                                                                                      tCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGa   <  2:9369534/93‑1 (MQ=255)
                                                                                          |                                                                                         
TAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGAC  >  NC_000913/3963951‑3964130

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