Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,962,227 T→C Y340C (TAC→TGC)  rhlB ← putative ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,962,2270TC99.6% 899.5 / ‑5.7 259Y340C (TAC→TGC) rhlBputative ATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (132/126);  total (133/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 267 total aligned reads displayed.

GCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACC  >  NC_000913/3962137‑3962316
                                                                                          |                                                                                         
 cAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAAt                                                                                        <  2:4112320/93‑1 (MQ=255)
 cAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAAt                                                                                        <  2:1324419/93‑1 (MQ=255)
  aGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATc                                                                                       <  1:2264802/93‑1 (MQ=255)
  aGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATc                                                                                       >  2:1243605/1‑93 (MQ=255)
   ggCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCt                                                                                      <  1:27200/93‑1 (MQ=255)
   ggCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCt                                                                                      <  1:954015/93‑1 (MQ=255)
   ggCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCt                                                                                      <  2:6220508/93‑1 (MQ=255)
     cAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTc                                                                                    <  1:3001701/93‑1 (MQ=255)
       aaTTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTcac                                                                                  <  2:2222285/93‑1 (MQ=255)
       aaTTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTcac                                                                                  >  1:5831086/1‑93 (MQ=255)
        aTTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTcaca                                                                                 <  2:2827630/93‑1 (MQ=255)
        aTTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTcaca                                                                                 >  2:6456811/1‑93 (MQ=255)
         ttCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAg                                                                                <  1:4925913/93‑1 (MQ=255)
           cAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTc                                                                              >  2:1939937/1‑93 (MQ=255)
            aaTGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCa                                                                             <  1:2558942/93‑1 (MQ=255)
              tGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGGTTGCGCCTGCGTGACCTGTACGACCAATACGGTGAACGCCATCTTCACAGTTATc                                                                           >  2:1879300/1‑93 (MQ=255)
              tGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATc                                                                           <  1:4559430/93‑1 (MQ=255)
               gCATACTCTTCACACGCCAGTCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCg                                                                          >  1:5219673/1‑93 (MQ=255)
               gCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCg                                                                          >  1:172075/1‑93 (MQ=255)
               gCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCg                                                                          >  2:540701/1‑93 (MQ=255)
                 aTACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCggg                                                                        >  2:4337708/1‑93 (MQ=255)
                  tACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGt                                                                       >  1:165422/1‑93 (MQ=255)
                   aCTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTa                                                                      >  1:2181201/1‑93 (MQ=255)
                   aCTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTa                                                                      >  1:3140766/1‑93 (MQ=255)
                    ctctTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTaa                                                                     >  2:5277266/1‑93 (MQ=255)
                    ctctTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTaa                                                                     >  1:2288828/1‑93 (MQ=255)
                       ttCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATc                                                                  <  2:5164161/93‑1 (MQ=255)
                         cacaCGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGt                                                                >  2:3072563/1‑93 (MQ=255)
                          acacGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTa                                                               >  1:3642925/1‑93 (MQ=255)
                          acacGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTa                                                               >  2:3288677/1‑93 (MQ=255)
                          acacGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTa                                                               >  2:384864/1‑93 (MQ=255)
                          acacGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTa                                                               <  1:2321483/93‑1 (MQ=255)
                          acacGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTa                                                               >  1:1098803/1‑93 (MQ=255)
                          acacGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTa                                                               >  2:4665068/1‑93 (MQ=255)
                          acacGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTa                                                               <  1:519423/93‑1 (MQ=255)
                           cacGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAg                                                              >  2:831235/1‑93 (MQ=255)
                           cacGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAg                                                              >  1:4848770/1‑93 (MQ=255)
                              gCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACTTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTa                                                           <  1:1334614/93‑1 (MQ=255)
                              gCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTa                                                           <  1:180054/93‑1 (MQ=255)
                                cAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTaaa                                                         >  1:3268664/1‑93 (MQ=255)
                                cAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTaaa                                                         <  2:1510487/93‑1 (MQ=255)
                                cAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTaaa                                                         <  1:2066178/93‑1 (MQ=255)
                                 aGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAg                                                        >  1:3528478/1‑93 (MQ=255)
                                 aGGCTGATAGAGTGACCGCGTGCGCCTGCGCGACCTGTACGACCAATACGGTGACCGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAg                                                        <  1:4746444/93‑1 (MQ=255)
                                  ggCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGa                                                       >  2:3028974/1‑93 (MQ=255)
                                  ggCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGa                                                       >  1:3315280/1‑93 (MQ=255)
                                   gCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGAc                                                      >  2:5997011/1‑93 (MQ=255)
                                   gCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGAc                                                      >  1:4444168/1‑93 (MQ=255)
                                    cTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACg                                                     <  1:1525621/93‑1 (MQ=255)
                                    cTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACg                                                     >  1:2939976/1‑93 (MQ=255)
                                     tGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGt                                                    <  2:475824/93‑1 (MQ=255)
                                     tGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGt                                                    <  1:5568805/93‑1 (MQ=255)
                                     tGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGt                                                    <  2:5352207/93‑1 (MQ=255)
                                      gATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTg                                                   >  2:6344996/1‑93 (MQ=255)
                                      gATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTg                                                   >  1:3404257/1‑93 (MQ=255)
                                       aTAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGc                                                  <  2:6110852/93‑1 (MQ=255)
                                         agagTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGt                                                <  1:5302311/93‑1 (MQ=255)
                                          gagTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTc                                               >  1:649770/1‑93 (MQ=255)
                                          gagTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTc                                               >  1:5459007/1‑93 (MQ=255)
                                          gagTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTc                                               <  1:5148757/93‑1 (MQ=255)
                                            gTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCAc                                             >  2:4257504/1‑93 (MQ=255)
                                              gACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGAGGTGCGTCACTg                                           >  1:4105026/1‑93 (MQ=255)
                                               aCCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGc                                          <  2:158748/93‑1 (MQ=255)
                                               aCCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGc                                          <  1:3506354/93‑1 (MQ=255)
                                                 cGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCg                                        <  1:84988/93‑1 (MQ=255)
                                                 cGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCg                                        >  1:3251495/1‑93 (MQ=255)
                                                   cTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGa                                      <  1:3868648/93‑1 (MQ=255)
                                                    ttGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGaa                                     <  1:4253851/93‑1 (MQ=255)
                                                     tGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGCCGTGCGTCACTGCCGGAat                                    >  2:1724742/1‑93 (MQ=255)
                                                     tGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAat                                    <  2:1200654/93‑1 (MQ=255)
                                                     tGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAat                                    >  1:2129158/1‑93 (MQ=255)
                                                      gcgcCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAata                                   <  1:6431996/93‑1 (MQ=255)
                                                      gcgcCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAata                                   <  2:1022856/93‑1 (MQ=255)
                                                       cgcCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAatat                                  <  2:1659853/93‑1 (MQ=255)
                                                        gcCTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATg                                 <  2:3455836/93‑1 (MQ=255)
                                                          cTGCGCGACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCa                               <  2:5621613/93‑1 (MQ=255)
                                                               cgACCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCa                          >  2:2757867/1‑93 (MQ=255)
                                                                 aCCTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCAcg                        >  1:6031776/1‑93 (MQ=255)
                                                                   cTGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCAcgcg                      >  1:1176948/1‑93 (MQ=255)
                                                                    tGTACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCAcgcgc                     >  2:164332/1‑93 (MQ=255)
                                                                      tACGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCgg                   <  2:50186/93‑1 (MQ=255)
                                                                       aCGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGc                  >  2:3371398/1‑93 (MQ=255)
                                                                       aCGACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGc                  <  2:5711707/93‑1 (MQ=255)
                                                                         gACCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCaa                >  1:2935423/1‑93 (MQ=255)
                                                                          aCCAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAAc               <  1:1789592/93‑1 (MQ=255)
                                                                           ccAATACGGTGAACTCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACg              >  2:6354351/1‑93 (MQ=255)
                                                                           ccAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACg              >  2:1683226/1‑93 (MQ=255)
                                                                           ccAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACg              <  2:5791494/93‑1 (MQ=255)
                                                                           ccAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACg              >  1:3188393/1‑93 (MQ=255)
                                                                            cAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGt             <  2:2506314/93‑1 (MQ=255)
                                                                            cAATACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGt             <  2:3606610/93‑1 (MQ=255)
                                                                             aaTACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTc            <  1:2487596/93‑1 (MQ=255)
                                                                              aTACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCg           >  2:5904236/1‑93 (MQ=255)
                                                                               tACGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCgg          >  2:3815977/1‑93 (MQ=255)
                                                                                aCGGTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGt         <  2:208521/93‑1 (MQ=255)
                                                                                  ggTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGACTATGCAAACCACGCGCGGCAACGTCGGTgg       >  2:2586116/1‑93 (MQ=255)
                                                                                  ggTGAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTgg       >  1:1923983/1‑93 (MQ=255)
                                                                                     gAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCaa    >  1:6008995/1‑93 (MQ=255)
                                                                                     gAACGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCaa    >  2:4489455/1‑93 (MQ=255)
                                                                                       aCGCAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAAcc  <  1:4884559/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCAGGCAAATTCAATGCATACTCTTCACACGCCAGGCTGATAGAGTGACCGCTTGCGCCTGCGCGACCTGTACGACCAATACGGTGAACGTAATCTTCACAGTCATCGGGTAAATCGTAGTTAAAGACGTGCGTCACTGCCGGAATATGCAAACCACGCGCGGCAACGTCGGTGGCAACC  >  NC_000913/3962137‑3962316

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