Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,791,507 T→C intergenic (‑182/‑107) htrL ← / → rfaD HtrL/ADP‑L‑glycero‑D‑mannoheptose‑6‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,791,5070TC99.7% 922.4 / ‑4.9 290intergenic (‑182/‑107)htrL/rfaDHtrL/ADP‑L‑glycero‑D‑mannoheptose‑6‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (245/44);  total (246/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.62e-01

Only 100 of 292 total aligned reads displayed.

AATATTGATACTAAAGCACTATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCATCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTCTCTGACTCCATA  >  NC_000913/3791417‑3791596
                                                                                          |                                                                                         
aaTATTGATACTAAAGCACTATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCAccc                                                                                         >  2:1012962/1‑93 (MQ=255)
 atatTGATACTAAAGCACTATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCt                                                                                        >  1:7072728/1‑93 (MQ=255)
  tatTGATACTAAAGCACTATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCtt                                                                                       >  1:4973894/1‑93 (MQ=255)
  tatTGATACTAAAGCACTATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCtt                                                                                       >  2:1914313/1‑93 (MQ=255)
    ttGATACTAAAGCACTATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGc                                                                                     >  1:4403638/1‑93 (MQ=255)
       aTACTAAAGCACTATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCAcc                                                                                  >  2:7581244/1‑93 (MQ=255)
         aCTAAAGCACTATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACctc                                                                                >  2:7167927/1‑93 (MQ=255)
          cTAAAGCACTATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACctct                                                                               >  2:5275905/1‑93 (MQ=255)
           tAAAGCACTATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTa                                                                              >  2:4733136/1‑93 (MQ=255)
           tAAAGCACTATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTa                                                                              >  1:7741870/1‑93 (MQ=255)
           tAAAGCACTATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTa                                                                              >  1:5298738/1‑93 (MQ=255)
           tAAAGCACTATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTa                                                                              >  2:7075657/1‑93 (MQ=255)
            aaaGCACTATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTAt                                                                             >  2:5857559/1‑93 (MQ=255)
                  cTATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAg                                                                       >  2:2931505/1‑93 (MQ=255)
                   tATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAgg                                                                      >  1:6885769/1‑93 (MQ=255)
                   tATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAgg                                                                      >  2:5950984/1‑93 (MQ=255)
                   tATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAgg                                                                      >  1:4076549/1‑93 (MQ=255)
                   tATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAgg                                                                      >  2:6600914/1‑93 (MQ=255)
                    aTTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAggg                                                                     >  2:1178266/1‑93 (MQ=255)
                    aTTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAggg                                                                     >  1:1351805/1‑93 (MQ=255)
                     ttCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGc                                                                    >  2:4899952/1‑93 (MQ=255)
                       cacaTGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTg                                                                  >  1:6909170/1‑93 (MQ=255)
                           tGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAgg                                                              >  2:4853760/1‑93 (MQ=255)
                           tGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAgg                                                              <  2:5764053/93‑1 (MQ=255)
                           tGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAgg                                                              >  2:7253635/1‑93 (MQ=255)
                              aaaaCCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGAt                                                           >  1:6302099/1‑93 (MQ=255)
                              aaaaCCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGAt                                                           >  2:3795495/1‑93 (MQ=255)
                              aaaaCCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGAt                                                           >  2:407441/1‑93 (MQ=255)
                               aaaCCAACATCCGCCATGGAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGAtt                                                          >  2:7742835/1‑93 (MQ=255)
                               aaaCCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGAtt                                                          <  1:1929006/93‑1 (MQ=255)
                                aaCCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTc                                                         >  1:2882966/1‑93 (MQ=255)
                                 aCCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCg                                                        >  2:1833824/1‑93 (MQ=255)
                                    aaCATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGAt                                                     >  2:4504449/1‑93 (MQ=255)
                                     aCATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGAtg                                                    >  2:419436/1‑93 (MQ=255)
                                     aCATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGAtg                                                    >  1:5430121/1‑93 (MQ=255)
                                      cATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGAtgt                                                   >  1:2812927/1‑93 (MQ=255)
                                       aTCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGAtgtg                                                  >  1:6663811/1‑93 (MQ=255)
                                       aTCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGAtgtg                                                  >  2:5986287/1‑93 (MQ=255)
                                       aTCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGAtgtg                                                  >  1:1544672/1‑93 (MQ=255)
                                        tCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGa                                                 >  1:5077859/1‑93 (MQ=255)
                                        tCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGa                                                 >  1:5924306/1‑93 (MQ=255)
                                        tCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGa                                                 >  2:206496/1‑93 (MQ=255)
                                          cGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATg                                               >  1:7424482/1‑93 (MQ=255)
                                          cGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATg                                               >  2:6449705/1‑93 (MQ=255)
                                            ccATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGt                                             >  1:4730237/1‑93 (MQ=255)
                                            ccATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGt                                             >  1:2144744/1‑93 (MQ=255)
                                              aTGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTAt                                           >  2:2311202/1‑93 (MQ=255)
                                               tGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATg                                          >  1:6092211/1‑93 (MQ=255)
                                                gAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATgt                                         >  1:1798339/1‑92 (MQ=255)
                                                gAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGa                                         >  1:2693565/1‑93 (MQ=255)
                                                  aGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGAtt                                       >  2:150178/1‑93 (MQ=255)
                                                  aGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGAtt                                       >  1:311580/1‑93 (MQ=255)
                                                   ggACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTa                                      <  1:576291/93‑1 (MQ=255)
                                                    gACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTaa                                     >  2:1381496/1‑92 (MQ=255)
                                                     aCTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACa                                    <  2:1356157/93‑1 (MQ=255)
                                                     aCTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACa                                    <  2:5632444/93‑1 (MQ=255)
                                                       tAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGa                                  >  1:2943626/1‑93 (MQ=255)
                                                       tAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGa                                  >  2:2585875/1‑93 (MQ=255)
                                                       tAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGa                                  >  2:3192284/1‑93 (MQ=255)
                                                          cTAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACAt                               >  1:926158/1‑93 (MQ=255)
                                                           tAAAACCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACAtt                              >  1:3814200/1‑93 (MQ=255)
                                                            aaaaCCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTc                             >  1:2329203/1‑93 (MQ=255)
                                                              aaCCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCgt                           >  2:3551051/1‑93 (MQ=255)
                                                               aCCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCgtg                          >  2:1303427/1‑93 (MQ=255)
                                                               aCCCAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCgtg                          >  2:554360/1‑93 (MQ=255)
                                                                 ccAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTc                        >  1:373817/1‑93 (MQ=255)
                                                                  cAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCt                       >  2:6233792/1‑93 (MQ=255)
                                                                  cAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCt                       >  2:6836481/1‑93 (MQ=255)
                                                                  cAAACTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCt                       >  2:4412253/1‑93 (MQ=255)
                                                                   aaaCTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTg                      >  2:3262814/1‑93 (MQ=255)
                                                                   aaaCTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTg                      >  2:2235552/1‑93 (MQ=255)
                                                                    aaCTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTga                     >  1:1428792/1‑93 (MQ=255)
                                                                     aCTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTgag                    >  2:3055181/1‑93 (MQ=255)
                                                                      cTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTgaga                   >  1:5046099/1‑93 (MQ=255)
                                                                      cTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTgaga                   <  1:4153584/93‑1 (MQ=255)
                                                                      cTAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTgaga                   <  1:6465753/93‑1 (MQ=255)
                                                                       tAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGAt                  >  1:6748981/1‑93 (MQ=255)
                                                                       tAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGAt                  >  2:7662547/1‑93 (MQ=255)
                                                                       tAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGAt                  <  2:6522610/93‑1 (MQ=255)
                                                                       tAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGAt                  >  2:2835949/1‑93 (MQ=255)
                                                                       tAGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGAt                  >  2:6911869/1‑93 (MQ=255)
                                                                        aGTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGAtt                 >  1:3464612/1‑93 (MQ=255)
                                                                         gTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTg                >  1:7581521/1‑93 (MQ=255)
                                                                         gTTTGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTg                >  2:7425001/1‑93 (MQ=255)
                                                                          tttGTTGCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGt               >  2:5153984/1‑93 (MQ=255)
                                                                           ttgttgCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGtc              >  2:1644061/1‑93 (MQ=255)
                                                                            tgttgCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGtct             >  1:3009435/1‑93 (MQ=255)
                                                                             gttgCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGtctc            >  2:2978766/1‑93 (MQ=255)
                                                                             gttgCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGtctc            <  1:7146010/93‑1 (MQ=255)
                                                                              ttgCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTctct           <  2:5289141/93‑1 (MQ=255)
                                                                              ttgCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTctct           >  2:2704028/1‑93 (MQ=255)
                                                                              ttgCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTctct           >  2:1792894/1‑93 (MQ=255)
                                                                                gCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTCTCTGa         <  1:2220648/93‑1 (MQ=255)
                                                                                gCAATTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTCTCTGa         >  1:3231919/1‑93 (MQ=255)
                                                                                  aaTTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTCTCTGACt       >  1:5695597/1‑93 (MQ=255)
                                                                                   aTTAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTCTCTGACTc      >  2:5672205/1‑93 (MQ=255)
                                                                                    ttAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTCTCTGACTcc     >  1:1230799/1‑93 (MQ=255)
                                                                                    ttAGCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTCTCTGACTcc     <  2:1203546/93‑1 (MQ=255)
                                                                                       gCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTCTCTGACTCCATa  <  1:1012962/93‑1 (MQ=255)
                                                                                       gCACCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTCTCTGACTCCATa  >  2:7868289/1‑93 (MQ=255)
                                                                                          |                                                                                         
AATATTGATACTAAAGCACTATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCATCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTCTCTGACTCCATA  >  NC_000913/3791417‑3791596

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