Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 259,241 A→G V287A (GTC→GCC)  phoE ← outer membrane porin PhoE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913259,2410AG99.5% 1673.0 / NA 405V287A (GTC→GCC) phoEouter membrane porin PhoE
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (222/180);  total (225/180)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 405 total aligned reads displayed.

CTGACATATTTTTGTTGAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGT  >  NC_000913/259143‑259338
                                                                                                  |                                                                                                 
cTGACATATTTTTGTTGAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCa                                                                                                 >  1:506867/1‑101 (MQ=255)
cTGACATATTTTTGTTGAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCa                                                                                                 <  2:506276/101‑1 (MQ=255)
 tGACATATTTTTGTTGAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCAt                                                                                                <  2:505924/101‑1 (MQ=255)
  gACATATTTTTGTTGAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATa                                                                                               >  2:506871/1‑101 (MQ=255)
   aCATATTTTTGTTGAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATaa                                                                                              >  1:506876/1‑101 (MQ=255)
     atatTTTTGTTGAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAAcc                                                                                            >  2:506883/1‑101 (MQ=255)
     atatTTTTGTTGAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAAcc                                                                                            >  1:506880/1‑101 (MQ=255)
      tatTTTTGTTGAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAAccc                                                                                           >  2:506890/1‑101 (MQ=255)
       atTTTTGTTGAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCa                                                                                          >  1:506891/1‑101 (MQ=255)
        tttttGTTGAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAg                                                                                         >  2:506896/1‑101 (MQ=255)
         ttttGTTGAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGc                                                                                        >  1:506902/1‑101 (MQ=255)
         ttttGTTGAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGc                                                                                        <  1:506438/101‑1 (MQ=255)
             gttgAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATg                                                                                    >  1:506909/1‑101 (MQ=255)
              ttgAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATgg                                                                                   <  2:506354/101‑1 (MQ=255)
              ttgAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATgg                                                                                   <  1:506236/101‑1 (MQ=255)
              ttgAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATgg                                                                                   <  2:506187/101‑1 (MQ=255)
               tgAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGa                                                                                  <  2:506081/101‑1 (MQ=255)
                gAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGAc                                                                                 >  2:506912/1‑101 (MQ=255)
                gAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGAc                                                                                 >  1:506916/1‑101 (MQ=255)
                 aaaTAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGc                                                                                        <  2:506154/93‑1 (MQ=255)
                  aataatACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGc                                                                               >  1:506923/1‑101 (MQ=255)
                  aataatACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGc                                                                               >  2:506921/1‑101 (MQ=255)
                   ataataCGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCa                                                                              >  1:506928/1‑101 (MQ=255)
                   ataataCGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCa                                                                              >  2:506930/1‑101 (MQ=255)
                    taataCGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAg                                                                             <  1:506594/101‑1 (MQ=255)
                    taataCGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAg                                                                             <  1:506305/101‑1 (MQ=255)
                       taCGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGAcc                                                                          >  2:506936/1‑101 (MQ=255)
                       taCGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGAcc                                                                          >  1:506939/1‑101 (MQ=255)
                          gTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCaaa                                                                       >  2:506946/1‑101 (MQ=255)
                           tAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAg                                                                      >  1:506948/1‑101 (MQ=255)
                           tAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAg                                                                      <  1:506383/101‑1 (MQ=255)
                             gCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTc                                                                    <  2:506312/101‑1 (MQ=255)
                             gCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTc                                                                    <  1:506215/101‑1 (MQ=255)
                              cACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCAAGCGATGGACGCAGACCAAAGTCa                                                                   >  2:506953/1‑101 (MQ=255)
                                ccGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCaaa                                                                 <  1:506660/101‑1 (MQ=255)
                                ccGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCaaa                                                                 >  2:506961/1‑101 (MQ=255)
                                 cGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAAc                                                                <  1:506243/101‑1 (MQ=255)
                                 cGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAAc                                                                <  2:506489/101‑1 (MQ=255)
                                  gACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACt                                                               <  2:506408/101‑1 (MQ=255)
                                  gACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACt                                                               >  1:506965/1‑101 (MQ=255)
                                  gACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACt                                                               <  1:506419/101‑1 (MQ=255)
                                     gTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGt                                                            <  1:506074/101‑1 (MQ=255)
                                      tCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTa                                                           <  2:506632/101‑1 (MQ=255)
                                      tCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTa                                                           >  2:506971/1‑101 (MQ=255)
                                       cGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTAt                                                          >  1:506974/1‑101 (MQ=255)
                                        gATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTAtt                                                         >  1:506978/1‑101 (MQ=255)
                                         atatAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTg                                                        >  2:506983/1‑101 (MQ=255)
                                          tataATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGa                                                       <  2:506472/101‑1 (MQ=255)
                                          tataATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGa                                                       >  1:506986/1‑101 (MQ=255)
                                           ataATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAg                                                      >  2:506990/1‑101 (MQ=255)
                                            taATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGc                                                     <  1:505996/101‑1 (MQ=255)
                                            taATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGc                                                     <  2:506251/101‑1 (MQ=255)
                                             agttGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGc                                                     >  1:506994/3‑100 (MQ=255)
                                              aTTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGa                                                   <  1:505843/101‑1 (MQ=255)
                                                tGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGAcc                                                 <  2:506335/101‑1 (MQ=255)
                                                tGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGAcc                                                 >  2:506999/1‑101 (MQ=255)
                                                  aCCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGc                                               >  1:507002/1‑101 (MQ=255)
                                                     aGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTc                                            <  2:506223/101‑1 (MQ=255)
                                                     aGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTc                                            >  2:507005/1‑101 (MQ=255)
                                                     aGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTc                                            <  2:506378/101‑1 (MQ=255)
                                                      gATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCa                                           >  1:507010/1‑101 (MQ=255)
                                                       aTCTTCATCACCGATCCCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCaa                                                                  <  1:506287/77‑1 (MQ=255)
                                                       aTCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCaa                                          >  2:507014/1‑101 (MQ=255)
                                                       aTCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCaa                                          <  2:506545/101‑1 (MQ=255)
                                                         cTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAg                                        <  2:506507/101‑1 (MQ=255)
                                                          ttCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGt                                       >  1:507021/1‑101 (MQ=255)
                                                          ttCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGt                                       >  2:507020/1‑101 (MQ=255)
                                                           tcatcaCCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGtt                                      <  1:506500/101‑1 (MQ=255)
                                                            catcaCCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAg                                        <  2:506258/98‑1 (MQ=255)
                                                            catcaCCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTc                                     >  1:507028/1‑101 (MQ=255)
                                                            catcaCCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTc                                     >  2:507029/1‑101 (MQ=255)
                                                              tcaCCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCtg                                   >  2:507034/1‑101 (MQ=255)
                                                               caCCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCtgt                                  >  1:507038/1‑101 (MQ=255)
                                                                aCCGATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCtgtg                                 <  1:506532/101‑1 (MQ=255)
                                                                   gATACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCt                              >  1:507043/1‑101 (MQ=255)
                                                                    aTACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCtt                             >  2:507045/1‑101 (MQ=255)
                                                                     tACCTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTa                            <  1:506367/101‑1 (MQ=255)
                                                                       ccTTCAATATCTTTCCCTTTCGATAAGGCATACCCCAGCGATGGACGCAGGCCAAAGt                                                                     <  2:506201/58‑1 (MQ=255)
                                                                       ccTTCAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTAtt                          >  2:507051/1‑101 (MQ=255)
                                                                           cAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTggta                      >  1:507056/1‑99 (MQ=255)
                                                                           cAATATCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCa                      >  2:507058/1‑101 (MQ=255)
                                                                              tatCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAg                   >  1:507062/1‑101 (MQ=255)
                                                                               atCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAgc                  <  2:506903/101‑1 (MQ=255)
                                                                               atCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAgc                  <  1:506508/101‑1 (MQ=255)
                                                                                tCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAgcc                 >  2:507067/1‑101 (MQ=255)
                                                                                tCTTTCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAgcc                 <  1:506462/101‑1 (MQ=255)
                                                                                   ttCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAgccgcc              <  1:506274/101‑1 (MQ=255)
                                                                                   ttCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAgccgcc              <  1:506349/101‑1 (MQ=255)
                                                                                   ttCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAgccgcc              <  2:506434/101‑1 (MQ=255)
                                                                                    tCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAgccgc               >  1:507072/1‑99 (MQ=255)
                                                                                    tCCCTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCa             <  1:506404/101‑1 (MQ=255)
                                                                                     cccTTTCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAg            >  2:507073/1‑101 (MQ=255)
                                                                                        tttCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTa         >  2:507080/1‑101 (MQ=255)
                                                                                         ttCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTAt        >  1:507081/1‑101 (MQ=255)
                                                                                         ttCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTAt        <  1:506156/101‑1 (MQ=255)
                                                                                          tCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTAtt       <  2:506583/101‑1 (MQ=255)
                                                                                          tCGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTAtt       <  1:506721/101‑1 (MQ=255)
                                                                                           cGATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTg      >  1:507088/1‑101 (MQ=255)
                                                                                           cGATAAGACATAACCCAGCGATAGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTg      >  2:507090/1‑101 (MQ=255)
                                                                                            gATAAGGCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTa         >  1:507093/1‑97 (MQ=255)
                                                                                                  |                                                                                                 
CTGACATATTTTTGTTGAAATAATACGTAGCACCGACGTCGATATAATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGT  >  NC_000913/259143‑259338

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