Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,787,375 T→A 7.3% pseudogene (979/2253 nt) ygaQ → uncharacterized protein; putative enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,787,3750TA 7.3% 231.8 / 5.9 67pseudogene (979/2253 nt)ygaQuncharacterized protein; putative enzyme
Reads supporting (aligned to +/- strand):  ref base T (42/20);  new base A (2/3);  total (44/23)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGTACCATCCATAATCCAATGGGAGATTTGCAGGATTTTATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAA  >  NC_000913/2787226‑2787524
                                                                                                                                                      |                                                                                                                                                     
tttGTACCATCCATAATCCAATGGGAGATTTGCAGGATTTTATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGACTTGGGGGTAGTGGTAAAAATGACCTGGATTATCCct                                                                                                                                                       >  2:1401482/1‑150 (MQ=255)
tttGTACCATCCATAATCCAATGGGAGATTTGCAGGATTTTATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCct                                                                                                                                                       >  2:42100/1‑150 (MQ=255)
tttGTACCATCCATAATCCAATGGGAGATTTGCAGGATTTTATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCct                                                                                                                                                       >  2:1719953/1‑150 (MQ=255)
  tGTACCATCCATAATCCAATGGGAGATTTGCAGGATTTTATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCctct                                                                                                                                                     <  1:1724276/150‑1 (MQ=255)
                     tGGGAGATTTGCAGGATTTTATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGa                                                                                                                                  <  2:256843/150‑1 (MQ=255)
                     tGGGAGATTTGCAGGATTTTATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGGGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGa                                                                                                                                  <  1:1285776/150‑1 (MQ=255)
                       ggAGATTTGCAGGATTTTATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATc                                                                                                                                >  1:1371084/1‑150 (MQ=255)
                       ggAGATTTGCAGGATTTTATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATc                                                                                                                                >  1:1352722/1‑150 (MQ=255)
                       ggAGATTTGCAGGATTTTATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATc                                                                                                                                >  1:1371097/1‑150 (MQ=255)
                       ggAGATTTGCAGGATTTTATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATc                                                                                                                                >  1:412558/1‑150 (MQ=255)
                       ggAGATTTGCAGGATTTTATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATc                                                                                                                                >  2:1640229/1‑150 (MQ=255)
                        gagaTTTGCAGGATTTTATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCa                                                                                                                               >  1:888258/1‑150 (MQ=255)
                         agaTTTGCAGGATTTTATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGACTTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAg                                                                                                                              >  1:747962/1‑150 (MQ=255)
                         agaTTTGCAGGATTTTATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAg                                                                                                                              >  1:1108554/1‑150 (MQ=255)
                          gaTTTGCAGGATTTTATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGt                                                                                                                             <  1:471630/150‑1 (MQ=255)
                             ttGCAGGTTTTTATTGCATTGGTTAAAGATGCTAAAATATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTaaa                                                                                                                          <  1:773303/150‑1 (MQ=255)
                               gCAGGATTTTATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGa                                                                                                                        >  2:446728/1‑150 (MQ=255)
                                cAGGATTTTATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGat                                                                                                                       >  1:1001219/1‑150 (MQ=255)
                                    aTTTTATTGCATTGGTTAAAGATGCTAAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGAtata                                                                                                                    >  1:1292937/1‑150 (MQ=255)
                                      tttATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGaa                                                                                                                 <  1:241932/150‑1 (MQ=255)
                                      tttATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGaa                                                                                                                 >  1:1286120/1‑150 (MQ=255)
                                         aTTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAggt                                                                                                              <  1:596477/150‑1 (MQ=255)
                                            gCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAggtggt                                                                                                           >  2:381772/1‑150 (MQ=255)
                                              aTTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTAt                                                                                                         >  1:1483686/1‑150 (MQ=255)
                                                   ttAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTg                                                                                                    >  1:1266945/1‑150 (MQ=255)
                                                       aGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTa                                                                                                >  1:1276872/1‑150 (MQ=255)
                                                           gCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAg                                                                                            >  2:255316/1‑150 (MQ=255)
                                                                aaaaaTTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAg                                                                                        >  1:1330075/1‑150 (MQ=255)
                                                                  aaaTTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGt                                                                                      <  2:1292937/150‑1 (MQ=255)
                                                                       tGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGt                                                                                 >  2:1522755/1‑150 (MQ=255)
                                                                           aTCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTg                                                                             >  2:1664692/1‑150 (MQ=255)
                                                                            tCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGa                                                                            >  2:1535800/1‑150 (MQ=255)
                                                                            tCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGa                                                                            <  1:1166738/150‑1 (MQ=255)
                                                                              gATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGAtt                                                                          >  2:722937/1‑150 (MQ=255)
                                                                              gATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGAtt                                                                          >  1:135549/1‑150 (MQ=255)
                                                                              gATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGAtt                                                                          >  1:135393/1‑150 (MQ=255)
                                                                              gATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGAtt                                                                          >  1:770284/1‑150 (MQ=255)
                                                                                tatCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTcc                                                                        <  1:434243/150‑1 (MQ=255)
                                                                                tatCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTcc                                                                        >  1:241670/1‑150 (MQ=255)
                                                                                 atCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCACTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCca                                                                       <  1:890852/150‑1 (MQ=255)
                                                                                     ttCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCcattca                                                                   >  2:1240991/1‑150 (MQ=255)
                                                                                          gATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTaaa                                                              <  1:1133713/150‑1 (MQ=255)
                                                                                          gATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTaaa                                                              <  1:1249065/150‑1 (MQ=255)
                                                                                          gATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTaaa                                                              <  2:888258/150‑1 (MQ=255)
                                                                                             tataCCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGTCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGAt                                                           >  2:215195/1‑150 (MQ=255)
                                                                                             tataCCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGAt                                                           >  2:1668157/1‑150 (MQ=255)
                                                                                               taCCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                         <  2:1001219/150‑1 (MQ=255)
                                                                                                  cTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCACTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCcat                                                                      <  2:117970/134‑1 (MQ=255)
                                                                                                  cTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCACTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCcat                                                                      >  1:117918/1‑134 (MQ=255)
                                                                                                  cTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCACTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCcat                                                                      <  2:117918/134‑1 (MQ=255)
                                                                                                  cTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCACTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCcat                                                                      >  1:117970/1‑134 (MQ=255)
                                                                                                        cTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACaaa                                                >  2:1435787/1‑150 (MQ=255)
                                                                                                        cTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACaaa                                                <  1:22493/150‑1 (MQ=255)
                                                                                                                 aaTTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTaa                                       >  2:1359663/1‑150 (MQ=255)
                                                                                                                     gggggTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATTGAAGGTGGTATCCCTGGCTATTAGAAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACaa                                   <  2:356610/150‑1 (MQ=255)
                                                                                                                      ggggTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACaaa                                  >  1:875194/1‑150 (MQ=255)
                                                                                                                         gTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatcc                               >  1:755149/1‑150 (MQ=255)
                                                                                                                              ggTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                          >  2:769674/1‑150 (MQ=255)
                                                                                                                               gTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCag                         <  2:1214966/150‑1 (MQ=255)
                                                                                                                               gTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCag                         >  2:745740/1‑150 (MQ=255)
                                                                                                                               gTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCag                         >  2:242295/1‑150 (MQ=255)
                                                                                                                                 aaaaaTGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCagaa                       <  2:1084101/150‑1 (MQ=255)
                                                                                                                                     aTGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGata                   >  2:594723/1‑150 (MQ=255)
                                                                                                                                          cTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAt              >  2:1307311/1‑150 (MQ=255)
                                                                                                                                           tGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATGAAAGATCCAGTAACAAAAGAACGGAAACAAATCCATCCAGAAGATATACAtc             >  2:1719657/1‑150 (MQ=255)
                                                                                                                                               ttATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCAc         >  1:1038775/1‑150 (MQ=255)
                                                                                                                                                 aTCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTg       <  1:1055159/150‑1 (MQ=255)
                                                                                                                                                   cccTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCa     >  2:945257/1‑150 (MQ=255)
                                                                                                                                                   cccTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCa     >  1:279328/1‑150 (MQ=255)
                                                                                                                                                    ccTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCaa    <  2:165360/150‑1 (MQ=255)
                                                                                                                                                    ccTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCaa    <  2:644116/150‑1 (MQ=255)
                                                                                                                                                     ctctGCTGATATATGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCaaa   <  1:97161/150‑1 (MQ=255)
                                                                                                                                                     ctctGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCaaa   <  1:626831/150‑1 (MQ=255)
                                                                                                                                                      tctGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCaaaa  <  2:550773/150‑1 (MQ=255)
                                                                                                                                                      tctGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCaaaa  <  2:1517400/150‑1 (MQ=255)
                                                                                                                                                      tctGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCaaaa  <  2:136509/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TTTGTACCATCCATAATCCAATGGGAGATTTGCAGGATTTTATTGCATTGGTTAAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAA  >  NC_000913/2787226‑2787524

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: