Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,765,794 (C)7→8 intergenic (‑18/‑124) yfjM ← / → rnlA CP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, RNase LS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,765,7871.C90.5% 202.8 / 14.0 63intergenic (‑11/‑131)yfjM/rnlACP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, RNase LS
Reads supporting (aligned to +/- strand):  ref base . (2/4);  new base C (26/31);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

CTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGA  >  NC_000913/2765654‑2765916
                                                                                                                                      |                                                                                                                                  
ctctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccc                                                                                                                              <  2:28511/139‑1 (MQ=255)
 tctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccccc                                                                                                                             >  1:924290/1‑139 (MQ=255)
       cGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccccccggta                                                                                                                       <  2:322665/139‑6 (MQ=255)
        gTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGa                                                                                                                      <  2:1175543/139‑1 (MQ=255)
           cTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccccccccggtagac                                                                                                                   >  2:1057552/1‑130 (MQ=255)
              tttCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATc                                                                                                                 <  2:736151/139‑1 (MQ=255)
                tCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGGATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCt                                                                                                               <  2:592627/139‑1 (MQ=255)
                  gACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGcc                                                                                                            >  1:422672/1‑139 (MQ=255)
                    cGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGccc                                                                                                           >  2:92396/1‑139 (MQ=255)
                      tAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccc                                                                                                                               >  2:114553/1‑116 (MQ=255)
                      tAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccc                                                                                                                               <  1:114553/116‑1 (MQ=255)
                                tCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                                                                >  2:734255/1‑139 (MQ=255)
                                tCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGt                                                                                               <  1:92396/139‑1 (MQ=255)
                                 cTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGt                                                                                               <  1:734255/139‑1 (MQ=255)
                                            gCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTCAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATGTAATGGACTGCCCCCCCCGGTAGCCGATCCTGCCATATAGTTGGAGCGCCGCCTATTCa                                                                                   >  2:159119/1‑139 (MQ=255)
                                                          aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTa                                                                     >  1:769867/1‑139 (MQ=255)
                                                          aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                                                    >  2:324914/1‑139 (MQ=255)
                                                            tATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                                                    <  1:1057552/139‑1 (MQ=255)
                                                            tATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                                                   >  2:115762/1‑139 (MQ=255)
                                                             aTCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCa                                                                  <  1:115762/139‑1 (MQ=255)
                                                               cGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGa                                                                                                 <  1:910019/106‑1 (MQ=255)
                                                               cGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGa                                                                                                 >  2:910019/1‑106 (MQ=255)
                                                               cGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cc                                                                                                                                 >  1:1175735/1‑73 (MQ=255)
                                                               cGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cc                                                                                                                                 <  2:1175735/73‑1 (MQ=255)
                                                               cGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                                                            >  2:1004563/1‑110 (MQ=255)
                                                               cGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                                                            <  1:1004563/110‑1 (MQ=255)
                                                                 aGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATc                                                              >  1:747073/1‑139 (MQ=255)
                                                                 aGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATc                                                              >  1:180493/1‑139 (MQ=255)
                                                                  gTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCaa                                                            <  1:324914/139‑1 (MQ=255)
                                                                    tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAAt                                                           <  1:1220841/139‑1 (MQ=255)
                                                                    tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                          <  2:422672/139‑1 (MQ=255)
                                                                     gctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccccc                                                                                                                            <  1:474569/72‑1 (MQ=255)
                                                                     gctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccccc                                                                                                                            >  2:474569/1‑72 (MQ=255)
                                                                      ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTg                                                                                                              >  1:266268/1‑86 (MQ=255)
                                                                      ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTg                                                                                                              <  2:266268/86‑1 (MQ=255)
                                                                      ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCa                                                                  >  1:978449/1‑130 (MQ=255)
                                                                      ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCa                                                                  <  2:978449/130‑1 (MQ=255)
                                                                      ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                                                         >  1:1027279/1‑139 (MQ=255)
                                                                         ctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACt                                                      >  2:1106291/1‑139 (MQ=255)
                                                                         ctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccccccc                                                                                                                           >  1:1137556/1‑68 (MQ=255)
                                                                         ctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccccccc                                                                                                                           <  2:1137556/69‑2 (MQ=255)
                                                                              tGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                <  1:537246/139‑1 (MQ=255)
                                                                               gagaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                <  1:90520/139‑1 (MQ=255)
                                                                               gagaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                <  1:159119/139‑1 (MQ=255)
                                                                               gagaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                <  2:1104132/139‑1 (MQ=255)
                                                                                agaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg                                               <  2:924290/139‑1 (MQ=255)
                                                                                 gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                          >  2:1050825/1‑127 (MQ=255)
                                                                                 gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                          <  1:1050825/127‑1 (MQ=255)
                                                                                  aTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                             <  2:769867/139‑1 (MQ=255)
                                                                                    cAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                           <  2:747073/139‑1 (MQ=255)
                                                                                     aGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtat                                          <  1:1106291/139‑1 (MQ=255)
                                                                                        aaGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtattat                                       <  1:383126/139‑1 (MQ=255)
                                                                                            ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCa                                                                  >  1:592302/1‑109 (MQ=255)
                                                                                            ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCa                                                                  <  2:592302/109‑1 (MQ=255)
                                                                                            ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGccc                                                                                                           <  1:1002517/67‑1 (MQ=255)
                                                                                            ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGccc                                                                                                           >  2:1002517/1‑67 (MQ=255)
                                                                                            ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccc                                                                                                                               <  2:151453/46‑1 (MQ=255)
                                                                                            ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccc                                                                                                                               >  1:151453/1‑46 (MQ=255)
                                                                                             gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTag                                  >  2:1052830/1‑139 (MQ=255)
                                                                                               tACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGTAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTagag                                >  1:119341/1‑139 (MQ=255)
                                                                                                aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATc                                                                                                                 <  1:876837/57‑1 (MQ=255)
                                                                                                aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATc                                                                                                                 >  2:876837/1‑57 (MQ=255)
                                                                                                aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                           <  1:1103356/127‑1 (MQ=255)
                                                                                                aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                           >  2:1103356/1‑127 (MQ=255)
                                                                                                aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGt                               <  2:1027279/139‑1 (MQ=255)
                                                                                                     aaCGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccc                                                                                                                              >  2:663867/1‑38 (MQ=255)
                                                                                                     aaCGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccc                                                                                                                              <  1:663867/38‑1 (MQ=255)
                                                                                                      aCGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCa                                                                                   >  1:1113450/1‑81 (MQ=255)
                                                                                                      aCGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCa                                                                                   <  2:1113450/81‑1 (MQ=255)
                                                                                                                 catcCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                  >  2:1023277/1‑103 (MQ=255)
                                                                                                                 catcCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                  <  1:1023277/103‑1 (MQ=255)
                                                                                                                  atcCCCCTTTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACCCTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTAt             <  1:192565/139‑1 (MQ=255)
                                                                                                                  atcCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTAt             <  2:180493/139‑1 (MQ=255)
                                                                                                                        cATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCa                                                                  >  1:322888/1‑80 (MQ=255)
                                                                                                                        cATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCa                                                                  <  2:322888/80‑1 (MQ=255)
                                                                                                                           tttGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAg    >  1:589779/1‑139 (MQ=255)
                                                                                                                             tgtgGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGa  >  2:533804/1‑139 (MQ=255)
                                                                                                                                ggACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcccc                                                                                          <  2:835762/48‑1 (MQ=255)
                                                                                                                                ggACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcccc                                                                                          >  1:835762/1‑48 (MQ=255)
                                                                                                                                      |                                                                                                                                  
CTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGA  >  NC_000913/2765654‑2765916

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: