Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,111,377 A→G 100% D172G (GAT→GGT)  opgH → OPG biosynthetic ACP‑dependent transmembrane UDP‑glucose beta‑1,2 glycosyltransferase; nutrient‑dependent cell size regulator, FtsZ assembly antagonist

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,111,3770AG98.2% 384.8 / ‑3.0 121D172G (GAT→GGT) opgHOPG biosynthetic ACP‑dependent transmembrane UDP‑glucose beta‑1,2 glycosyltransferase; nutrient‑dependent cell size regulator, FtsZ assembly antagonist
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (57/60);  total (61/60)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

Only 100 of 121 total aligned reads displayed.

TCTGTTGATCCTGACGCTCGCGCAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGATATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGATCCTCTTTGCGGTACTGTTCTGTTG  >  NC_000913/1111282‑1111473
                                                                                               |                                                                                                
tCTGTTGATCCTGACGCTCGCGCAAACTGTCGTCGCGATCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTAtg                                                                                              <  1:394470/100‑1 (MQ=255)
tCTGTTGATCCTGACGCTCGCGCAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTAtg                                                                                              <  2:868813/100‑1 (MQ=255)
tCTGTTGATCCTGACGCTCGCGCAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTAtg                                                                                              >  1:2800058/1‑100 (MQ=255)
  tGTTGATCCTGACGCTCGCGCAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTAtggt                                                                                            <  1:2656078/100‑1 (MQ=255)
   gTTGATCCTGACGCTCGCGCAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGAGCGCTGATTAATCCTATGGGTAtggtt                                                                                           >  1:613585/1‑100 (MQ=255)
     tGATCCTGACGCTCGCGCAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTAtggttgg                                                                                         <  1:2500721/100‑1 (MQ=255)
     tGATCCTGACGCTCGCGCAAACTGACGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTAtggttgg                                                                                         >  2:405458/1‑100 (MQ=255)
      gATCCTGACGCTCGCGCAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTAtggttggt                                                                                        <  1:362330/100‑1 (MQ=255)
      gATCCTGACGCTCGCGCAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTAtggttggt                                                                                        >  1:1306315/1‑100 (MQ=255)
      gATCCTGACGCTCGCGCAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTAtggttggt                                                                                        >  2:2791251/1‑100 (MQ=255)
      gATCCTGACGCTCGCGCAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTAtggttggt                                                                                        <  2:167921/100‑1 (MQ=255)
         ccTGACGCTCGCGCAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAg                                                                                     <  1:2581944/100‑1 (MQ=255)
           tGACGCTCGCGCAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGa                                                                                   >  1:902581/1‑100 (MQ=255)
               gCTCGCGCAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTtg                                                                               <  2:1286383/100‑1 (MQ=255)
                 tCGCGCAAACTGTCGTCGCGACCTGGTATGTGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTtgtg                                                                             <  1:2261693/100‑1 (MQ=255)
                 tCGCGCAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTtgtg                                                                             <  2:901880/100‑1 (MQ=255)
                  cgcgCAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTgg                                                                            >  2:2637771/1‑100 (MQ=255)
                  cgcgCAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTgg                                                                            >  2:2955209/1‑100 (MQ=255)
                  cgcgCAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTgg                                                                            >  1:2481192/1‑100 (MQ=255)
                  cgcgCAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTgg                                                                            >  2:636284/1‑100 (MQ=255)
                   gcgcAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTggg                                                                           >  1:979779/1‑100 (MQ=255)
                   gcgcAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTggg                                                                           <  1:1149666/100‑1 (MQ=255)
                   gcgcAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTgg                                                                            <  2:1097491/99‑1 (MQ=255)
                        aaCTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTcc                                                                      >  1:575083/1‑100 (MQ=255)
                        aaCTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGATATGGTTGGTCAGGATTTGTGGGTTTcc                                                                      >  2:921344/1‑100 (MQ=255)
                         aCTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCt                                                                     <  2:1406015/100‑1 (MQ=255)
                          cTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCtt                                                                    <  2:1085891/100‑1 (MQ=255)
                          cTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCtt                                                                    <  2:2874478/100‑1 (MQ=255)
                            gtcgtcGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTa                                                                  >  2:2186071/1‑100 (MQ=255)
                             tcgtcgCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTAt                                                                 <  2:2172752/100‑1 (MQ=255)
                               gtcgCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGc                                                               <  2:3003623/100‑1 (MQ=255)
                               gtcgCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGAGGGTTTCCTTTATGc                                                               >  1:3003656/1‑100 (MQ=255)
                                   cgACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCt                                                           >  1:454676/1‑100 (MQ=255)
                                   cgACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCt                                                           >  1:222349/1‑100 (MQ=255)
                                    gACCTGGTATATGAAGACCATTCTTTCTTATCAGTGTTGGGCGTTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCtt                                                          <  1:1716656/100‑1 (MQ=255)
                                     aCCTTGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGTTGATTAATCCTATGGGTATGGTTGGTCAGTATTTGTGGGTTTCCTTTATGCAGCTTc                                                         <  2:1828624/100‑1 (MQ=255)
                                     aCCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTc                                                         >  1:2174248/1‑100 (MQ=255)
                                     aCCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTc                                                         >  1:1028670/1‑100 (MQ=255)
                                       cTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTg                                                       <  1:2709318/100‑1 (MQ=255)
                                        tGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGc                                                      >  2:1928065/1‑100 (MQ=255)
                                        tGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGc                                                      >  1:1326288/1‑100 (MQ=255)
                                            atatGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTta                                                  <  1:1958469/100‑1 (MQ=255)
                                             tatGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCCTCTGCCTtat                                                 >  2:2463719/1‑100 (MQ=255)
                                             tatGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTAGTCAGGATTTGTGGGTTTCCTTTATGAAGCTTCTGCCTtat                                                 >  1:2578186/1‑100 (MQ=255)
                                               tGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTatat                                               >  1:422185/1‑100 (MQ=255)
                                               tGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTatat                                               >  2:1449606/1‑100 (MQ=255)
                                                  aGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATAtgct                                            <  2:1193585/100‑1 (MQ=255)
                                                   gACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATAtgctg                                           <  1:395443/100‑1 (MQ=255)
                                                   gACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATAtgctg                                           <  1:1066256/100‑1 (MQ=255)
                                                   gACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGCTTTCCTTTATGCAGCTTCTGCCTTATAtgctg                                           <  1:2887350/100‑1 (MQ=255)
                                                    aCCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATGCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATAtgctgc                                          <  1:2130406/100‑1 (MQ=255)
                                                    aCCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATAtgctgc                                          <  1:1191822/100‑1 (MQ=255)
                                                     ccATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGTTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCa                                         <  2:2268972/100‑1 (MQ=255)
                                                     ccATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATAtgctgc                                          >  2:7019/1‑99 (MQ=255)
                                                      naTTCTTCCTTATCAGGGTTGGACGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTACTTTATGCAGCTTCTGCCTTATATGCTGCa                                         >  2:645363/2‑99 (MQ=255)
                                                      cATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCaa                                        >  2:2508911/1‑100 (MQ=255)
                                                      cATGCTTCCTTATCAGGGTTGGGCGTTGATTAATACTATAGGTATGGTTGATGAGGATTGGTGGGTTTGCTTTCTACAGCTTGTGCTTTATAtgct                                            >  2:1809284/1‑96 (MQ=255)
                                                       aTTCTTCTTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATAtgctgc                                          <  1:536568/97‑1 (MQ=255)
                                                       aTTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCaaa                                       >  2:2773857/1‑100 (MQ=255)
                                                       aTTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCaaa                                       >  1:2248001/1‑100 (MQ=255)
                                                         tcttcCTTATCAGGGTTGGGCGCTGATTAATTCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAAcc                                     <  1:1817369/100‑1 (MQ=255)
                                                         tcttcCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAAcc                                     >  2:2714692/1‑100 (MQ=255)
                                                         tcttcCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAAcc                                     <  1:2752063/100‑1 (MQ=255)
                                                          cttcCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCg                                    >  2:2788800/1‑100 (MQ=255)
                                                           ttgcTTATCAGGGTTGGGCGGTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCgg                                   <  2:1576026/97‑1 (MQ=255)
                                                           ttcCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGTGTTTCCTGTATGCAGCTTCTGCCTTATATGCTGCAAACCgg                                   <  2:616355/100‑1 (MQ=255)
                                                           ttcCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCgg                                   <  2:1789070/100‑1 (MQ=255)
                                                           ttcCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCgg                                   <  2:2526385/100‑1 (MQ=255)
                                                            tcCTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGATTCCTTTATGCAAATTCTACCTTATATGCTGCAAACCGGt                                  >  1:1043859/1‑100 (MQ=255)
                                                            tcCTTATCAGGGTTGGGCGCTGATGAATTCTATGGGTATGGTTGGTCAGGATTTGTGGTTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCgg                                   <  2:612349/99‑1 (MQ=255)
                                                             ccTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTa                                 <  1:2132776/100‑1 (MQ=255)
                                                              cTTATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTAt                                >  2:2658558/1‑100 (MQ=255)
                                                               ttATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATc                               >  2:1247646/1‑100 (MQ=255)
                                                                tATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATcc                              >  1:9738/1‑100 (MQ=255)
                                                                tATCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATcc                              >  1:2497037/1‑100 (MQ=255)
                                                                  tCAGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTAATATGCTGCAAACCGGTATCCTg                            >  1:981404/1‑100 (MQ=255)
                                                                    aGGGTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGAt                          <  2:2817559/100‑1 (MQ=255)
                                                                      ggTTGGGCGGTGAGTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCGTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGATcc                        <  2:2641071/100‑1 (MQ=255)
                                                                      ggTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGATcc                        <  2:1037166/100‑1 (MQ=255)
                                                                      ggTTGGGCGCTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGCATCCTGATcc                        >  2:2640556/1‑100 (MQ=255)
                                                                           agcgcTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGAGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGATCCTCttt                   >  2:538282/2‑100 (MQ=255)
                                                                            gcgcTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGTGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGATCCTCTTTg                  <  1:1002250/100‑1 (MQ=255)
                                                                            gcgcTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGATCCTCttt                   >  2:1587942/1‑99 (MQ=255)
                                                                            gcgcTGATTAATCCTATGGGGATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGATCCTCTTTg                  <  1:1002897/100‑1 (MQ=255)
                                                                             cgcTGATTAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGATCCTCTTTGc                 <  2:2225492/100‑1 (MQ=255)
                                                                                   ttAATCCTATGGTTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGATCCTCTTTGCGGTAct           >  2:3002974/1‑100 (MQ=255)
                                                                                    tAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGATCCTCTTTGCGGTActg          >  1:1814540/1‑100 (MQ=255)
                                                                                    tAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGATCCTCTTTGCGGTActg          <  2:1043852/100‑1 (MQ=255)
                                                                                    tAATCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGATCCTCTTTGCGGTAct           >  2:2236386/1‑99 (MQ=255)
                                                                                       tCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCATACCGGTATCCTGATCCTCTTTGCGGTActgttc       <  2:1815434/100‑1 (MQ=255)
                                                                                       tCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGATCCTCTTTGCGGTActgttc       <  2:2469448/100‑1 (MQ=255)
                                                                                       tCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGATCCTCTTTGCGGTActgttc       <  2:1323579/100‑1 (MQ=255)
                                                                                       tCCTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGATCCTCTTTGCGGTActgttc       >  1:1861326/1‑100 (MQ=255)
                                                                                        ccTATGGGTATGGTTGGTCAGGATTTGTGGGTTTACTTTATGCAGCTTCTGCCTTATATGCTGCAAAGCGGTATCATGCTCCTCTTTGCGGTActattc       >  2:848537/1‑99 (MQ=255)
                                                                                        ccTATGGGTATGGTTGGTCAGGATCTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCTGTATCCTGTTCCTCTTTGCGGTActgttc       <  2:1000600/99‑1 (MQ=255)
                                                                                         cTATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCTGTATCCTGATCCTCTTTGCGGTActgttc       <  2:1274428/98‑1 (MQ=255)
                                                                                          tATGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGATCCTCTTTGCGGTActgttctgt    <  1:2463748/100‑1 (MQ=255)
                                                                                          naTGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGATCCTCTTTGCGGTActgttctgt    >  1:2138984/2‑100 (MQ=255)
                                                                                           aTGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCATCTTCTGCCTTATATGCTGCAAACCAGTATCCTGATCCTCTTTGCGGTActgttctgtt   >  1:1176082/1‑100 (MQ=255)
                                                                                            tGGGTATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGATCCTCTTTGCGGTACTGTTCTGTTg  <  1:1817855/100‑1 (MQ=255)
                                                                                               |                                                                                                
TCTGTTGATCCTGACGCTCGCGCAAACTGTCGTCGCGACCTGGTATATGAAGACCATTCTTCCTTATCAGGGTTGGGCGCTGATTAATCCTATGGATATGGTTGGTCAGGATTTGTGGGTTTCCTTTATGCAGCTTCTGCCTTATATGCTGCAAACCGGTATCCTGATCCTCTTTGCGGTACTGTTCTGTTG  >  NC_000913/1111282‑1111473

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