Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,997,977 T→C intergenic (‑64/‑390) b2865 ← / → b2866 putative lipoprotein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,997,9770TC98.5% 657.9 / ‑1.5 204intergenic (‑64/‑390)b2865/b2866putative lipoprotein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base C (124/77);  total (124/80)
Fisher's exact test for biased strand distribution p-value = 5.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 209 total aligned reads displayed.

TTTTTTATTCAGGCGTCCCGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCTCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTCGCTGTAGAGG  >  NC_000913/2997887‑2998066
                                                                                          |                                                                                         
ttttttATTCAGGCGTCCCGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATccct                                                                                         >  1:3675439/1‑93 (MQ=255)
ttttttATTCAGGCGTCCCGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATccct                                                                                         <  1:4032377/93‑1 (MQ=255)
ttttttATTCAGGCGTCCCGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATccct                                                                                         <  1:3373105/93‑1 (MQ=255)
ttttttATTCAGGCGTCCCGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATccct                                                                                         <  1:684732/93‑1 (MQ=255)
  ttttATTCAGGCGTCCCGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCAATCCCTTGc                                                                                      <  1:1349012/93‑1 (MQ=255)
       ttCAGGCGTCCCGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCt                                                                                  >  1:358021/1‑93 (MQ=255)
       ttCAGGCGTCCCGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCt                                                                                  >  2:3202160/1‑93 (MQ=255)
         cAGGCGTCCCGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTgg                                                                                <  1:3979761/93‑1 (MQ=255)
           ggCGTCCCGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGaa                                                                              >  2:1959409/1‑93 (MQ=255)
            gCGTCCCGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGaaa                                                                             >  1:909901/1‑93 (MQ=255)
             cGTCCCGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGaaaa                                                                            >  2:1448250/1‑93 (MQ=255)
              gTCCCGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGaaaaa                                                                           <  2:1707389/93‑1 (MQ=255)
              gTCCCGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGaaaaa                                                                           >  1:1526505/1‑93 (MQ=255)
               tCCCGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAAc                                                                          <  1:812945/93‑1 (MQ=255)
               tCCCGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAAc                                                                          >  2:357156/1‑93 (MQ=255)
                 ccGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCTCTTGCTTCTGGAAAAACGt                                                                        <  2:1569847/93‑1 (MQ=255)
                 ccGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGt                                                                        >  2:3256348/1‑93 (MQ=255)
                  cGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTg                                                                       >  2:3421398/1‑93 (MQ=255)
                  cGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTg                                                                       >  2:3888780/1‑93 (MQ=255)
                  cGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTg                                                                       <  1:169786/93‑1 (MQ=255)
                  cGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCATTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTg                                                                       <  1:4289265/93‑1 (MQ=255)
                   gCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTgg                                                                      <  2:1319301/93‑1 (MQ=255)
                    cACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTgga                                                                     >  1:716115/1‑92 (MQ=255)
                    cACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGt                                                                     <  1:1829448/93‑1 (MQ=255)
                    cACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGt                                                                     <  2:3115596/93‑1 (MQ=255)
                     aCTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTa                                                                    >  1:2391122/1‑93 (MQ=255)
                      cTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTaa                                                                   <  1:1597957/93‑1 (MQ=255)
                      cTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTaa                                                                   >  2:2405561/1‑93 (MQ=255)
                       tCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAAt                                                                  >  1:2422218/1‑93 (MQ=255)
                        cAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATa                                                                 <  1:3301000/93‑1 (MQ=255)
                        cAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATa                                                                 >  1:3362781/1‑93 (MQ=255)
                        cAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATa                                                                 <  1:2783353/93‑1 (MQ=255)
                         aataTCTTTACTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATag                                                                <  1:3342833/90‑1 (MQ=255)
                         aaaaTCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATag                                                                >  1:843185/1‑93 (MQ=255)
                          aaaTCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATaga                                                               <  2:1028319/93‑1 (MQ=255)
                          aaaTCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATaga                                                               >  2:421438/1‑93 (MQ=255)
                           aaTCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATagag                                                              >  1:2136298/1‑93 (MQ=255)
                           aaTCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATagag                                                              >  2:740334/1‑93 (MQ=255)
                           aaTCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATagag                                                              >  1:2284705/1‑93 (MQ=255)
                           aaTCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATagag                                                              >  1:1667133/1‑93 (MQ=255)
                            aTCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAgaga                                                             >  1:3936191/1‑93 (MQ=255)
                            aTCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAgaga                                                             >  1:157340/1‑93 (MQ=255)
                            aTCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAgaga                                                             <  2:3652456/93‑1 (MQ=255)
                            aTCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAgaga                                                             <  2:1178666/93‑1 (MQ=255)
                             tCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGaa                                                            >  1:4164084/1‑93 (MQ=255)
                             tCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGaa                                                            >  1:3852746/1‑93 (MQ=255)
                             tCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGaa                                                            >  2:461745/1‑93 (MQ=255)
                             tCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGaa                                                            <  1:3214193/93‑1 (MQ=255)
                               tttCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAag                                                          >  2:2176919/1‑93 (MQ=255)
                                ttCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAaga                                                         <  1:634706/93‑1 (MQ=255)
                                ttCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAaga                                                         >  2:3790261/1‑93 (MQ=255)
                                   cTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGt                                                      >  2:3307273/1‑93 (MQ=255)
                                       cTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAAt                                                  >  2:1050296/1‑93 (MQ=255)
                                        ttAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAAtt                                                 >  1:594831/1‑93 (MQ=255)
                                        ttAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAAtt                                                 >  2:1164968/1‑93 (MQ=255)
                                            aaaaCAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCAt                                             <  2:1883193/93‑1 (MQ=255)
                                             aaaCAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCAtt                                            >  2:2052878/1‑93 (MQ=255)
                                              aaCAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCAttt                                           >  1:4391304/1‑93 (MQ=255)
                                                 aaGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGc                                        >  2:4243482/1‑93 (MQ=255)
                                                  aGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCa                                       >  1:2854667/1‑93 (MQ=255)
                                                   aaCGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGGTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCaa                                      <  1:836034/92‑1 (MQ=255)
                                                        cTTATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGa                                 >  1:1095558/1‑93 (MQ=255)
                                                         ttATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGaa                                >  1:3607855/1‑93 (MQ=255)
                                                          tATCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAg                               >  2:3522278/1‑93 (MQ=255)
                                                            tCATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGCTAATTTCATTTGGCAAAGGGAAGTa                             >  1:4235367/1‑93 (MQ=255)
                                                             cATAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTat                            >  1:4186125/1‑93 (MQ=255)
                                                              aTAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTatc                           >  1:2158597/1‑93 (MQ=255)
                                                               tAGCAGCCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTatca                          <  2:673600/93‑1 (MQ=255)
                                                                agcagcCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTatcat                         <  1:1022375/93‑1 (MQ=255)
                                                                 gcagcCAGCCTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTatcatc                        >  1:2326729/1‑93 (MQ=255)
                                                                    gccagccTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGc                     >  1:4512807/1‑93 (MQ=255)
                                                                    gccagccTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGc                     >  2:2539303/1‑93 (MQ=255)
                                                                       agccTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCaaa                  >  2:2489091/1‑93 (MQ=255)
                                                                        gccTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCaaaa                 >  2:4146640/1‑93 (MQ=255)
                                                                         ccTCTTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAat                <  1:1751902/93‑1 (MQ=255)
                                                                          ctctTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAata               <  2:2401076/93‑1 (MQ=255)
                                                                          ctctTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAata               <  1:3512852/93‑1 (MQ=255)
                                                                           tctTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAatat              >  2:2842012/1‑93 (MQ=255)
                                                                            ctTGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATAtt             <  1:3407154/93‑1 (MQ=255)
                                                                             ttGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTc            >  2:87143/1‑93 (MQ=255)
                                                                             ttGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTc            >  2:1756433/1‑93 (MQ=255)
                                                                             ttGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTc            >  1:509754/1‑93 (MQ=255)
                                                                              tGATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTCg           >  1:2051342/1‑93 (MQ=255)
                                                                               gATACCCAATCCCTTGCTTCTGGAAAAACGTGGTCATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTCGc          <  2:3937803/93‑1 (MQ=255)
                                                                               gATACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTCGc          >  2:563313/1‑93 (MQ=255)
                                                                                 tACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTCGCtt        >  2:4032375/1‑92 (MQ=255)
                                                                                 tACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTCGCTg        >  1:878901/1‑93 (MQ=255)
                                                                                 tACCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTCGCTg        <  1:534077/93‑1 (MQ=255)
                                                                                  aCCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTCGCTGt       >  1:4485935/1‑93 (MQ=255)
                                                                                  aCCCAATCCCTTGCTTCTGGAAAAACGTGGTAATAAAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTCGCTGt       >  2:4676406/1‑93 (MQ=255)
                                                                                   cccAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTCGCTGTa      >  1:158419/1‑93 (MQ=255)
                                                                                    ccAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTCGCTGTag     >  2:2734491/1‑93 (MQ=255)
                                                                                     cAATCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTCGCTGTaga    <  1:4562644/93‑1 (MQ=255)
                                                                                      aaTCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTCGCTGTagag   >  1:263511/1‑93 (MQ=255)
                                                                                      aaTCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTCGCTGTagag   <  2:2603787/93‑1 (MQ=255)
                                                                                      aaTCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTCGCTGTagag   >  1:1485144/1‑93 (MQ=255)
                                                                                       aTCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAGAGGGAAGTATCATCAGCAAAATATTCGCTGTAGAgg  >  2:4750126/1‑93 (MQ=255)
                                                                                       aTCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTCGCTGTAGAgg  >  1:3010394/1‑93 (MQ=255)
                                                                                       aTCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTCGCTGTAGAgg  <  1:851815/93‑1 (MQ=255)
                                                                                       aTCCCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTCGCTGTAGAgg  >  2:3011098/1‑93 (MQ=255)
                                                                                          |                                                                                         
TTTTTTATTCAGGCGTCCCGCACTCAAAATCTTTCCTCGCTTAAAAAACAAGACGCCTTATCATAGCAGCCAGCCTCTTGATACCCAATCTCTTGCTTCTGGAAAAACGTGGTAATAGAGAAAGAGGTTAATTTCATTTGGCAAAGGGAAGTATCATCAGCAAAATATTCGCTGTAGAGG  >  NC_000913/2997887‑2998066

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