Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,099,661 (A)5→6 intergenic (‑19/+164) yggN ← / ← b2959 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,099,6561.A98.8% 939.2 / 2.5 249intergenic (‑14/+169)yggN/b2959hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base A (111/135);  total (114/135)
Fisher's exact test for biased strand distribution p-value = 9.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

Only 100 of 274 total aligned reads displayed.

TGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCT  >  NC_000913/3099564‑3099745
                                                                                             |                                                                                          
 gACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGa                                                                                           <  2:5641751/93‑1 (MQ=255)
   cGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaa                                                                                         <  2:4081313/93‑1 (MQ=255)
   cGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaa                                                                                         <  1:1658710/93‑1 (MQ=255)
    gCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaa                                                                                       <  2:1708483/93‑1 (MQ=255)
      tGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaa                                                                                     >  1:3705566/1‑92 (MQ=255)
       gCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaat                                                                                    >  2:2607533/1‑91 (MQ=255)
         aCTGGTAGGCGGCGTGAGCGGTCATGGCCGTCACTGAAAGTGCTGCCGCCAGCAGGATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaatca                                                                                  <  1:1754070/93‑5 (MQ=255)
           tGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAgg                                                                                 <  1:5552502/93‑1 (MQ=255)
           tGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAgg                                                                                 >  2:2140992/1‑93 (MQ=255)
           tGGTAGTCGGCGTGAGCGGTCATTGCCGGCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAgg                                                                                 <  2:1839275/93‑1 (MQ=255)
             gTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGcc                                                                               <  2:5589484/93‑1 (MQ=255)
              tAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCa                                                                              <  2:3327401/93‑1 (MQ=255)
              tAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCa                                                                              >  1:2671237/1‑93 (MQ=255)
                gTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGc                                                                            <  1:6142940/93‑1 (MQ=255)
                gTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGc                                                                            <  1:5228442/93‑1 (MQ=255)
                 tCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCa                                                                           <  1:2128720/93‑1 (MQ=255)
                 tCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCa                                                                           <  1:1285513/93‑1 (MQ=255)
                 tCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCa                                                                           >  2:2718144/1‑93 (MQ=255)
                   ggCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAtt                                                                         <  1:1698257/93‑1 (MQ=255)
                   ggCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAtt                                                                         <  1:4616832/93‑1 (MQ=255)
                     cGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTg                                                                       <  2:2459986/93‑1 (MQ=255)
                      gTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGc                                                                      <  1:2793589/93‑1 (MQ=255)
                      gTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGc                                                                      >  1:5379296/1‑93 (MQ=255)
                       tGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCt                                                                     >  2:4606573/1‑93 (MQ=255)
                       tGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCt                                                                     <  1:4136110/93‑1 (MQ=255)
                       tGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCt                                                                     >  2:6116039/1‑93 (MQ=255)
                       tGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCt                                                                     <  2:5233758/93‑1 (MQ=255)
                        gAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTg                                                                    <  2:3634754/93‑1 (MQ=255)
                        gAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTg                                                                    <  1:4408059/93‑1 (MQ=255)
                           cGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGcc                                                                 <  2:1134260/93‑1 (MQ=255)
                           cGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAGAATCAGGCCAGCATTTGCTGGcc                                                                 >  1:3486109/1‑93 (MQ=255)
                             gTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGcccc                                                               >  2:2380154/1‑93 (MQ=255)
                             gTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGcccc                                                               >  1:3382676/1‑93 (MQ=255)
                             gTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGcccc                                                               <  2:3704669/93‑1 (MQ=255)
                              tCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCa                                                              <  1:986008/93‑1 (MQ=255)
                              tCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCa                                                              >  1:2891928/1‑93 (MQ=255)
                                aTTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGa                                                            >  2:3239937/1‑93 (MQ=255)
                                aTTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGa                                                            >  2:3942283/1‑93 (MQ=255)
                                 ttGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGAt                                                           >  2:3422999/1‑93 (MQ=255)
                                 ttGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGAt                                                           <  1:3203011/93‑1 (MQ=255)
                                 ttGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGAt                                                           >  2:3574225/1‑93 (MQ=255)
                                  tGCCGTCACTTAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGAtt                                                          <  2:2773768/93‑1 (MQ=255)
                                      gTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATa                                                      <  1:2511716/93‑1 (MQ=255)
                                       tCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATaa                                                     >  1:1198023/1‑93 (MQ=255)
                                       tCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATaa                                                     >  2:5828154/1‑93 (MQ=255)
                                       tCACTGAAAGTGCTGCCCCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATaa                                                     <  2:2790274/93‑1 (MQ=255)
                                            gAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg                                                <  1:5313797/93‑1 (MQ=255)
                                            gAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg                                                >  2:5531598/1‑93 (MQ=255)
                                             aaaGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGt                                               <  1:4198610/93‑1 (MQ=255)
                                              aaGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTa                                              <  2:4374757/93‑1 (MQ=255)
                                              aaGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTa                                              <  2:2420812/93‑1 (MQ=255)
                                               aGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAc                                             <  2:3076145/93‑1 (MQ=255)
                                               aGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAc                                             <  1:477221/93‑1 (MQ=255)
                                               aGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAc                                             <  1:3057847/93‑1 (MQ=255)
                                                gTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcg                                            <  2:603236/93‑1 (MQ=255)
                                                  gctgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTAAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcgcg                                          <  2:4456333/93‑1 (MQ=255)
                                                   ctgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGt                                         >  2:148593/1‑93 (MQ=255)
                                                     gccgccAGCAGCATTTTGCGCATCATAGTCTTCGCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGttt                                       >  1:567016/1‑93 (MQ=255)
                                                     gccgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGttt                                       >  2:1963540/1‑93 (MQ=255)
                                                       cgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTc                                     <  1:721115/93‑1 (MQ=255)
                                                        gccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTcc                                    >  1:145197/1‑93 (MQ=255)
                                                        gccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTcc                                    <  2:4290450/93‑1 (MQ=255)
                                                         ccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCa                                   >  2:5016167/1‑93 (MQ=255)
                                                         ccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCa                                   <  1:4705860/93‑1 (MQ=255)
                                                           agcagcATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATg                                 >  1:3004756/1‑93 (MQ=255)
                                                           agcagcATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATg                                 >  1:1109539/1‑93 (MQ=255)
                                                              agcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCg                              >  2:826224/1‑93 (MQ=255)
                                                               gcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCgg                             >  2:5722843/1‑93 (MQ=255)
                                                                caTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGa                            >  2:2535285/1‑93 (MQ=255)
                                                                caTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGa                            <  1:4884571/93‑1 (MQ=255)
                                                                 aTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGAt                           >  2:903819/1‑93 (MQ=255)
                                                                  ttttGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATg                          >  2:4500677/1‑93 (MQ=255)
                                                                  ttttGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATg                          >  2:1546906/1‑93 (MQ=255)
                                                                  ttttGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATg                          >  2:4753488/1‑93 (MQ=255)
                                                                       cgcATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcggcg                     >  1:227002/1‑93 (MQ=255)
                                                                        gcATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcggcgc                    >  1:5027301/1‑92 (MQ=255)
                                                                        gcATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcggcgc                    >  1:1519828/1‑92 (MQ=255)
                                                                         catcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcggcgca                   >  2:3187421/1‑91 (MQ=255)
                                                                         catcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTa                   <  2:92917/93‑1 (MQ=255)
                                                                          atcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGCaa                  >  1:3310038/1‑93 (MQ=255)
                                                                           tcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTaaa                 >  2:2253047/1‑93 (MQ=255)
                                                                           tcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGCaaa                 <  2:5368994/93‑1 (MQ=255)
                                                                           tcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCAAtctgg                              >  2:1218752/1‑76 (MQ=255)
                                                                            catAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGCAAac                >  1:3865508/1‑93 (MQ=255)
                                                                            catAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGCAAac                <  2:4186517/93‑1 (MQ=255)
                                                                             atAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTaaaac               >  1:1907267/1‑91 (MQ=255)
                                                                             atAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAaca               >  1:2856163/1‑93 (MQ=255)
                                                                             atAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGCAAaca               >  1:4509366/1‑93 (MQ=255)
                                                                              tAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAacac              >  1:2728444/1‑93 (MQ=255)
                                                                                gTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCt            >  1:214860/1‑93 (MQ=255)
                                                                                gTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCt            >  2:5132281/1‑93 (MQ=255)
                                                                                 tCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCtt           <  1:3587699/93‑1 (MQ=255)
                                                                                  cTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTa          >  2:5771561/1‑93 (MQ=255)
                                                                                  cTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGCAAACACCTTa          <  1:1017048/93‑1 (MQ=38)
                                                                                   ttCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGCAAACACCTTAt         <  1:5844033/93‑1 (MQ=255)
                                                                                    tCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGCAAACACCTTATc        >  1:2088979/1‑93 (MQ=255)
                                                                                    tCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGCAAACACCTTATc        >  2:2072259/1‑93 (MQ=255)
                                                                                     cccTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATcc       >  1:3138277/1‑93 (MQ=255)
                                                                                     cccTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATcc       >  1:5149078/1‑93 (MQ=255)
                                                                                     cccTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATcc       >  2:735679/1‑93 (MQ=255)
                                                                                             |                                                                                          
TGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCT  >  NC_000913/3099564‑3099745

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