Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,856,683 C→T 8.3% A90V (GCG→GTG)  yidI → inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,856,6830CT 8.3% 228.2 / 9.4 72A90V (GCG→GTG) yidIinner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (35/31);  new base T (3/3);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

TGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCT  >  NC_000913/3856534‑3856832
                                                                                                                                                     |                                                                                                                                                     
tGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATtgctgc                                                                                                                                                       >  2:799009/1‑150 (MQ=255)
   aaGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCt                                                                                                                                                    >  1:259558/1‑150 (MQ=255)
      cATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACt                                                                                                                                                 >  2:674088/1‑150 (MQ=255)
      cATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACt                                                                                                                                                 >  2:1372764/1‑150 (MQ=255)
                  cAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTAt                                                                                                                                     >  2:985787/1‑150 (MQ=255)
                    aaaCGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTg                                                                                                                                   >  1:952502/1‑150 (MQ=255)
                     aaCGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTgg                                                                                                                                  >  2:84548/1‑150 (MQ=255)
                      aCGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGTCATTGCTGCGCTACTTTGCGCATTTATTGGc                                                                                                                                 >  1:1366752/1‑150 (MQ=255)
                      aCGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGTCATTGCTGCGCTACGTGGCGCAGTTATTGGc                                                                                                                                 >  1:459783/1‑150 (MQ=255)
                      aCGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGc                                                                                                                                 <  1:849562/150‑1 (MQ=255)
                       cGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCg                                                                                                                                <  1:1356369/150‑1 (MQ=255)
                       cGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCg                                                                                                                                <  2:1402426/150‑1 (MQ=255)
                        ggCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCgg                                                                                                                               <  2:530673/150‑1 (MQ=255)
                         gCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGt                                                                                                                              >  2:588615/1‑150 (MQ=255)
                           tGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTAt                                                                                                                            >  2:396305/1‑150 (MQ=255)
                               aTTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGt                                                                                                                        >  2:1253298/1‑150 (MQ=255)
                                ttAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTa                                                                                                                       >  1:554918/1‑150 (MQ=255)
                                ttAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTa                                                                                                                       >  1:554945/1‑150 (MQ=255)
                                 tAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTaa                                                                                                                      <  2:1118448/150‑1 (MQ=255)
                                  aaaaGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAg                                                                                                                     <  1:1258397/150‑1 (MQ=255)
                                  aaaaGAGGTCTCCTCGCCGGAATAAAAGGCGAAAAAAGTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAg                                                                                                                     <  1:716303/150‑1 (MQ=255)
                                    aaGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGa                                                                                                                   >  1:1365021/1‑150 (MQ=255)
                                    aaGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGa                                                                                                                   <  1:524918/150‑1 (MQ=255)
                                        ggTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAAt                                                                                                               >  1:148209/1‑150 (MQ=255)
                                        ggTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAAt                                                                                                               >  1:857689/1‑150 (MQ=255)
                                         gTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATc                                                                                                              >  1:1324622/1‑150 (MQ=255)
                                            ttcgCGCCGGAAGAAAAGGCGACAAAATGACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCCTTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCt                                                                                                           <  1:1333412/150‑1 (MQ=255)
                                               tcGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTggg                                                                                                        >  1:66659/1‑150 (MQ=255)
                                                cGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTgggg                                                                                                       <  2:993696/150‑1 (MQ=255)
                                                 gCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGc                                                                                                      >  2:998908/1‑150 (MQ=255)
                                                  ccGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCa                                                                                                     >  2:279193/1‑150 (MQ=255)
                                                         aaaaGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCgcgg                                                                                              >  1:1268174/1‑150 (MQ=255)
                                                          aaaGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCgcggc                                                                                             <  2:393968/150‑1 (MQ=255)
                                                                gAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTg                                                                                       >  1:1160365/1‑150 (MQ=255)
                                                                 aaaaaaTTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTgg                                                                                      >  2:1201147/1‑150 (MQ=255)
                                                                   aaaaTTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGtt                                                                                    <  2:1341862/150‑1 (MQ=255)
                                                                      aTTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGttttt                                                                                 <  2:1108921/150‑1 (MQ=255)
                                                                      aTTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGttttt                                                                                 <  2:567566/150‑1 (MQ=255)
                                                                           caccGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTggcgg                                                                            <  1:9241/150‑1 (MQ=255)
                                                                           caccGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTggcgg                                                                            <  1:784222/150‑1 (MQ=255)
                                                                             ccGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGgcggcg                                                                          >  1:703441/1‑150 (MQ=255)
                                                                             ccGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGgcggcg                                                                          >  2:1286419/1‑150 (MQ=255)
                                                                                    gaAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCAGTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCt                                                                   <  2:1260883/150‑1 (MQ=255)
                                                                                      aaaaaGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTgg                                                                 <  2:568690/150‑1 (MQ=255)
                                                                                          agagGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCttt                                                             <  1:343083/150‑1 (MQ=255)
                                                                                              gACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGTGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCAt                                                         >  1:494918/1‑150 (MQ=255)
                                                                                              gACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGTGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCAt                                                         >  1:495372/1‑150 (MQ=255)
                                                                                              gACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGTGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCAt                                                         >  1:494935/1‑150 (MQ=255)
                                                                                               aCGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATa                                                        >  1:1361777/1‑150 (MQ=255)
                                                                                                  tAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCt                                                     <  1:1291365/150‑1 (MQ=255)
                                                                                                  tAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCt                                                     <  1:908983/150‑1 (MQ=255)
                                                                                                          gATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCttttt                                             <  1:998908/150‑1 (MQ=255)
                                                                                                           aTAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGTGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTg                                            <  2:495372/150‑1 (MQ=255)
                                                                                                           aTAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGTGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTg                                            <  2:494918/150‑1 (MQ=255)
                                                                                                           aTAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGTGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTg                                            <  2:494935/150‑1 (MQ=255)
                                                                                                           aTAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTg                                            <  1:467484/150‑1 (MQ=255)
                                                                                                           aTAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTg                                            <  1:467498/150‑1 (MQ=255)
                                                                                                            tAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTgg                                           >  1:1350620/1‑150 (MQ=255)
                                                                                                              aaaTCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTa                                         <  2:47605/150‑1 (MQ=255)
                                                                                                               aaTCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTAt                                        >  2:1252544/1‑150 (MQ=255)
                                                                                                               aaTCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTAt                                        <  1:279193/150‑1 (MQ=255)
                                                                                                               aaTCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTAt                                        <  2:4640/150‑1 (MQ=255)
                                                                                                                aTCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATc                                       <  1:1020284/150‑1 (MQ=255)
                                                                                                                   cTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGc                                    >  1:821139/1‑150 (MQ=255)
                                                                                                                    tcaatcaaAGCGGTATCGCACTGGGCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCa                                   >  2:560142/1‑150 (MQ=255)
                                                                                                                        tcaaAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCgt                               >  1:778130/1‑150 (MQ=255)
                                                                                                                          aaaGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCgtgt                             >  1:920135/1‑150 (MQ=255)
                                                                                                                          aaaGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCgtgt                             >  1:829627/1‑150 (MQ=255)
                                                                                                                          aaaGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCgtgt                             >  2:810286/1‑150 (MQ=255)
                                                                                                                             gCGGTATCGCACTGGCCAGTGCTGCGCTACTTTGCGCATTTATGGGCGGTATGCGTAAGGAACATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTgca                          <  2:507445/150‑1 (MQ=255)
                                                                                                                                         tGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGa              >  1:1065798/1‑150 (MQ=255)
                                                                                                                                            ccATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGatat           <  1:398198/150‑1 (MQ=255)
                                                                                                                                            ccATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGatat           <  1:1189261/150‑1 (MQ=255)
                                                                                                                                                tgctgcGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCt       <  1:1375780/150‑1 (MQ=255)
                                                                                                                                                tgctgcGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCt       <  1:956113/150‑1 (MQ=255)
                                                                                                                                                tgctgcGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCt       <  1:645788/150‑1 (MQ=255)
                                                                                                                                                    gcgcTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATc   >  1:675564/1‑150 (MQ=255)
                                                                                                                                                     cgcTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCt  >  2:643904/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGATAAAATCCTCAATCAAAGCGGTATCGCACTGGCCATTGCTGCGCTACTTTGCGCATTTATTGGCGGTATGCGTAAGGAAAATCGCTGGGGCATTCGCGGCGCGCTGGTTTTTGGCGGCGGGACGCTGGCTTTCCATACCTTGCTTTTTGGTATCGGCATCGTGTGCAGCATTTTGCTGATATTCCTTATCT  >  NC_000913/3856534‑3856832

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